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Bemerkungen zur Stabilität eines physikalischen Pendels 


Zweite Mitteilung 
Von Wolfhart Haacke in Braunschweig 


In dieser zweiten Mitteilung werden die Stabilisierungsphänomene des physikalischen Pendels bei be- 
Tiebiger periodischer ebenen Bewegung des Aufhängepunktes untersucht, falls die Schwerkraftrichtung in 
dieser Bewegungsebene liegt. Bei Benutzung der Theorie der Mathieuschen Differentialgleichung, besonders 
der Struttschen Parameterkarte, kann systematisch und mit geringem Rechenaufwand entschieden werden, 
welche stabilen und labilen Lagen das Pendel einnehmen kann. Ausführlich werden die harmonische ellip- 
tische Oszillation und ihre Spezialfälle diskutiert. 


In.this second communication, ihe phenomena of stabilizing a physical pendulum are investigated, if 
the point of suspension is moving periodically, in an arbitrary manner, in a plane that contains the direction 
of gravity. The theory of the Mathieu differential equation, especially Strutt’s parameter chart permit to 

} discern, systematically and with little computation work, the stable and unstable positions of the pendulum. 
The harmonic elliptical oscillation and its special cases are discussed in detail. 


Dans cette deuwieme communication on ewamine les phenomenes de stabilisation d’un pendule physique, 
dont le point de suspension subit un mouwvement plan periodiquwe arbitraire, dans le cas, que la direction 
de la gravitation est situee dans le plan de ce mouvement. En employant la theorie de l’&qualion differentielle 
deMathieu et particulierement la carte parametrigue de Struit, on peut determiner systömatiquement et 
sans grande peine calculatoire les positions stables et labiles du pendule. L’oscillation harmonique elliptique 
et ses cas particuliers sont discutes en detail. 

B 3T0M BTOPoM cooÖmeHnun HECHeAyIOTCA ABICHUA CTAÖMINSBAUUN duaNyecKoroO MAATHUKA 
Hpu J11000M IIEPHONHYEeCKOM INIOCKOM ABH3KCHHUH TOUYKH TIONBeCA, ECIH HAIPaBIeHNE CHJIEL 
TAIKeCTH JIe3KHT B 3TOU INIOSKOCTU ABuskenns. Menonssya Teopum AudepenumansHoro 
ypasuenmı Marse, B 0C0OÖeHHOCTH KapTy mapameTpoB ÜTPYTTA, MO3KHO CHCTEMATHYeCKH 
H NYTeM IPOCTBIX BbIymcHeHna PeIIHTB BONPOC, KAakMe YCTOÄyuBble U HEYCTOÜYHBBIE 
HONOSKEHHUA MOTYT ÖBITb 3AHATEI MASITHUROM. 1lo1po6H0o uccnenyeTcH raPMOHHYECKO® FZILIMIITH- 
yJecKoe KOJIeÖaHHe H ETO YACTHEIe CIYYan. 


1. Einführung . 


In einer früheren Untersuchung unter dem gleichen Titel wurden die Stabilisierungs- 
phänomene eines physikalischen Pendels bei vertikaler harmonischer Erschütterung des Auf- 
hängepunktes untersucht [1]?). Nun soll die Fragestellung verallgemeinert und beliebige ebene 
periodische Oszillationen des Aufhängepunktes zugelassen werden. Allerdings wird verlangt, daß 
der Vektor der Schwerkraft in der Oszillationsebene liegt. Wird durch die Bewegung des Auf- 

hängepunktes und die Schwerkraftrichtung eine Ebene eindeutig bestimmt, so stellt sich aus 
Symmetriegründen ein räumliches Pendel in diese Ebene ein und schwingt als ebenes Pendel. 
Diese Eindeutigkeit besteht in jedem Fall außer dem der vertikalen Erschütterung. Unter- 
suchungen über dieses Problem sind schon von Hirsch [2], Erdölyi [3] und Klotter [4] durch- 
geführt worden. In der vorliegenden Arbeit wird mit Hilfe der Theorie der Mathieuschen Diffe- 
rentialgleichung das Pendelproblem neben einigen neuen Ergebnissen systematischer und mit 
geringerem Rechenaufwand gelöst. Dadurch kann man auf einen Teil der früher benötigten Ver- 
nachlässigungen verzichten und die Auswirkungen der übrigen besser beurteilen. Teilweise ist 
in den früheren Arbeiten von einer periodischen Bewegung des Aufhängepunktes gesprochen und 
ihre Frequenz berechnet worden. ‘Aus den folgenden Betrachtungen wird sich ergeben, daß eine 
fastperiodische, aber nicht periodische Bewegung Voraussetzung für eine Stabilisierung ist. 
Bezüglich der Sätze aus der Theorie der Mathieuschen Differentialgleichung und die Aus- 
wertungsmethode der Struttschen Karte schließen wir unmittelbar an [1] an. 


er 2. Die Bewegungsgleichung 


Wir betrachten ein physikalisches Pendel mit der Masse m, das um einen Aufhängepunkt A 
_ in einer @,y-Ebene schwingen kann. Dieser Aufhängepunkt werde periodisch in dieser Ebene er- 
 schüttert. Dabei seien die Erschütterungsamplituden in der z- und in der y-Richtung klein gegen 
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den Abstand s des Aufhängepunktes A vom Schwerpunkt S des Pendels. Die Schwerkraft wirke 


in negativer y-Richtung, ® sei der Winkel der positiven y-Richtung mit AS. Sind a, y die Ko- a 
ordinaten von A und &, und y, diejenigen von 8, so gilt 


% =20—ssin®, 


1, =y+4scos ® 
Mit der kinetischen Energie 


1 { pe 
eds D® — ms D (x cos B+5sind)+ (+3) 


und der potentiellen Energie 
U=mg(y-+scos ®) 
ergibt sich die Lagrangesche Gleichung 
9, dB —ms[ücos d+ (Y+gsind]=0..(l). 


Dabei ist ©, das Trägheitsmoment um A. Wir setzen 
nun | 


A, % Y Dr u > 

l= —.—, = —,, — 

Aa ms’ S 2} 7 wi R ‚An®n?] oe 

undu=2zrnt. Die EURE nach u werden durch Striche bezeichnet. Dann ergibt sich aus (1) % 
— fer cos DA mE’ RSin Dj Or (2) 

Wie in [1] wollen wir nach K Existenz kleiner Sn a um ee fragen, die durch or 


am! [ol | 


i>» 


definiert werden. Setzen wir D=9--a, so können wir rl<Z n_ annelituen. Dann läßt ac 0) 
linearisieren und wir erhalten 


+ [&"sina— (m + h)cosa]p=£"cosa-+ (N h)sin a ne 
Betrachten wir eine allgemeine elliptische Oszillation des Aufhängepunktes, so ist 
&=asinwcosu-+ bcos wsin u, 
N= — 0005 cos W+ bsin wsin u 


O<bsa< und DE). 


Über N bzw. n setzen wir zunächst noch ut voraus. Dann erhalten wir aus (3) 
9" +A+YycsW)p=co-+ BIERE (ua) u.2 Lan: 

Dabei bedeuten Fe 

A=—hcoso, v1 ra ar 


Dhsine. rer und 
; e? sin 2(x — ®) 
 Y4arb?-+ (a? — b2)? sin? 2 (@—0o) ©) 


co — 


3. Die Stabilitätsbedingung 


(4) ist ei eine inhomogene Mathieusche Ditferentialgleichung. Die 
Gl. (4), hat nach dem Flo ae Ana Gb & in: folgende G 


N 
RT Ye ‚ 
Dr Dr 
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Dabei sind a,unda, ee tlanskonstanten, FW) =flu+2r) und ureell, doch u= 5 En —o,1, ‚3% 


da wir nur stabile Bewegungen betrachten, bei denen o klein bleibt. Bei dieser Einschränkung für 
w existiert eine in 2 periodische Partikularlösung von (4) ([6] S. 145). Doch kleine Schwingungen z 


_ um eine Mittellage ergeben sich höchstens dann, wenn il 
| h o— lim — 1 (aan dr oder rn, [ot de) IR 
t>o m” 

erfüllt ist. Es gilt j 4 
Ss | 
lim 2 one) RER RR BE re A ER 1]. u 

u>o U v e 
Da die Partikularlösung 9, in 2 periodisch ist, folgt aus 2 
u 2 

1 | % 
lim 1 [ec dıe= is im | 910 (z)ar — 0 er 

un U d B3: 


“0, wenn man p,(w) in ihre Fourierreihe 


' + ö Ban 3% 
Bot: r a EIER TEE NN HN C (6) "a 
te ee Er 
B 
” entwickelt. Dies ist die notwendige Bedingung für das Auftreten kleiner Schwingungen. Aus 6,0, vs 
. werden die „möglichen Mittellagen“ & bestimmt. Dann geht man mit dem nun bekannten «a in - Ba. 
die homogene GI]. (4), ein und diskutiert die „Stabilität“ ihres Fundamentalsystems in Abhängig- Be 
keit von a, b, nund w an Hand der Struttschen Karte (Figur 2in [1]). Im allgemeinen sollen 73 
2 a, bund w vorgegeben und die Erschütterungsfrequenzen bestimmt werden, für die sich Sta- Bi: 
&  bilität ergibt?). u 
ee; " Setzt man (6) in (4) ein, so erhält man ein unendliches Gleichungssystem für die c,: BB 
r . \ - we; 4 
Be, ww! Ent. 
>. ati +,0ı= A Be . 
2 ‚ Bu 
Bor. . } z : >. 
Bi. HER were NS b- 
, EN a ade a 
w 6, ya 5a ») (y—ı ar. +1) = ne 0 für | > 1 ß $ y wi 
a - 
Die Koeffizientendeterminante C konvergiert dann nach der Kochschen Theorie [7]. B sei die we) 
_ Determinante,. die sich aus C Be wenn man die nullte ee durch die rechte Seite ersetzt: Se 
a ey N ö rn E: 
20-9) u : 
N Er | 
Br EN ra 13 0 = 
a 0 | 9 o 36 0 E 
re f > Dr ? HR 0 a) 
er y art wa | Ri 
z r 
Br 
a 
4 an die de Warn a die zu instabilen Punkten der Struttschen Karte führen, R We ” 
S. heißen re ns ’ f e- ; B;. 
x 23* All 
% I 


ist für |A|< 1 analytisch. Also bis auf O (y*): 


.. „Im allgemeinen wird man sich auf größere Erschütterungsfrequenzen höschränkene sodBß 
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Nach Koch konvergiert B ebenfalls und es ist c, = Falls (6) überhaupt existiert, folgt 


Da} j va 
aus c,—=0 auch B=0. Im Ausnahmefall liegen weder stabile noch labile Lagen vor. Es 
ist nun 


+0 
92 n 1 { 2 oycosx 
ee DI RE 


In der eckigen Klammer steht eine für A = 0 reguläre Funktion. y und o haben dieselbe Größen-. 
ordnung. Über Awollen wir nur voraussetzen, daß |A|<1 ist, d. h. n®1>25 cm/see”?, Bei einem 
Zylinderstabpendel von 25 cm etwa könnte n dann bis zu 1 Hz herunter Drau werden. Es 
gilt 
ae Be] Nat 3S 177086 
ee ee 


Also ist 
SE NE ENER 
„4-9 A— (+1)? 


vyv=— 


1 il ee 1 | | 
KR: = 
ik > ran, 1 2yACRH1) , say , el _ 2mcotayA ® 
| Werra j—» YA+r+1 Al = Yrd—aa), 5 
Dann ist no 
Di Yy ncot YA bee 1 OYCOSR | ga ” rcotm A K 1 
Bo he Te 


el ee] 1 1 mCotg a | 
h sin Al + 9 [a2 cos? (& Oz sin? («—w)] N eos U AR ERE 5 Thcosa(i-Hähcos 5 
_.e&sin2(& —w) 23 
TAU Ihe (M). 


Die Erfülltheit von (7) ist eine notwendige Bedingung dafür, daß i in der Richtung & Stabilität 
herrscht. (7) soll die Stabilitätsbedingung genannt werden. 


4. Diskussion der Stabilitätsbedingung 


man neben h,a? und 5? die Größen ha? und hb? vernachlässigen kann. Dann können wir auch Sy - 
1-+-hcos 21 setzen. Wenn man n nicht in’ dieser Form einschränken will, muß man uf 7) 
zurückgreifen. Es ist zu beachten, daß z als Funktion von A für || <0,8 dem Betrage nach 
unterhalb 1,5 bleibt. Es ist z(0) # — 0,3. Mit wachsendem 4 strebt z monoton gegen Null. 
Für A<0 oszilliert z zunächst, doch strebt z beiA > — 1 über alle Grenzen. Man wird im all- 
gemeinen (7) mit Hilfe eines Iterationsverfahrens lösen müssen. — Da nun z für |A| < 0,8 in der 
Größenordnung 1 bleibt, vereinfacht sich (7) bei den gegebenen Vernachlässigungen zu 


4hsna=esin2(a—®) ..... Re A Be 
Dah>0Ound ®> Oist, kann es wegen der Vorzeichen der trigonometrischen Funktionen - Br 
stensin folgenden Sektoren Mittellagen geben: Im Intervall 0 sası kommt nur w <a <5 at | 
undinz sa s2znur z <a <a-+ w und = +wsas2nin ns Be Be 


Im Falle zirkularer Oszillation ist 
E=acosu und 


Dann wird 
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Aus (8) erhalten wir 
sina=0. Füra=0 ist A= —h, y=a und für a=r ist A=h, y=a (8a). 


Da h und a klein sind, ist bei zirkularer Oszillation die natürliche Gleichgewichtslage stabil, da 
die Umgebung des Nullpunktes der Struttschen Karte im ersten Quadranten einen stabilen 
Be ? stabil [1]. Da- 
bei ist a der durch die reduzierte Pendellänge Z dividierte Radius des Oszillationskreises. Der 
Vergleich mit [1] ergibt, daß für große n und kleine a die Stabilisierung durch vertikale Er- 
schütterung und durch zirkulare Oszillation die- 
selben Phänomene hervorbringt. Schränken wir 
h nicht ein, so müssen wir auf (7) zurückgreifen. 
Wir erhalten wieder sina—=0. Die Diskussion 
für beliebige A < 1 wird wie bei der vertikalen 
Erschütterung durchgeführt [1]. Doch daneben 


2 
‘ tritt jetzt noch die Gleichung 1 —+ 5 z2=0 auf. 


Da |z| für || <0,8 unterhalb 1,5 bleibt, kann 
diese Gleichung bei unseren Voraussetzungen 
”) über « keine weiteren Wurzeln ergeben. — Nun 
N "wollen wir eine lineare Erschütterung in der 
Richtung w betrachten. Aus (8) erhalten wir jetzt 


\ Ahsin a= a2sin (a —ıw) . „u. 18). 


Daraus ergibt sich folgender Satz: Das Stabili- Bild 2 
sierungsphänomen eines durch elliptische Oszilla- 

tion des Aufhängepunktes erschütterten Pendels läßt sich immer auch durch eine lineare Er- 
schütterung erreichen, falls die Oszillationsfrequenz hinreichend groß und die Ellipsenachsen 
gegen die reduzierten Pendellängen hinreichend klein sind. Diese lineare Erschütterung muß in 
Richtung der großen Ellipsenachse erfolgen und als Amplitude die lineare Exzentrizität der 
Ellipse haben. — Der Spezialfall dieses Satzes für zirkulare Oszillation ist schon oben gesondert 
behandelt worden. Dieses Ergebnis hat bereits Erd&lyi ausgesprochen [3]. Doch scheint der 
hier benutzte Weg zu seiner Ableitung wesentlich kürzer zu sein. Die Abb. 2 hat für den Fall 
der linearen Erschütterung schon Klotter angegeben [4]. Die Unterschiede zwischen den Phäno- 
menen, die durch elliptische und durch lineare Oszillation erzeugt werden, ergeben sich erst, wenn 
wir n kleiner werden lassen, wie ein Vergleich von (7) mit der aus (7) durch 5 = 0 entstehenden 
Gleichung zeigt. Die Unterschiede sind also von höherer Ordnung klein. Eine graphische Dar- 


Bereich darstellt. Die obere Gleichgewichtslage & = 0 ist für n > = 


- p) ? 
stellung von (Sb) in einer - ‚&-Ebene mit w als Parameter hat Klotter angegeben [4]. Sie 


soll hier nicht wiederholt werden. Die Entscheidung, ob ein in dieser Ebene gefundener Parameter 
> - zu einer stabilen oder labilen Lage gehört, fällt man mit Hilfe der Struttschen Karte [1]. Bei 
linearer Erschütterung ist 


A=—hcosa. mie @lcos (aa) la 


Ist A > 0, also 5 SuS 2 und «& zugleich Wurzel von (8b), so ergibt die Struttsche Karte, 


x ; daß diese Parameter für hinreichend kleine h und «a sicher zu stabilen Lösungen Anlaß geben. 
y Für kleine a und große n gibt es also in den beiden unteren Quadranten keine labilen Lagen. 


_ Wie die Struttsche Karte zeigt, gilt dies nicht mehr für kleine n. Für n? > ER) ergibt 


sich sicher noch Stabilität, wie man mit Hilfe einer Abschätzung der Grenzkurve ce, zeigen 
kann [1]. 


> on | 2 
Kin ist 1 <.0, also == <a +5 ‚so ist Al < En für Stabilität notwendig. Setzt man in 


Rn 


” 4 ‚diese Bedingung (5b) und (8b) ein, so erhält man 
5, } 


N: 2h gi 
ee" rien 7 Iimaj< 05 ( —@) Re, a et (9). 
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Dies ist die hinreichende Bedingung für die Stabilität in den oberen Quadranten, nachdem a durch 
(8b) bestimmt ist. 


2 
daher ist in (7) kein Absolutzeichen bei cos (« — w) notwendig. Man könnte annehmen, daß 
die Gleichgewichtslagen von der Drehrichtung der elliptischen Oszillation abhängen. Doch 8 
zeigt, daß dies nicht der Fall ist. 


Setzt man = (0), so ergibt sich die in [1] ausführlich behandelte vertikale Erschütte- 
rung. Ergänzend dazu soll nur noch gesagt werden, daß sich aus (3b) dann neben sin «—=0 noch 


Nach Abb. 2 ist dann nur = sa —wo<s—o und 0<sa—w ERLEN möglich; 


203. — S ergibt. Doch (9) zeigt, daß es sich dabei um zwei zur Vertikalen symmetrische labile 


Lagen handelt. Die Parameterwerte liegen auf einer Grenzkurve der Struttschen Karte, die 
ja noch zum labilen Bereich gehört. Mit fallenden n schieben sich diese Lagen immer näher an 
die stabile Lage «= 0; beim Zusammenfallen (kritische Frequenz) dieser beiden Richtungen mit 
der vertikalen stabilen Lage wird diese nun labil. | 


Bei horizontaler Erschütterung (0-2) ist 


2 
sino@hA+acosa)=0 ....:... . (8c). ® 
Ist «= 0, so finden wir = — hl, y=0;istta=m,soisti=h Hk) y=0. (4), ist jetzt eine - | F 
Differentialgleichung mit konstanten Koeffizienten. Die Lage « = ist daher labil, die Lage 
«—r stabil. (Hier ist Hirsch [2] ein Irrtum unterlaufen.) Ist cos« = nn ‚soist = = | 
und |y|= a |sin al. Wegen der Grenzkurve ce, muß bei Stabilität A< a u | sein. Also herrscht 


g? 
2 nl? a® (1 — 2a) 
gewichtslage x = symmetrische stabile Lagen in den unteren Quadranten. 


sicher Stabilität für n® > Dies sind dann zwei zur natürlichen Gleich- 


5. Niehtharmonische Oszillation des Aufhängepunktes 


. Zum Schluß soll noch kurz der Fall der nichtharmonischen Oszillation des Aufhängepunktes 
betrachtet werden. Wir gehen von Gleichung (3) aus, die noch allgemein gilt. Sie ist jetzt eine 
allgemeine inhomogene Hillsche Differentialgleichung. Das Floquetsche Theorem gilt auch 
jetzt. Die Struttsche Karte ist jetzt je nach der Koeffizientenfunktion verschieden. Für einige _ 
Funktionen, z. B. Treppenfunktion [8] und f(w)=A-+ yı cos u-+ y, cos 2u [9] sind dieGrenz- 
kurven der Parameterebene diskutiert. Weitere findet man in der umfangreichen Literatur über 
die Hillsche Differentialgleichung verstreut ([5] Literaturverzeichnis). Die Existenz der 


periodischen Partikularlösung [6] gilt ebenfalls ungeändert. Entscheidungen über 1-4 


(k —=0,1,2,...) müssen für jede Koeffizientenfunktion gesondert durchgeführt werden. Sonst 
ist die oben angegebene Methode im vollen Umfang anwendbar. e . 
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Elementare Lösungen der Wellengleichung isentropischer 
Gasströmungen 
Herrn Joseph Lense zum 60. Geburtstag am 28. 10. 1950 gewidmet 


| von Robert Sauer in München 


Die Wellengleichung der eindimensionalen nichtstationären isentropischen Strömungen idealer Gase 
läßt sich bekanntlich lösen mit Hilfe einer hypergeometrischen Funktion als Riemannsche Funktion; im 
Falle n-atomiger Gase sind die Lösungen durch elementare Funktionen gegeben. Es wird gezeigt, daß die- 
selben Lösungen auch für Gase mit allgemeinerer Adiabatengleichung gelten und daß man analoge Ergebnisse 
auch für die Wellengleichung der zweidimensionalen stationären isentropischen Strömungen erhält. 


The wave equation of one-dimensional unsteady isentropie flow of ideal gases can be solved, as it is well 
known, by means of a hypergeometrie function as Riemann function; in the case of n-atomic gases, the solu- 
tions are given by elementary functions. In this paper it is shown that the same solutions hold for gases 
obeying more general adiabatic laws and, furthermore, that one obtains analogous results for the wave equation 

- of iwo-dimensional steady isentropic flow. 


L’equation d’ondes des. Ecoulements non stalionnaires isentropiques a une dimension des gaz ideaux 


m 


est resolue, comme il est bien connu, & l’aide d’une fonction hypergeometrique comme fonction de Riemann;. 


dans le cas des gaz a n atomes, il existe des solutions representees par des fonctions Elömentaires. Il est de- 
montre, dans la note presente, que les m@mes solutions sont valables aussi dans le cas de gaz ob£issant ü des 
lois d’adiabate plus generales, et qu’on obtient des resultats analogues pour l’equation d’ondes des Ecoulements 
siationaires isentropique a deum dimensions. 


E We BonHoBoe ypaBHeHHe ONHOMEPHOTO HECTANHOHAPHOTO NSCeHTPOIHYECKOrO HOTOKA HNEAAIB- 
> HOTO ra3a, MOSKeT ÖBITb, KAK USBECTHO, PeImeHo IIPMH HOMOIMM THIEPTeOMeTPHyeckuX DyHRuNH 
z kak hyurıua Pumana. B cayyae n-aroMHOTO ra3a pemmeHnsa IpencTaBIBIOTCH HIeMEHTAPHEIMH 


F $Pyueımamu. B aToit padore nokassIBaerca, YTO AA TAB0B, HONIMHAWINUXCH BoNee ODINEMy 
3AaKOHy anuadaTbl, ACHCTBH TEIBEBI TE-e caMae pelmeHng, a YTO B CILYYAE BOAHOBOTO VPaBHeHHA 
ABYyMEPHOTO CTAUHOHAPHOTO HBEHTPOHNYECKOTO IOTOKA HONYYAWTCH AHAJIOTHYHEIE PEe3YJIBTATE. 


$1. Einleitung 

Die Wellengleichung für die eindimensionalen nichtstationären isentropischen Strömungen 
idealer Gase mit der Druck-Dichte-Beziehung p = const. o*läßtsich bekanntlich!) in den Spezial- 
fällen <= ni (n—= —1,0,1,2,...) explizit lösen. Im folgenden wird gezeigt, daß die 
expliziten Lösungen nicht nur für die idealen Gase sondern auch für hypothetische Gase mit 
allgemeineren Druck-Dichte-Beziehungen gelten ($2) und daß sich analoge explizite Lösungen 
er auch für die Wellengleichung der zweidimensionalen stationären isentropischen Strömungen 
| gewisser hypothetischer Gase angeben lassen ($ 3). 
’ Anschließend werden die Sonderfälle na = 0 und n = — 1 diskutiert, wobei auch auf die 
= praktische Anwendung dieser Fälle zur Approximation von Strömungen realer Gase mit belie- 


biger Druck-Dichte-Beziehung eingegangen wird ($4). Zum Schluß wird auf Zusammenhänge 


IB: mit der Theorie der Differentialgleichungen mit geradlinigen Charakteristiken hingewiesen ($ 5). 


$ 2. Nichtstationäre eindimensionale Strömungen 
Bei der eindimensionalen nichtstationären isentropischen Strömung eines Gases mit der 


- Druck-Dichte-Beziehung p = p (o) und der dadurch bestimmten Schallgeschwindigkeit a = v* 


BL Jäßtsich ein Strömungspotential f (A, «) einführen?) als Funktion der charakteristischen Zustands- 
u variablen 


p 

21 5 ap 

2 +u+4 mit A= Jas ( ) 

Ä « bedeutet hierbei die Gasgeschwindigkeit und es wird ee daß A nicht konstant, 
I% also a eine Funktion a(A) ist. Durch Auflösung von (1) ergibt sich 


, u=A—u, A=A+tu, also a=a(li+u). 
Nach Ermittlung des Strömungspotentials f (A, #) ist der Strömungsverlauf (% = Ortskoordinate, 
En, Zeit) durch 


Er ; a— u a+u 1‘ 
| ah ee ARE) 
Yollständig festgelegt?). 
2) I. Hadamard, Teen sur la propagation des ondes. Paris 1903, 8. 163, 168. — G. Darboux, 


surfaces, II. Paris 1915, S.65. — Love und Pidduck, Phil. Trans, Roy. Soc. London, A 222, 
— 226. 

'gl. z. B. BR. Sauer, Thöorie des ee des fluides compressibles. Verlag B6ranger, Paris 1951. 
Er ’ 


340 Sauer, Elementare Lösungen der Wellengleichung a % ai 
Das Strömungspotential genügt der 


ee > ch +0 


| 
wobei der Strich die Ableitung der rer die Druck-Dichte-Beziehung p(e) gegebenen Funktion 2a 
a(A-+-u) nach A-+-u bedeutet. | 

Ein Spezialfall der Wellengleichung (3) ist die Darbouxsche Differentialgleichung \ 


v 
| Aut, Ari —0 - (4) | 
mitn=—1,0,1,2,.... Diese besitzt die expliziten Lösungen?)?) 1 
S=La)+MW für n—0 | 
4 , r 
S-Ia+MW— Fra +MW für n=—1 Be 
RO Er a M = Ede ER 
| EU EN 2 8 RER lo 
2 mit den hinreichend oft differenzierbaren, sonst aber willkürlichen Funktionen L(}), M(n). CR 
8 Wir fragen nun nach den allgemeinsten Druck-Dichte-Beziehungen »(o), für welche die 4 
2 "Wellengleichung (3) sich zur Darbouxschen Gl. (4) spezialisiert, d.h. nach denjenigen Gasen, 
Be deren eindimensionale nichtstationäre Strömungen durch die Gl. (5) explizit gegeben sind. 
a Die gesuchten Druck-Dichte-Beziehungen ergeben sich aus der Bedingung 
1—ad nn 
20 2 At-u 
Die Lösung dieser Differentialgleichung ist 
4 an | 
re (6) 


mit der Integrationskonstanten (Ü. 
Um die durch Gl. (6) festgelegten Druck-Dichte-Beziehungen zu finden, stellen wir die 


Dichte o und den Druck p ebenfalls als Funktionen von A dar. Aus 
RL REAL EERZNN 


ao doao 0 


erhält man | 
do _ dA_ ae dA 
e a aaR, 
C+ er A 
also j 
j > o(o -H Ber, Auen) mit ÖC=cont.+0 ...(7). 


Hiernach liefert 
Altn \2 
ee Fears 
| dp= ado= © (04 Da 
sofort 
A?n+3 


PT Glan Onf3) Ir 5, es, 


Durch Elimination von A aus den Gl. (7) und (8) folgt lea “ ey Druck-Dichte- 
Beziehung PEN, 


er : r EN s ? (ee # 


= En 0% Da 
Pt rt Firm® + 
mit Be 
2n+3 
„Et see 2 
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Dichte-Beziehung vor, für welche dieselbe Wellengleichung (4) mit demselben Lösungen (5) 
wie bei den idealen Gasen gilt. 

Zusatzbemerkung: Für nicht ganzzahlige Werte von n hat Riemann für die Lösungen 
der Wellengleichung (4) eine Integraldarstellung mit hy pergeometrischen Funktionen angegeben 
und durch diese Lösungen die Ausbreitung endlicher Druckwellen in idealen Gasen beschrie- 
ben?). Die vorangehenden Untersuchungen zeigen, daß die Riemannschen Lösungen nicht 
nur für die idealen Gase, sondern ebensowohl für die durch Gl. (9) en hy pothe- 
tischen Gase gelten. 


$ 3. Stationäre zweidimensionale Strömungen 
Wir führen nun analoge Untersuchungen durch für die zweidimensionale stationäre isen- 
tropische Gasströmung. Hier kann man ein Strömungspotential f(A, «) einführen?) als Funktion 
der charakteristischen Zustandsvariablen 
2al 


2 e ea zu 
EN mit W= fe N 


w ist der Betrag und ® der Winkel des BD windipkeitspektors und es wird vorausgesetzt, 

daß W nicht konstant, also w eine Funktion w(W) von w ist. Außerdem soll w>a sein, d.h. 

wir beschränken uns auf Strömungen mit Überschallgeschwindigkeit; dann sind die charak- 
7 teristischen Variablen reell. Durch Auflösung von Gl. (10) ergibt sich 


d—=)—u, W=i-+u also v—=w(-+u). 
Nach Ermittlung des Strömungspotentials f(A, u) ist der Strömungsverlauf (x, y = Ortskoor- 
dinaten) durch 


22= „008 An) (hHf) — 4 sn mh) 


(11) 
148, 1 
2y= „sin %— u) (fHf) +, os Au) (fi) 
vollständig festgelegt. 
Das Strömungspotential genügt der Wellengleichung 
vw" +w 
Sau — a N Pr Nr, 


wobei die Striche Ableitungen der durch die Druck-Dichte-Beziehung p(o) gegebenen Funktion 
w(A-+-u) nach A-+u bedeuten. 

Wir fragen jetzt wieder nach den allgemeinsten Druck- Dichte- Beziehungen p(o), für wel- 
‚che sich die Wellengleichung (12) zur Darbouxschen GI. (4) spezialisiert. Sie ergeben sich aus 
der Bedingung 


BER, AT ER RE Er EN: 


aus der man (vgl. E. Kamke: Difierentialgleichungen, Lösungsmethoden und Lösungen, 1, 1942, 
S. 425, 2.105) 


' P) Ar A L 
ol) cos (A+ u)+ (0, Fewerreen ) sinA+u) für n20 (14) 


|e= [la -+u+C,)sin(A+uw)+cos(A+W] für n-—ı 


erhält. Durch Gl. (14) ist in ähnlicher Weise wie in $2 die Druck-Dichte-Beziehung »(e) be- 
stimmt. Wir beschränken uns darauf, diese im Sonderfall n = 0 in $4 herzuleiten. 
Während die bei der nichtstationären eindimensionalen Strömung eingeführten, durch 
Gl. (9) gekennzeichneten hypothetischen Gase als Spezialfall die idealen Gase enthielten, sind 
-  dieidealen Gase unter den bei der stationären zweidimensionalen Strömung eingeführten, durch 
© - Gl. (14) gekennzeichneten hypothetischen Gasen (— vom Sonderfall n = 0 abgesehen —) nicht 
enthalten. 
 Zusatzbemerkung: Für die hypothetischen Gase, welche durch die Differentialglei- 
chung (13) mit nicht ganzzahligen Werten von n gekennzeichnet sind, lassen sich die zweidimen- 
S Panlan stationären Strömungen mit Hilfe der Riemannschen Lösungen?) der Darbouxschen 
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$ 4. Diskussion der Sonderfälen=0 undn= —1 


Im Fall n = 0 erhält man bei den eindimensionalen nichtstationären Strömungen an Stelle 
der Gin. (6), (7) und (8) 


a=4+Q, o=04+0Q, g=p+S(4+0) 


und daher an Stelle der Gl. (9) die Druck-Dichte-Beziehung 


\r=9+ 08° | re ed en RE Te (15). 


Hier reduziert sich die Zahl der Parameter von 3 auf 2, es ergibt sich also im wesentlichen nur die 
Druck-Dichte-Beziehung. der idealen Gase mit «=3. 
Bei den zweidimensionalen stationären Strömungen spezialisiert sich Gl. (14) im Fall 


==) i 
4 " v=(,csA+u)+(,;sinA+u)=w,cos(A+u—d) -» .»....(16) 
mit wiederum nur zwei Integrationskonstanten w,, ®,. Um hieraus die Druck-Dichte-Be- 
ziehung herzuleiten, bilden wir zunächst nach Gl. 1. und (16) 


Yu: ey a? 


Dann ergibt sich aus der Bernoullischen Gleichung 


= PA (wi — u?) . 3 ® 


5 3 
a 
für die Dichte de dp do dp 1 vwdu  uldw h 
ee ap oa an ee 
also 
e= ww 


und für den Druck 


2 
dp=ade-(") (do Org ee 


2, (oe E 09% 
also Ba wen Ne ae Na Sn 
2 oO Ss E Be Te Fe . daz 
P=Po+t Co (a —aretgo) (10 
Im Falln— —1 erhält man bei den eindimensionalen nichtstalionären Strömungen in > 
analoger Weise P A 
‘| a—= A(CA—]1), „_ C(0A—1) dp = 0(0A— 1) dA > 


4A ’ 


und hieraus die Druck-Dichte-Beziehung 


PM — bei 0105 
P= ° s Be (18) 
re. bei U +0 7 Br - 
+ (2, +0(,0 Tr 
Bei den zweidimensionalen stationären Strömungen mit n—= — 1 lassen sich die zur Aufstellung 
der Druck-Dichte-Beziehung erforderlichen Eliminationen nicht explizit durchführen. ; 
Die Sonderfälle n=0undn= Fe an man dazu verwenden, um ar realer % 


mungen mit numerischen Be des Einzelfalls ne Zur Gewinnung er er 
Einsichten ist es daher sehr wünschenswert explizite Darstellungen von Strömungen zur Verfü- 5 
gung zu haben, die sich in verhältnismäßig einfacher Weise den vorliegenden Anfangs- und. Ran. 12 
bedingungen anpassen lassen, Die Sonderfälle n=0 und n = — 1 liefern solche ee 


la BR A A Ben u. es la BE DE se 
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Parameter kann jedoch eine gewisse Approximation der realen Druck-Dichte-Beziehung er- 
reicht werden. Die Druck-Dichte-Beziehungen (15) bzw. (17) des Falles n = 0 enthalten 2 Pa- 
rameter und ermöglichen dadurch eine Approximation erster Ordnung, d.h. eine Übereinstim- 
mung der realen.und hypothetischen Druckkurven p = p(e) in einem Punkt und der Tangente in 
diesem Punkt. Die Druck-Dichte-Beziehung (18) des Falles an = —1, C = 0 enthält 3 Parameter 
und ermöglicht dadurch sogar eine Approximation zweiter Ordnung, d. h. eine Übereinstimmung 
der Druckkurven in der Krümmung. 

Die Druck-Dichte-Beziehungen (15), (17) und (18) sind in anderem Zusammenhang schon 
mehrfach bei gasdynamischen Untersuchungen verwendet worden. Die Beziehung (15) findet 
sich in Arbeiten von Bechert und Pfriem, die Beziehung (17) wurde gelegentlich von Prandtl. 
Munk und P£r6s, die Beziehung (18) von Sauter benützt. 

$ 5. Beziehungen zur Theorie der Differentialgleichungen mit geradlinigen 
Charakteristiken 

‚Im Vorangehenden wurden explizite Darstellungen von Gasströmungen dadurch gefun- 
den, daß die Wellengleichung (3) bzw. (12) in die Form der Darbouxschen Differentialgleichung 
(4) gebracht wurde. Eine andere Möglichkeit, zu expliziten Darstellungen von Gasströmungen 
zu gelangen, bietet die Theorie der Differentialgleichungen mit geradlinigen Charakteristiken®). 
Diese bezieht sich auf lineare Differentialgleichungen 


A&n)fe+t BE tlt&EmMm=0 2.2... (N), 


die in den Ableitungen weiter Ordnung homogen sind, und auf die durch die Legendre-Trans- 
formation — ei 
E=ff, H=f, F=ES4nH—f, E=R, 9=M 


_ zugeordneten quasilinearen Differentialgleichungen 


ET AR E)Eas C(F5, Fu)Fr2,=0N. (00% 


Wenn die Charakteristiken der Gl. (19) geradlinig sind oder ein Rückungsnetz mit 2 Scharen- 
paralleler Kongruenter Kurven bilden, sind die Charakteristikennetze der Gl. (20) Rückungs- 
netze bzw. Geradennetze. In beiden Fällen lassen sich die Lösungen der Differentialgleichungen 
(19) und (20) mit Hilfe zweier willkürlicher Funktionen explizit angeben‘). 

Man kann nun explizite Darstellungen von Gasströmungen dadurch erhalten, daß man die 
Wellengleichungen (3) und (12) durch Auswahl geeigneter Druck-Dichte-Beziehungen in die 
Gl. (19) zu transformieren sucht, derart, daß die Charakteristiken in der &, n-Ebene ein Geraden- 
netz oder Rückungsnetz bilden. Die allgemeine Durchführung dieser Aufgabe soll einer späteren : 
Arbeit vorbehalten bleiben. Wir wollen ohne auf die Beweise einzugehen, lediglich auseinander- 
setzen, daß man als Sonderfälle von neuem die in $4 diskutierten Strömungen mit n = 0 und 
n=—-1l erhält: 

Im Fall n = 0 der zweidimensionalen stationären Strömung, der durch die Druck-Dichte- 
Beziehung (17) gekennzeichnet ist, kann man als &,n- -Ebene die u, v-Hodographenebene (u, v—=car- 
tesische Geschwindigkeitskoordinaten) und als #,H-Ebene die x, y-Strömungsebene neh- 
men. Die Charakteristiken in der w, v-Ebene bilden ein Rückungsnetz, das durch Rotation eines 
Kreises um einen Punkt desselben erzeugt wird. Die Gleichung dieser Charakteristiken in Po- 
larkoordinaten ergibt sich aus .Gl. (16), indem man unter Berücksichtigung von A—-u=din 
Gl. (16) an Stelle von A+ uentweder 2 A — d oder 2 u + Veinführt und A = const bzw. u =const 
setzt. 

Im Fall n= 0 der eindimensionalen nichtstationären Strömung, der durch die Druck- 
Dichte-Beziehung (15) 2 Ei kann man als &,n-Ebene die u, g-Zustandsebene 


(u= Gasgeschwindigkeit, 1=— — Be =) und als, H-Ebene die x, {-Ebene wählen. Die 


v 
Charakteristiken in der «, g-Ebene bilden ein Rückungsnetz, das von kongruenten und parallelen 


- Parabeln erzeugt: wird und als Einhüllende eine parallele Parabel doppelter Größe besitzt?). 


Im Fall n=— 1 der eindimensionalen nichtstationären Strömung, der durch die Druck- 


Pa | Dichte-Beziehung (18) gekennzeichnet ist, hat man die Unterfälle O= 0 und Ü=+0 zu unter- 


scheiden. Bei O = O0 kann man als &, n-Ebene und Z, H-Ebene wieder die «, q- und z, t-Ebenen 
‚benützen und erhält dann in der vw, q-Ebene ein geradliniges, von den Tangenten einer Parabel er- 


s- zeugtes Pe nsikenretz). Bei Ü +0 ergibt sich dasselbe Charakteristikennetz, jedoch hat 
s Ber jetzt 


f je 2a 
es = (dm), &=2—(1+4) u, H=— 77 t zusetzen. 


’ 


TER. Sauer, Z. angew. Math. Mech. 25/27, 1947, 8. 151153. 
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Der Übergang vom ebenen Spannungszustand zum ebenen 


Formänderungszustand im breiten gebogenen Balken‘) 
Von @. Sonntag in München 


Unter Verwendung bekannter Lösungsansätze für Platten wird der Übergang vom ebenen Spannungs- 
zustand zum ebenen Formänderungszustand beim breiten gebogenen Balken allgemein streng angegeben und 
für einige Beispiele berechnet. 


By use of the known solution for plates, the transition from the plane state of stress to Ihe plane stale 
of deformation is given, generally and ewactly, in Ihe case of Ihe bending of a broad beam, and illustrated 
by some examples. 


On fait usage des solutions bien connues pour des plaques pour donner exaclement et de maniere 
generale, dans le cas d’une poutre large courb£e, la transition d’un Etat de tension plane dans un £tat de 
deformation plane; quelques examples sont donnes. 

Ilepexon OT INIOCKOTO HANPSIKEHHOTO COCTOAHHA K_INIOCKOMY COCTOAHHIO MebopMmaumi B 
IINPOKHX ÖAJIKAX, HATPY>KCHHEIX HA H3TUÖ, uecHenyeTca BoOmeH u CTPoroä dopMe pm HOMOLIM 
USBECTHBIX pemenni ANA MISCTHH H MINOCTPHPYeTcH HA HECKONBKUX YACTHEIX CIIYYAAX. 


I. Einführung 


Wir betrachten nach Bild 1 einen gewöhnlichen rechteckigen Balken der Stützweite 2 a. 
Seine Breite 2b sei klein gegenüber seiner Stützweite 2a. Unter einer Belastung p=f(x) wird 
sich der Balken durchbiegen, die Biegespannung überschreite nicht die Elastizitätsgrenze des 
isotropen Werkstoffes. Wenn die Querdehnungszahl »+0 ist, erfahren die Druckfasern eine 


Bild 2. Querverfor- 


ee Bild 3. Plattenstreifen, zweiseitig 


Spannungszustand . frei drehbar gelagert 


Bild 1. Schmaler, zweiseitig frei 
drehbar gelagerter Balken 


Querdehnung, die Zugfasern eine Querkontraktion, solange diese Querverformung nicht be- - 


hindert ist. Lange schmale Balken sind in ihrer Querdehnung unbehindert, ihr Querschnitt 


verformt sich nach Bild 2. Spannungen treten nur in der x-,2-Richtung auf, die y-Richtung ist’ 


spannungsfrei; deshalb sprechen wir von einem ebenen Spannungszustand. 


” Wenn die Breite 2b des Balkens immer mehr zunimmt, liegt zuletzt der unendlich aus- 
gedehnte Plattenstreifen nach Bild 3 vor. Die Belastung pseivon Y unabhängig. Von 
den x-z-Ebenen ist keine irgendwie ausgezeichnet; sie sind alle vor und nach einer Belastung 
p=f(x) parallel ohne Änderung ihres Abstandes. Denkt man sich den Plattenstreifen in 
le Balken nach Bild 2 aufgeschnitten, so können diese wegen ihres Zusammenhanges 
keine Querverformung erfahren, sie bleiben eben nach Bild 4, man spricht vom ebenen 
Formänderungszustand. Dieser wird jedoch erzwungen von Span- 
nungen in %-Richtung. Es läßt N Re RR, daß hier die 
Beziehung o,—=»o, und M,=»M, 


Ersatzscherkräfte 05—=0). Es liegt dort der ebene Spannungszustand vor. 

Bild a. Momente Mur Mit zunehmendem Abstand vom freien Rand werden sich die jetzt auf- 
Pr \ « Yy 

Erzwingung des ebenen tretenden Spannungen o, den Spannungen des ebenen Formänderungs- 


Kormänderungszustandes Zustandes nähern. Bevor die strenge Lösung gegeben wird, soll noch eine => 


weitere Betrachtung vorgenommen werden. 
Eine völlig freie Platte werde an ihren Rändern &= --adurch ein gleichmäßiges Moment; 
belastet (Bild 5). Die Platte verformt sich dann zu einer Sattelfläche mit den ‚Krümmun 


') Die Untersuchung wurde durchgeführt nach Anregung von Herrn Prof. L. Pr andtlauf Eimpf pfehlu 
von Herrn Prof. L. Föp pl, wofür ich beiden Herren verbindlichst danke. k "X r 


Wenn der breite Balken jetzt iR en Breite 2b hat, so ist an 
den Rändern „= +b, o,—=0 und M,=0 (selbstverständlich sind auch die 
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1 1 M, 
3 a= und u der Größen n=TJ 090 = —V0j0, wobei der Index 0 bedeutet, daß 


\ - 2 
es sich um die völlig freie Platte handelt. Wenn jetzt die Platte an den Rändern 2=-+«a 
nicht mehr frei, sondern gelenkig gelagert ist, werden die Krümmungen behindert und zwar 
um so mehr, je länger die gelagerte Seite im Verhältnis zur Stützweite ist. Für den Platten- 
FR: streifen, ba —= x, ist die Krümmung 0,=0. Es 
zeigt sich nun, daß, je mehr die Krümmung os 
durch Kräfte am gelagerten Rand herabgesetzt 
wird, auch die Krümmung o, verringert wird; 


Bild 5. Sattelfläche des völlig freien durch M,„ Bild 6, Steifigkeit M „Ja des breiten Balkens als Funktion der 
) belasteten breiten Balkens Krümmung o0,. 03 = 0 entspricht dem ebenen Formänderungs- 
zustand 


was somit als eine Erhöhung der Plattensteifigkeit in Erscheinung tritt. Es soll nun die 
- Plattensteifigkeit M,/o, als Funktion von 0,/0, bzw. — vo,,/o, dargestellt werden. Es gilt 
01 = 010 + ? (020 — 03), wodurch der Einfluß der Krümmungsverringerung 09 —0,, infolge der 
Lagerung, auf die Hauptkrümmung o,, ausgedrückt wird. Es folgt nach einfacher Umformung 


& —_ 122% und die Darstellung (Bild 6) ENTER für den Bereich 


010 Q10 _ [41 1 — »? — 903/010 
02 = O1 = — VO bis 90. 

Nach dieser Vorbetrachtung besteht die Aufgabe nun darin, die Krümmung o, bzw. das 
Moment M,, als Funktion von y zu bestimmen, um daraus zu erkennen, wie und in welchem 
Maße der Übergang vom ebenem Spannungszustand am Rand y=-+b der Platte zum ebenen 
Formänderungszustand innerhalb der Platte erfolgt. Nachfolgend wird die strenge Lösung 
dieser Aufgabe gegeben. 


II. Die allgemeine Lösung 


Zunächst wird entsprechend einem unendlich breiten Balken nach Bild 3 für eine beliebige 
Normalbelastung p=f(x) nach der gewöhnlichen Balkenformel M, berechnet. Für die Durch- 


’ biegung wäre wegen der behinderten Querdehnung E/(1 — »v?) an Stelle von E zu setzen. Aus 
‘gleichem Grunde tritt, wie zuvor näher ausgeführt, M,=»M, auf. Schneidet man jetzt aus 
> - diesem Plattenstreifen einen Balken endlicher Breite 2b heraus, so muß dem Schnittmoment M,, 
} - ein gleiches Moment umgekehrten Vorzeichens überlagert werden, um den Rand y= + lastfrei 


zu erhalten. 

Die Aufgabe ist somit auf die Zusatzlösung der längs der Seiten &—=-+a« drehbar ge- 
lagerten oder eingespannten Platte zurückgeführt, deren freie Ränder y=-+b durch gleiche 
Momente M,= f(x) belastet werden (s. Bild 7). 


Der allgemeine Ansatz lautet?): 


w— 53 A, cos vAx (Ay Sin eAy-+n;ojeiy) | 
i=0,1,2,8 +. 06 2 WE HE 
+ 2 ‚B, sin JA (jAy Sin jay+n, Coijäy)| 
i=2,4,6..- 


worin w — Durchbiegung, A 32: A, B, n = Beiwerte. 
73 


Es läßt sich zeigen, daß der Ansatz die Bipotentialgleichung AAw=0 erfüllt. 
Baar 2 j 
\ g m)" Die Beiwerte A und B werden aus den Auflagerbedingungen und der 


"Nachfolgende Funktionen sind zur Lösung von Plattenproblemen längst bekannt. Nach vielseitiger 
ntierung des Verfassers scheint die hier behandelte Sonderfrage noch nicht durchgeführt zu sein. 


Due 2 
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Randbelastung M,— f(x) bestimmt. Die Beiwerte n gibt die Forderung nach Verschwinden ER h 
der Ersatzscherkräfte On a an den freien Rändern y„=-+b. (Die Ersatzscherkräfte resultieren 
aus der Zusammenfassung der Quer- 

kräfte und Torsionsmomente.)?) 


Der Ansatz (1) wird unbequem für 
große Werte bja. Für b> 3a verwendet 
man zweckmäßig den unendlich ausge- h 
dehnten Platten-Halbstreifen nach Bild 8. 
Es folgt entsprechend (1) 


u= 23 A, cos el —ny)eir 


2 end Ber -jAy 
a a Bild 8. Plattenhalbstreifen unter a B; sin JA@a(l —n,Y) eT 
an rtpelsstnngder bien Momentbelastung des freien j=2,4,6°.. 
Ränder andes (2). 


Damit kann die Aufgabe für beliebige Randbelastung M,—=f, («) bzw. M,=fs ?@) als gelöst be- 
trachtet werden. 


III. Beispiele - 


a) Der gelenkig gestützte breite Balken mit gleichmäßiger Balasrnas ©) 

Das Beispiel gleichmäßiger Belastung und gelenkiger Stützung wurde gewählt, weil sich B 
hierbei das Moment M, in genügender Annäherung durch cos (rz/ 2a) ausdrücken läßt, so daß 
man mit dem ersten Glied der anzusetzenden Reihe auskommt. Bei einer Belastung der Platte 
durch ein reines Moment M, wäre auch am freien Rand ein konstantes Moment M„=vM,=const. 
anzubringen, wofür die Rechnung sehr viel komplizierter wird. Dieser Fall wird später noch 
ausführlicher besprochen. 

Für den unendlich breiten Balken gleichmäßiger Belastung p gilt bekanntlich, wenn 
Koordinaten nach Bild 3 zugrunde gelegt werden: 


FE RT 


x 
M,= My; 1-3); Mu=WwM 12 Men ee 
Bir: ler 1 | \ 
A a w—, M,0 = 19° (ee, a? a1) < N ER re (4), 
’ : Br worin | N Sr HRZE nz Balkenhöh. 
Be | TI: RE EEE 


, Um die Rechnung zu vereinfachen, setzen wir mit geringem Fehler anstatt (3) 


M,= M,cosAs mit I=5n; Mer Mu. a 5), 
damit wird 
Er er ji ! M, 
V-—iay cos Az EN RN 0 6 (6), 


Für %,.. beträgt Bis Abweichung nach (6) gegenüber 4) zw — 2,5%. i 
Zu unserem neuen Ansatz (5), (6) gehört streng eine Belastung p—=c: cosAz. Obwohl 
wir eingangs eine gleichmäßige Belastung p angenommen hatten, dürften die Ansätze (5), (6) 
unsere Frage nach dem Übergang vom ebenen Spannungszustand zum ebenen Formänderung- 
zustand befriedigend beantworten, da sowohl die Momentenverteilung M, als auch die Durch- 
biegung w dem Fall gleichmäßiger Last weitgehend angenähert wurden. « EM 
Gehen wir jetzt zum Balken endlicher Breite 25 über, so folgt für die anzubringende fe 
Randbelastung M, nach (1): | 


w=Acos/As (AySinAytn Gofay) ee a ee ze 
Die Stützbedingung am Rand <= a, w=0 und M,=0 ist erfüllt. A bestimmen, wir au 
Bedingung (M,) „-+»—= M,, cosAx. Es folgt mit den A BRIZ en h x 
unit 
M,= — N{w"+ vu") = — NAcosAz{Ay SinAy-+A2(2 +n)&o Cof y 
rrl—AySnAy—Anbofay) R 5 ) 


®) Vgl. etwa die Ausführungen in C. B. De und R. Sram Techn, 
rn 1939, 8.436 oder Beton-Kalender (1944), S. 129. 
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Daraus für y=+tbund s=0: 


y 1, ui: 
PTBA Sin bLP FRA ICN 0). 
Die Ersatzscherkräfte 
= — No + @— ur] 
= —NAcosAaz{AyCofiy+ (3 +n) SinAy RAN) 
neh + an) [—Aybofiy—R(l+n)SinAyl} R 


verschwinden für „=0 wie aus Symmetriegründen erforderlich, aber auch für y=--b muß 
Q,=0 sein, daraus n wie folgt 


„= U—mAy Cory Hr) Sinay_ —AbGtgAb-+1 ach ESTER ER 


= . 1 —») Siniy 


Nachdem für einen gegebenen Wert b bzw. b/a die Beiwerte n nach (11) und A nach (9) ermittelt 
sind, kann der Verlauf von M, nach (8), die Biegung w nach (7) oder jede beliebige Schnittgröße 
nach den bekannten Differentialformeln der Plattentheorie berechnet werden. 


Plattenstreifen a)b=0 


Schmaler Balken 
b —o b 

2 

Bild 9. Verlauf von M, über die Breite 25 des breiten zweiseitig frei drehbar gelagerten. Balkens mit gleichmäßiger Belastung in 2 
Abhängigkeit von der Stützweite 2a. My, =vM, „ entspricht dem ebenen Formänderungszustand f 


In Bild 9 wurde für Abmessungsverhältnisse b/a—=0,25; 0,5; 1,0; 2,0 und co der Verlauf ; 
des Quermomentes M, : (M,) = » im Mittelschnitt &—=0 dargestellt (Rechenschiebergenauigkeit |) / 
> für v=0,5. In der Darstellung wurde die Auflagerbreite 2 5 konstant gelassen und die Stütz- 
j "weite 2 a variiert. Für andere Schnitte & = const ergibt sich hier das gleiche Bild, wenn M,,, jetzt 
allgemein das Moment des ebenen Formänderungszustandes bedeutet. 


w b) Breiter Balken mit reiner Momentenbelastung M,„= const. 


Wie im Beispiel a) bereits erwähnt, muß der zu überlagernde Ansatz der Durch- 
biegung w (1) so gewählt werden, daß sich am freien Rand ein konstantes Moment M,, ergibt. 
Dieser Bedingung M,= const für y= + b genügt der Ansatz?) entsprechend (8): 
W=—D At coside (RPy Siniay+ (24 m) Eofiäy | 

i ; 


| . (12) | 
4r [A ıdy Sin say — 122m, Cofiryl} N 


miti=l —3, +5, —7. 
A; und n, werden für ade Glied für sich ee entsprechend (9) und (11) wenn darin 


 ansta jetzt ©A gesetzt wird. Für den freien Rand y=+« folgt dann nach Voraussetzung x i 
M,= Mu‘ = [esan— zo Bio+ cos 5A2— 4.00 7ict'- - | M,= ost‘) f 
— — 


A 9 Val. Hütte, I, 27. Aufl., 8.207. 


. der gleich der Stützweite des Plattenhalbstreifens ist. 


. ein wesentlich schnellerer Übergang zum ebenen Formänderungszustand zu 


\ } „7 . ar R Ein Rah Bar hr a in Den 
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Die Differentialbeziehung zwischen M, und w läßt den Ausdruck für w leicht erkennen, es folgt 


4x 1 1 : EBEN, 3 
BD A En (13) 


Gl. (8) läßt erkennen, daß bei gleichmäßiger Belastung (oder streng für die Belastung 
p=c.cosAz, dem Gl. (8) entspricht) die Verteilung von M,— f(y).cos A& für alle Schnitte 
parallel zum freien Rand y= + b ähnlich ist. Das trifft bei dem zunächst einfach erschei- . 
nenden Fall konstanter Momentbelastung nicht mehr zu, da die Glieder der Reihe (12) in y ver- u 
schieden schnell abklingen. : 


| 
mise 28 RAN, E- 
j Br 


| c) Der zweiseitig gelenkig gestützte Plattenhalbstreifen 
mit gleichmäßiger Belastung 


Wenn 5b> 2a wird, können wir den Einfluß des zweiten freien Randes vernachlässigen. 
Wir verwenden an Stelle von (7) den entsprechenden Ansatz nach (2): 


4 
E 
: 
Fir 


w—= AcosA2.(—nYy)e Ar na 


wobei jetzt nach Bild8 y vom Rande aus zählt. 
Es folgt entsprechend (8): 


M,=— NAcosAn: e-W{R (IL —ny)+2A—r®li—ny)}: - - - - (IB). 
Daraus folgt für y=0: 


M, ’ 
- : 5 I Be LEN 
2 [A (1—»)+2n]AN (16) 
Die Ersatzscherkräfte 

y—= — NAcosAn- et 2 {Al ny)—3n+R—)Al—ny)+tnl} (iM 


verschwinden für y= 0, wenn 


ee EP N A Se 


Die einfache Ausrechnung des Verlaufes von M, und w vom freien Rande zum Inneren des 
unendlichen Halbstreifens fortschreitend, wurde für 20 und für die Querdehnungszahl »=0,5 
in Bild 10 dargestellt. 


3a 


ebener 


Verformungs =, . 
zustand c 


1 RN a % 
Bild 10. Übergang des ebenen Spannungszustandes (am freien Rand) zum ebenen Formänderungszustand beim zweiseitig fr 
gelagerten Plattenhalbstreifen mit gleichmäßiger Belastung. Die Kurve W/A gibt die zusätzliche Durchbiegung in R 
r . an. M,,, entspricht dem ebenen Formänderungszustand a a a a 

j x 


Der ebene Formänderungszustand wird hier in einem Abstand vom freien 
u 
Beim eingespannten breiten Balken ist wegen des Vorzeichenwechsels der Bieg 


fache Rechnung wird hier verzichtet. a 
Eingegangen am 21. November 1950. 
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Die Temperaturverteilung in Schmierschichten zwischen 
parallelen wärmedurchlässigen Wänden 
Von @. Vogelpohl in Göttingen 


Vf. zeigt, daß für die Couelte-Strömung einer Flüssigkeit, deren Zähigkeit nicht von der Temperatur 
abhängt, für die Berechnung der Mittelwerte der Temperatur die Strömungsgeschwindigkeiten quer zur 
Schicht gemittelt werden dürfen, und daß dazu die Wärmeleitung quer zur Schicht vernachlässigt werden darf. 


The author, considering the Couette flow of a fluid, the viscosity of which does not depend on the tempe- 
rature, proves lhat, for the computation of the mean values of the temperature, it is possible lo use Ihe mean 
values of the velocity of the. flow taken across the layer and, for that purpose, to neglect the heat conduction 
in that direction. 


L’auteur traite l’&coulement de Couette pour un fluide, dont la viscosileest independanle de la temp£- 
ralbure; il montre qu’on peut, pour un tel fluide, Evaluer les valeurs moyens de la temperature en employant 
les valeurs moyens de la vitesse du flux & travers la couche, et qu’il est permis, & cette fin, de negliger la 
‚sonduction de la chaleur enn cette direction. 

ABrop uecaenyer NOTOR KyaTrTa MIA SKHAKOCTH, BABKOCTB KOTOPOH He 34BHCHT OT TEMIIE- 
P&Typbl, AH TMOKAsBIBAeT, YTO TIPM BbIyHC/TeHuM CpeAHeh TeMIEepaTypEIl MOUCHO TONB3OBATBCH 
CPeAHHM 3HAYEHHEM CKOPOCTH IIOTORA B HOTMEPYHOM GEeYEHHNH CJIOA H YTO IIPH ITOM NOLYCTUMO 
IpeHeöperaTb TEeIMNOIPOBONHOCTBIO B TIOHEPEeYHOM HALPABJICHHM. 


1. Aufgabestellung 


> In der hydrodynamischen Theorie der Schmiermittelreibung wird zur Vereinfachung der 
mathematischen Behandlung die Temperatur und damit die Zähigkeit des Schmiermittels meist als 
konstant angenommen. Das ist ein Grenzfall, der entweder langsame Bewegung oder sehr große 
Wärmeleitfähigkeit des Schmiermittels voraussetzt. Nun haben aber alle Flüssigkeiten eine 
geringe Wärmeleitfähigkeit, so daß es bei größeren Geschwindigkeiten, wenn die Erwärmung 
durch innere Reibung nicht mehr vernachlässigbar klein ist, sinnvoll erscheint, den anderen 
Grenzfall zu betrachten, bei dem die gesamte entwickelte Wärme im Schmierfilm verbleibt und 
dessen Temperatur erhöht. Dafür wer- 
den wärmeundurchlässige Wände vor- 
ausgesetzt, die es in Wirklichkeit nicht 
gibt. Andererseits stehen aber bei hohen 
Geschwindigkeiten für den Wärmeüber- 
gangin die Gleitflächen nur außerordent- 
lich geringe Zeiten zur Verfügung. Zum 
Beispiel hat Buske bei Gleitgeschwin- 
digkeiten von 73 m/sec. Temperatur- 
erhöhungen bis zu 180°C gemessen (La- 
ger für hohe Anforderungen, Bericht, 


% 


(y)=h.z=0. feste Wand 


(y)=0,z=1, bewegte Wand 


S.11 der Lilienthal-Gesellschaft, 1940, x=0 
7) AN, Abb. 38), wobei die Durchflußzeit Bild 1. Schmierspalt und Geschwindigkeitsverteilung 
Fa des Öles etwa 0,002 sec. beträgt. 
#) Vom Verfasser war der ebene Schmierkeil mit ganz in Öl verbleibender Wärme unter der 


Annahme einer mittleren Temperatur und demnach mittleren Zähigkeit quer zur Schicht be- 
handelt worden, vgl. F.H. 386, S. 19/23t). Prüft man die Zulässigkeit dieser Annahme und 
betrachtet dazu einen parallelen Schmierspalt nach Bild 1, so ergibt sich unter Vernachlässigung 
der Wärmeleitung aus der weiter unten folgenden Gl. (5) an der festen Wand ein unendlich großer 
Temperaturgradient. Herr Blok in Delft wirft daher mit Recht die Frage auf, ob man bei der 
Behandlung der nicht isothermen Strömung mit Mittelwerten der Temperatur quer zur Schicht 
rechnen darf und stellt folgende Aufgabe: 

„Mit den in F.H.425?) verwendeten Methoden ist die Temperaturverteilung einer Strömung 
zwischen athermanen Wänden unter den Bedingungen von Bild 1 zu bestimmen und fest- 
zustellen, ob der Mittelwert der Temperatur quer zur Schicht endlich bleibt, und wenn ja, ob 
dieser Mittelwert gegen denjenigen konvergiert, den man aus der Näherungsrechnung im Sinne 
von F. H. 386 erhält, bei der die Strömungsgeschwindigkeit und die Temperatur vorher gemittelt 
wurde. Dieses Ergebnis würde für die Zulässigkeit des im F.H.386 benutzten Verfahrens ent- 

scheidend sein.“ ARE 
Für den Fall einer Flüssigkeit, deren Zähigkeit von der Temperatur unabhängig ist, läßt 
‚sich das Problem vollständig lösen. Er stellt den Grenzfall der Temperaturerhöhung dar; für 


ze > Aal 


) F.H. 386 — VDI-Forschungsheft 386, Beiträge zur Kenntnis der Gleitlagerreibung; Berlin 1937. 
) F.H. 425 — VDI-Forschungsheft 425, Der Übergang der Reibungswärme von Lagern aus der 
schicht in die Gleitflächen; Düsseldorf 1949. 
‚2 he | 
Bi 24 
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ir 
3. 
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Flüssigkeiten, deren Viskosität mit zunehmender Temperatur abnimmt, werden nicht so hohe 
Temperaturen erreicht. 


2. Grundgleiehung und Partikularintegral 
Als Ausgang wird die Energiegleichung für die Strömung einer zähen Flüssigkeit [F.H.425, 
S.4, Gl.(2)] genommen, wobei die wirklichen Werte der Koordinaten (x), (y), Geschwindig- 
keiten (u) und Temperaturen (9) durch Einklammern gekennzeichnet sind. Unter Vernach- 
lässigung der Wärmeleitung in der X-Richtung lautet die Gleichung [ebda. S. 6, G1.(10)]: 


| ode): N (1) 


i 


mit dem spezifischen Gewicht y, im technischen Maßsystem in kg/m?, der spezifischen Wärme c 
in mkg/kg °C, der Zähigkeit n in kgs./m? und der Wärmeleitfähigkeit A in mkg/m. s °C. 

Es werden dimensionslose Koordinaten &, % durch (x<)= hz und (y) = hy eingeführt, 
und die Geschwindigkeitsverteilung linear in der Form der Couette- Strömung 


4 Be 
Ba in 2 une Bauml 20 0 Zn u u 


7; Y)\ _ U FRE, 

(u) = of A Ul—y), also: N 
’ 
ö angenommen, ferner die Zähigkeit „—=n, als konstant und von (9) unabhängig vorausgesetzt. 


Mit der Beziehung 
a A RE EN 


wird die dimensionslose Größe ® statt der Temperatur (9) eingeführt; nimmt man noch die 
BR; Pe&cletsche Zahl 


et N. 0? 
Gl De BE 


; pas Eh year ea N (8) 
| A 
‚hinzu, so wird nach Einführung dieser Größen die Gl. (1) 

x 0 989 

N tal ES RNER SC act" 45 . (4) 

. oder mit z=1—y | 

od 99 2 

Pers za! Re EN u A 


Die Randbedingungen sind 2(9)/2(y) = 0 für (y)=0 und (y)= h, oder 
: 
BI für 2=0 und z2=1, 
und die Anfangsbedingung 9 —=0 für 20. 

. Es läßt sich ein Partikularintegral 9, in der Form 


r- 


Bi 9%, = Ac+ 04 2 + %2?+ 032° 


angeben; die Randbedingungen und die Differentialgleichung (5) liefern 
2 18 1%" | | 
d=n,: IR 6 ie 1 Dr 2 En Ju 
also | 2% 1 0 N 
- d,=n, +45 BR — 22) a ass RB SB ee (6). "> sie 2 


\ Ar 


Die Konstante c, wird aus den Anfangsbedingungen bestimmt, für 2&—0 soll 9-0 sein. Für ie 


die Partikularlösung ist das nicht möglich. Fordert man, vorbehaltlich einer späteren Kontrolle, es 
daß für 20 oder Mittelwert 9, —=0 werden soll, so wird [vgl. F.H.425, S.10, G1.(35)] wege 1. 
i BR Be 


De 
Ask 


7 1 E 


en . 2x 7 I 
3 2) _— | m —— ER 1 ze 

Mi N 2:9, dz Pe +6 co hier | co ©’ N 
-3#+2#) Dr, DER 


Le: 
vr) 2 


also ist 


Z 
9. Pe Fi 6 (a 
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Um von dieser Annahme unabhängig zu sein, wird die endgültige Bestimmung von (©, erst in 
Ziffer 4 durchgeführt und vorerst mit G1.(6) gearbeitet. 


3. Die Lösung der Differentialgliederung 
Die Lösung von Gl.(5) wird in der Form 9—= 9, +9, gesucht, wobei 9,,„ die allgemeine 
Lösung der zu G1.(5) gehörenden homogenen Gleichung Ä 
BO OEN 020 
Pe dx == 922 A A super Le (7) 
E ms 
bedeutet. Der Ansatz 9,» —Cy(z)e ”*, führt auf die Gleichung 
IL SE RS x, (8) 
deren ‘vollständige Lösung durch 
9 E 
y— zil2 Ez 3 ma) +@%J_1 5 mar) a En N (9) 
gegeben ist [vgl. F.H.425, S.7, Gl.(22)]. . 
| Die Randbedingungen fordern dy/dz—=0 für z=0 und z=1; durch Differentiation von 
7 G1.(7) folgt unter Beachtung der für Besselsche Funktionen geltenden Beziehungen 
ni d 2 2 3 
E = mi ade 6 mer) Jaıs 3 „or)| NR OR N ehe (10). 
4 
2 
j Nun ist J3/,(0)—0, aber J_g7 (0)—=%; daher kann das Produkt 2J_s 3 m2>l?| für z=0 
einen endlichen Wert haben. Das ist auch der Fall, wie aus der für Besselsche Funktionen 
allgemein gültigen Entwicklung 
ar 22 ie) 
,O-STCHT ehr Aa 2) (Ar-LAy ni i 
| mit = Email hervorgeht, es folgt 
F 14 TE 
x MC J_ Emm) = onr20 
also eine Konstante; daher muß 0, —= 0 sein. Für z= 1 folgt aus G1.(10): 
: dp 2 
h E—masnlgm)=0. 
> Diese Gleichung ist erfüllt, wenn m— > m, gesetzt wird, wobei mit m, die Nullstellen. der 
e Funktion Js73(x) bezeichnet werden. Die Lösung der homogenen Gleichung wird also 
2 
n ,) 2 


Drım nn B3 C, „112 J_1j3 (m,, 2312) e E 
a1 ; 


und die allgemeine Lösung der G1.(), 9= 9, -+ Bym - 
A Bi 3 2 © 
err ri t oo — 5 Mu ve, 
Nee ea nt SO, unmone ee 
Bi E j ; Pe 6 al i 
2 ? \ 4. Die Erfüllung der Anfangsbedingung 
Zur Erfüllung der Bedingung 9—=0 für = 0 muß 
Dr ner. - * © 
DR: «; BEN ERT_ lm) ee. e. (12) 
f } sei 


raus c, und die C,, zu bestimmen sind. Zu diesem Zweck wird die linke Seite der Gl.(12) 
nach Besselschen Funktionen fortschreitende Reihe entwickelt (vgl. F. H.425, S.8.). 
als o diese Gleichung mit 2 und der Eigenfunktion, also mit 2°? J_17 (m; 2?) für 


. ar, RE 


-i 
Dur 
% 


h 
2 

\ 
Ei 
u 


a 
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i—=1,2,...© zu multiplizieren und über 2 von O bis 1 zu integrieren. Dabei treten folgende 


2/3 1/3 
Integrale auf, die mit m; 2? =o, 2= 2 ae 2 (1) do auf die lineare Ver- 
mM; 3M; mM; 


änderliche im Argument der Besselschen Funktion umgeformt werden: 


10 . m; 

[3 J_173 (m;2?/?) de = 2 les-ın (0)de= M, 
' mi 
1 


mM: 


4 ı 
[: 292 J_175 (m;2?®) d2= ErTeIT [ow J_ıs(eo)de=N; 
() u) 
1 
s 6, 2/2 J_17 (m; 2") d2 = 60, |??? Ja) (m; 2)], — 0 ES SER, 
ö 
1 
’ fs C;2?J_13 (m, 2) I-ı3 (m; @P)d2e=0 für u#t ; 
0 
1 m; 
4 ; 
[617-1 map ae 5 | EW-in Ode = 2 F-An mol 
Ur 
0 


0 


Damit ergibt sich für (,, 

M„—N, 
>71: (m,)P 
Die Integrale M,, und N, lassen sich nicht analytisch ausdrücken, sie wurden durch numerische 
Quadratur ermittelt (vgl. Ziffer 5). 

Zur Bestimmung der Konstanten c, muß Gl.(12) mit z und der zur ersten Nullstelle m, —=0 
gehörigen Eigenfunktion y,—=1 multipliziert und dann integriert werden, wie bei P. Schaf- 
heitlin, Theorie der Besselschen Funktionen, Leipzig und Berlin, S. 73/74 näher ausgeführt ist. 

Bei der Integration der mit 2 multiplizierten G1.(12) fällt die rechte Seite fort, weil alle 


Integrale vom Typ der Gl.(13) sind. Die Auswertung der linken Seite liefert c,—7/60; das ist 
derselbe Wert, der durch die Überlegung der Mittelwertbildung beim Partikularintegral unter 


q,= RANG 


Ziffer 2 gefunden wurde. Außerdem wird für die Mittelwertbildung und die Bestimmung von c, _ 
in der Entwicklung in die Fourier-Bessel-Reihe dasselbe Integral verwendet, obwohl es sich 


um zwei ganz verschiedene Vorgänge handelt. 
Jetzt sind alle verfügbaren Konstanten bis auf die zahlenmäßige Berechnung bestimmt, 


die Lösung erfüllt die Differentialgleichung und die Rand- und Anfangsbedingungen, die Funk- 


tionenfolge konvergiert gleichmäßig gegen die gesuchte Lösung. Näheres darüber vgl. F. H.425, 
S.9. | 


5. Die zahlenmäßige Auswertung 


Für die Berechnung der Lösung wurden die Tafeln von Watson, Bessel-Functions, _ 


Oxford 1922, benutzt, die dort tabulierte Funktion Y,/s mußte auf J_ı,; umgerechnet werden. 
Die erste Nullstelle von Js); ist von Emde aufgegeben, Z. math. Physik (3) Bd. 24 (1916) 


S. 239/250 mit m, = 3,375 6106. Die weiteren Nullstellen konnten mit den üblichen Näherungs- 
formeln auf fünf Dezimalen bestimmt werden. Ebenso wurden die Werte von Jı7s(m,) durch 
Interpolation berechnet. Für die Fourier-Bessel-Entwicklung erwiesen sich die ersten fünf 


Tabelle 1. Zahlenwerte für die fünf ersten Nullstellen 


u 1 2 3 4a Bi 
my 3,3755 6,5303 9,6766 12,8206 | 15,9637 
Mu —0,15 072 | 0,013310 | —0,009 296 | 0,002 7802 | —0,002 5301 | 
Ny —0,10123 | 0,019 503 | —0,007 301 | 0,003 6142| —0,002089 
J_us(mu)| —0,43 260 | 0,31189 | —0,25 636 0,22 277, 120, 1DDeR a 
Al 0,1322 0,03 183 0,01 518 0,008 402 | 0,0055 
($ m,)* 25,637 95,949 210,682 369,828 IR FRTECNE 


£ 
E 
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Glieder wegen der stark zunehmenden Argumente der Exponentialfunktion als ausreichend. 
Die wichtigsten Zahlenwerte sind in Tabelle 1 zusammengefaßt. Der Vollständigkeit halber soll 
auch noch die Lösung wiedergegeben werden: 


_22 ,1[(7 = —, 25,64 2 
I=-5.+7 5 324 22) — 0,1322 YzJ_113 (3,376 Y$) ©” Po 


_0,03183 Ya J_us (6,530 Ya) e "Pe 

— 0,01518 YzJ_1s (9,677 VD) Pe 2. (BB). 
— 0,008402 Y2 I_415 (12,82 /) or 

_. 0,005529 Ya I_41s (15,96 Ya) e * Pe.. 


03% Die Ergebnisse der Zahlenrechnung sind 
in Bild 2und 3in Parameterdarstellung wieder- 
gegeben. Für die in F. H.425 gerechneten Bei- 
spiele mit U = 20, 40 und 60m/s., A= 0,Imm 
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feste Wand 
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Bild 2 Bild 3 


Ba — 10-:m und einer Spaltlänge von O,lx cm, entsprechend dem halben Umfang eines Zapfens 


von 200 mm Dinr. sind in Bild 4, 5 und 6 die wirklichen Temperaturen ermittelt und perspek- 
_ tivisch wiedergegeben. Die Pe&cletschen Zahlen sind 25920, 51840 und 77760; sie gehen hier 
als Zahlenwerte ein; wegen des fehlenden Wärmeüberganges kommt ihre physikalische Bedeutung 
_ nur bedingt zur Geltung. : | 

Bemerkenswert ist für den Vorgang, daß das von der bewegten Fläche mitgenommene 
-Öl bis ins Unendliche hinein eine niedrigere Temperatur an der bewegten Wand mit kon- 
er Differenz gegenüber der festen Wand verursacht. In Bild 5 und 6 stellen die ge- 
"Kurven in der letzten Temperaturverteilung rechts diesen Verlauf dar, der durch 
ularintegral nach dem zweiten Glied der rechten Seite, von Gl.(6a) gegeben ist. 
t der Unterschied gegenüber dem bereits erreichten Temperaturverlauf so gering, 
t mehr in Erscheinung tritt. 
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6. Die mittleren Temperaturen 
verte für die Temperaturverteilung (d) in strömenden Flüssigkeiten mit der Ge- 
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gebildet. Nach Einführung der dimensionslosen Größen y bzw. z und d ergibt sich hieraus 
für 9 die schon in Ziffer 2 genannte Formel: 


1 
9—=2|20 da. 
0 


Wird Gl.(11) mit e,— 7/60 mit 2 z multipliziert und von 0 bis 1 integriert, so fallen alle Glieder 
mit Ausnahme von 2x/Pe weg, weil die Integrale mit J_,,s die Form der Gl.(13) haben und 
6, so bestimmt war, daß die Glieder mit 2? und 2? ausgeglichen werden. Die Mittelwerte an 
der Stelle x sind also 
= 2% -  nÜU: 2(2)n.Ü 
9=5, oder Glen Nm). . 
Das ist aber der gleiche Wert, den man aus Gl.(1) erhält, wenn der Einfluß der Wärmeleitung, 
also das Glied mit A vernachlässigt und für (w) der Mittelwert U/2 eingesetzt wird. Die dimen- 


sionslose Form für 9 folgt entsprechend aus Gl.(5) mit z= 1/2. 

Damit ist gezeigt, daß für die Couette-Strömung und eine Flüssigkeit, deren Zähigkeit 
nicht von der Temperatur abhängt, für die Berechnung der Mittelwerte der Temperatur die 
Strömungsgeschwindigkeiten quer zur Schicht gemittelt werden dürfen; ebenso darf dazu die 
Wärmeleitung quer zur Schicht vernachlässigt werden. j 

Offen ist noch, ob diese Beziehungen bei anderen Geschwindigkeitsverteilungen, also bei 
der Hagen-Poiseuille-Strömung und bei mit der Temperatur veränderlicher Zähigkeit auch 
gültig bleiben. Diese Frage muß eine besondere Untersuchung entscheiden. In wirklichen 
Schmierschichten bestimmt aber der Couette-Anteil weitgehend den Charakter der Strömung, 
ebenso wird die veränderliche Zähigkeit in erster Linie die Zahlenwerte hinsichtlich der Höhe 
der erreichten Temperaturen beeinflussen, so daß an dem grundsätzlichen Ergebnis kaum etwas 
geändert werden dürfte. 

Für die Durchführung der Zahlenrechnungen bin ich Herrn Dipl.-Physiker U. Rost und 
Frau E. Rost zu besonderem Dank verpflichtet. . 


. (16). 


Eingegangen am 11. Dezember 1950. 


Über die Fehlerabschätzung beim Extrapolationsverfahren 
von Adams | 


I. Gleichungen 1. Ordnung 
Von P. Matthieu in Zürich 


Es wird eine neue Methode der Fehlerabschätzung entwickelt, die mit größter Leichtigkeit im Anschluß 
an die erfolgte Integration einer Gleichung vorgenommen werden kann. Die Fehlerschranke beträgt im 
allgemeinen etwa das Anderthalbfache des wirklichen Fehlers. Die Methode wird für die beiden wichtigsten 
Formen des Adamsschen Verfahrens zur Integration von Gleichungen erster Ordnung näher ausgeführt. 


. Anew method of estimating the error is presented that can be very easily performed as soon as Ihe equation ” 
E is integrated. The bounds of the error are generally about 1,5 times larger than Ihe real error. T’he method 
is explained in Ihe detail for ihe two most important forms of the method of Adams used for the integration 
of differential equations of Ihe first order. 
Une methode nowvelle pour Vestimation d’erreur est developp£e, applicable avec grande facilite, apr&s 
que V’equation fut iniegree. La limite d’erreur est en general le sescuple, & peu pres, de l’erreur reel. 
methode est exposee d’une manitre detaillte pour les deux formes les plus importantes du procede d’ Adams 
d’integration d’equations du premier ordre. 
Paspa6arsıaerch HOBBIÜÄ METON OMCHKH NOTPEIIHOCTH, KOTOPbIH MOKEeT ÖbITp € DoNBmoh 
JIETKOCTBIO IIPHMEHEH TIOCAIe HHTETPUPOBAHNA ypasHeHua. Ilperest morpemmocTn OÖBIYHO HPH- 
ÖJIHSNTENBHOB Ba Pasa TIPeBEIMIAET MeÜCTBHTEIBHYIO norpemmoctp. Ilpen saraemsıä meToQ 
us1araetcı bomee TOAPOOHO MIA ABYX OCHOBHEIX hopm nHTerpmpoBauna AuhepenmmansHbıx Be, 
ypasHeHüt mepBoro mOopsaRa 10 cnocoßy Anamca. Wer: 


1. Einleitung Bu“ 

Es ist heute nicht mehr bestritten, daß unter allen Methoden zur angenäherten numerischen 
Lösung gewöhnlicher Differentialgleichungen das Extrapolationsverfahren von Adams das weit 
aus bequemste und vorteilhafteste ist. Das Verfahren kann) dabei in verschiedenen Gestal 
durchgeführt werden. Über diese ist in der Literatur so häufig berichtet worden, daß hier n 
mehr näher darauf eingegangen werden soll!). Neuerdings spielt das Verfahren von Ad 


1) Vgl.z.B. E.Kamke: Differentialgleichungen, 2.Aufl. Leipzig 1943, S.143—148; oder Fr.A. { 
Methoden der praktischen Analysis, 2. Aufl. Berlin 1950, 8. 348-361. BE ee” we: 
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wiederum eine wichtige Rolle, weil es diejenige Methode ist, nach welcher Differentialgleichungen 
vorzugsweise vermittelst der großen digitalen Rechenmaschinen gelöst werden. Auch die vor- 
liegenden Untersuchungen können unmittelbar und sogar besonders gut auf den Fall der maschi- 
nellen Integration angewandt werden, weil die dabei verwendeten Intervalle meist sehr klein sind. 
Wie bei jeder Näherungsmethode bildet naturgemäß auch beim Extrapolationsverfahren 
die Frage der Fehlerschätzung ein wichtiges Problem. Dieses ist erst verhältnismäßig spät in 
Angriff genommen worden. Erst seit etwa zwei Jahrzehnten ist von verschiedenen berufenen 
Autoren eine Reihe von Arbeiten über diesen Gegenstand veröffentlicht worden, die großenteils 
' In dieser Zeitschrift erschienen sind?). Indessen ist zu sagen, daß die in diesen Arbeiten ent- 
wickelten Methoden, so wertvoll sie auch sind, bei der praktischen Durchführung doch große 
Schwierigkeiten bereiten. Außerdem ist der wirkliche Fehler oft vie] kleiner als die durch die 
Abschätzung gegebene Schranke. Aus diesen Gründen soll hier eine andere Methode der Fehler- 
abschätzung entwickelt werden, die sich in diesen Hinsichten günstiger verhält. Diese ist ohne 
wesentlichen Unterschied auf alle Differentialgleichungen anwendbar. In der vorliegenden Arbeit 
beschränken wir uns jedoch auf den Fall der Gleichungen erster Ordnung, während diejenigen 
höherer Ordnung in einer späteren Arbeit behandelt werden sollen. 


2. Prinzip der Lösung 


Dem hier verfolgten Gedankengange liegt ein allgemeines Prinzip zugrunde, das nicht in 
unmittelbarem Zusammenhange mit dem Extrapolationsverfahren steht, das sich jedoch in 


) besonders sinngemäßer Weise auf das letztere anwenden läßt. Es sei eine reelle Differential- 
| gleichung erster Ordnung in der Form 
gegeben und außerdem als Anfangsbedingung eine Beziehung 
U U EN Te 


wobei z,und y,gegebene Zahlwerte sind. Gesucht sei die durch (1) und (2) bestimmte Lösung y(#) 
für 2 x,. Von der Funktion f(x,%) setzen wir zunächst nur voraus, daß sie im Rechteck 


u SE Smt u TE RE REN.: 
stetig und nach % stetig differenzierbar sei. Ferner gelte innerhalb des Rechtecks (3): 
ir € DEE RN a a a FREE 


undes bezeichne % die kleinere der beiden Größen a und b/K. Dann verläuft y(«) für u <e <utk 
sicher innerhalb des Dreiecks 


VHSEe<mtk My Kl al N ua 
3 Im Unterschied zur exakten Lösung y(x) bezeichne y*(x) eine irgendwie erhaltene Näherungs- 
F lösung, von der wir nur voraussetzen, daß sie der Anfangsbedingung (2) genüge, innerhalb des 
Bo %: Dreiecks (5) verlaufe und stetig sowie stückweise stetig differenzierbar sei. Setzen wir ferner 
ru wobei ö(z) den unbekannten Fehler der Näherung y*(x) bezeichnet, so muß (6) der Gl.(1) 
) genügen, was auf die Beziebung 
5 IE N Le SR a u) 
führt. Der Fehler ö(x) genügt dabei gemäß den gemachten Voraussetzungen über y*(x) der 
Beziehung 
| a a LS 


- Die rechte Seite von (7)entwickeln wir nach Taylor und erhalten bei Verwendung der Lagrange- 
_ schen Form des Restgliedes in leicht verständlicher Schreibweise 


Es ’ Ä vr + far) +3 my + N ne 


en ap { . . . . x . . 
Pa) ist dabei eine Funktion von x, von der nur bekannt ist, daß ihre Funktionswerte zwischen 
den Grenzen 0 und 1 liegen. Setzen wir zur Abkürzung 


er Nr ° | : ° i 
= mr), na=—% (nyt+ 9). d=fl&yN)—y" (10), 
=) Es handelt sich vor allem um folgende Arbeiten: Mises, R. v.: Zur numerischen Integration von 
ntialgleichungen. Z. ang. Math. Mech., Bd. 10. — Schulz, G.: Interpolationsverfahren zur numerischen 

on gewöhnlicher Differentialgleichungen. Z. ang. Math. Mech., Bd. 12. — Tollmien, W.: Über die 
schätzung beim Adamsschen Verfahren zur Integration gewöhnlicherDifferentialgleichungen. Z. ang. 
1., Bd. 18. — Fricke, A.: Über die Fehlerabschätzung des Adamsschen Verfahrens zur Integration 
er Ditferentialgleichungen erster Ordnung. Z. ang. Math. Mech., Bd. 29. 
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so schreibt sich (9) in der Form f 
m = ee ee 


Das ist eine lineare Differentialgleichung für 6, gemäß welcher ö unter Berücksichtigung von (8) 
in der Gestalt | 


er] f pıda 
ö(a)=e ” RE a 


dargestellt werden kann. Hierbei ist jedoch zu bemerken, daß der Koeffizient p,(x) nicht genau 
bekannt ist, weil in ihm die unbekannten Funktionen ö(x) und 9(x) auftreten. Hingegen läßt 
sich (12) leicht abschätzen. Besteht im Dreieck (5) die Uneleichugz 

m<— (ay<M. a a 
so ergibt sich gemäß (12) für u, <x <x,+ k sofort | \ 


T . 
In sem [Iglet de ee, AL. 


Falls weiterhin im obigen Intervall die Ungleichung 


| 
USD RN NN re © 


öl < Zar men ea 


besteht, folgt aus (14): 


Die Ungleichungen (14) und (16) lassen sich in ganz ähnlicher Weise praktisch verwerten wie die 
Abschätzungsformeln in den unter ?)genannten Arbeiten. In vielen speziellen Fällen lassen sie sich 


noch verschärfen und vereinfachen. Falls z.B. m>0, also 2 a y) im Dreieck (5) überall negativ 
ist, ist der Ausdruck 


SE Spla)dz 
a | e* N 


eine monoton wachsende Funktion,’und aus (12) ergibt sich sofort 


S % 5 
|ö(«)| </ lg] de < Lama)... een (18). 
{ Ganz ähnlich ist im Falle M <0 (17) eine monoton abnehmende Funktion, für die die Beziehung 
a 
P [mad Br 
> 4 < e” B= em (2—%) Be a A Tu : EP fe (19) 3 X E 
besteht. Daraus folgt gemäß (12) sofort >= 
{ | al sermen gde<i@—a)emtm no. 2 Or 
2 E7 
' . Alle vorstehenden Formeln beziehen sich auf das Intervall u <x <z,+k. ! - 


Die eigentliche praktische Bedeutung dieser Betrachtungen beruht nun aber darauf, daß 
i in allen konkreten Fällen ö(«) stets so klein ist, daß man in den Formeln (11) und (12) pı(@) 
he durch die bekannte Funktion p(x) ersetzen kann. Das kommt darauf heraus, daß man in der 
Taylorschen Entwicklung von (7) nur noch daslineare Glied in ö berücksichtigt, während man 
die höheren Glieder vernachlässigt. Auch kann aus der exakten Darstellung (12) stets sehr 
leicht abgeschätzt werden, welche Auswirkung die Ersetzung von 9, durch p auf die Größe a iR 
Fehlers hat. Dazu sind zwei Dinge erforderlich. Erstens ist die ungefähre Größe des Fehlers zu 
schätzen; zweitens ist abzuschätzen, welchen Einfluß die Ersetzung von p, due in (12) 
hat. Damit erhalten wir nun an Stelle von (11) die Gleichung PAR 


"pn A 
Wir bezeichnen sie im folgenden als die zu (1) und y* (2) gehörende Fohlerglichung, Au 
ergibt sich der Fehler ö(x) außerordentlich AnEeRau EN in der Korn u 


— 


KR 
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3. Anwendung auf das Verfahren von Adams bei unmittelbarer Extrapolation 


Bei der Anwendung auf das Verfahren von Adams ist zu beachten, daß y*(x) sich aus 
lauter ganzen rationalen Funktionen der gleichen Ordnung zusammensetzt. Für jedes Intervall 
wird ja eine neue derartige Funktion gebildet. An den bisherigen Überlegungen kann dies jedoch 
nichts ändern. Wir machen vorerst über die Funktion f(z,y) noch die zusätzliche Annahme, 
daß sie sowie ihre sämtlichen partiellen Ableitungen bis zur 6. Ordnung stetig und nach jeder 
Veränderlichen stetig differenzierbar seien®). Die weiteren Entwicklungen sind nun je nach der 
verwendeten Gestalt des Adamsschen Verfahrens verschieden. 

Wir beginnen mit dem Fall der unmittelbaren Extrapolation und entwickeln die sich 
ergebenden Methoden am Beispiel der sehr häufig verwendeten Fxtrapolation vermittelst Poly- 
nomen vom vierten Grade. Zur Durchführung des Verfahrens ist es zunächst nötig, das Inte- 


‚grationsintervall, ausgehend von 2£—=x,, in eine Anzahl gleicher Teilintervalle einzuteilen. Die 


Begrenzungspunkte seien mit %,, &, 23, ..... bezeichnet und die gemeinsame Intervallänge mit h. 
Ferner ist es, damit das eigentliche Verfahren in Gang gebracht werden kann, nötig, die ange- 
näherte Lösung y*(x) für die drei ersten Intervalle, insbesondere für die vier Werte &,, %}, X, 23 
auf irgendeine Weise möglichst genau vorauszubestimmen. Wir bezeichnen diese Funktionswerte 
durch y*(2,)=yn (0 <n<3; ya —Y,) sowie ferner die sich daraus ergebenden Größen f(x,, %%) 
durch fz. Dann stellt man in bekannter Weise!) das Differenzenschema: : 


Log Y fs 


” Fr * Aft 2% * 
1 Yı fi if f? wert 
L rg 3 = 
%g Y5 Er ir d’fz Fe Aifs 5 
a» u KA If} 73 S A 5 ach - 
Afi d’fs 


2 yi Fi Af3 


af: 


auf®). Die Größen y; und fr haben dabei für n>3 die gleiche Bedeutung wie für n<3. Der 


Extrapolationsprozeß besteht darin, daß man, nachdem y*(x) bereits in den ersten n—1 Inter- 
vallen (n>3) bekannt ist, zunächst mit 


Un—1 = —n De DE ER De ke VER ar Ute (24) 
das Interpolationspolynom 
A er 4° = n 
Br (z) N ER + er : Un—1 H f : Un—1 (un-ı + 1) | 
; a . (25) 


73 


ze A we Un—1 (in-ı + 1) (Un—ı u. | 


bildet. Daraus ergibt sich dann die im n-ten Intervall zur Verwendung gelangende Näherungs- 
funktion Y7(x) in der Form 


Ya)=ytat | Frıla)da ea RT 


n—1 


aus der man weiterhin den nächsten Tabellenwert 7 sofort in der Gestalt 


M=lia)=ytat | Fuala)de | 
n—1 i AL 5 } (27) 
= y-MiHth BE aA * Br ras 5% g Are (n >53) | 


erhält. Die vollständige Näherungsfunktion y*(z) erscheint also in dieser Weise zusammen- 
gesetzt aus der für die drei ersten Intervalle vorausbestimmten Anfangslösung und für n >3 
aus den einzelnen biquadratischen Polynomen (26). Bei der Durchführung des Verfahrens ist 


‘es nicht nötig, die genaue Gestalt (26) dieser letzteren zu kennen. Es genügt vielmehr, die 
Werte y% direkt gemäß (27) zu berechnen. Für unsere Zwecke ist es jedoch gerade wichtig, den 


3) Über die Änderungen, die eintreten, wenn diese Bedingung nicht erfüllt ist, vgl. S. 362. t 
- 4) Die vierten Differenzen sind zur Integration von (1) nicht nötig, werden jedoch bei der Fehlerabschät- 
gebraucht. 


r 
- 


sich zeigen, daß unter sehr allgemeinen Voraussetzungen |g(z)| in jedem einzelnen Intervall 
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Gesamtverlauf der Funktionen (26) in den zugehörigen Intervallen in Betracht zu ziehen. Dieser 
ergibt sich nach durchgeführter Integration der Gl.(1) aus dem Schema (23) gemäß (25) und (26) 
sofort in der Form 


Ya(a) = ya + h|fanı Un +6 Afsa +34 "a1 +24 a) nn 
| (28). 
en an (n>3) 


Nachdem nunmehr der Gesamtverlauf der Ra. y*(x) auf Grund von (28) bekannt 
ist, kann der Fehler ö(x) gemäß (22) berechnet werden. Es ist grundsätzlich festzuhalten, daß 
dies mit beliebiger Genauigkeit geschehen kann, wenn man nur genügend viele Punkte der 
Kurven (28) verwendet. Für die praktische Anwendung handelt es sich jedoch darum, solche 
Rechenmethoden anzugeben, die möglichst einfach im Gebrauch sind. Wir machen dazu noch 
die in. konkreten Fällen immer leicht zu erfüllende Annahme (von der man sich übrigens nötigen- 
falls ohne jede Schwierigkeit befreien kann), daß die in den drei ersten Intervallen angenommene 
Ausgangsfunktion dort nicht merklich von der genauen Lösung y(x) abweiche. 

Da die beiden Exponentialfaktoren in (22) stets positiv sind, erhält man eine Schranke 
für den Fehlerbetrag |ö(z)| wenn man in (22) q(x) durch irgendeine Majorante des Betrages la()| 
ersetzt. Im folgenden wird gezeigt, wie eine slche Majorante stets in einfachster Weise gebildet 
werden kann, womit dann gleichzeitig auch unser Problem der Fehlerabschätzung im wesentlichen 
gelöst ist. Zunächst lassen sich in jedem Intervall zwei Werte der Funktion g(x) sofort angeben, 
nämlich die Werte in den beiden Intervallendpunkten. Gemäß (10), (26) und (28) hat g(8) i im 
n-ten Intervall (n >3) die Gestalt 


1) = Sa Vila) — FR) BR, 
= (m, Ya(@)) flo yo) - yo) | | 
Weiterhin folgen aus (26) und (27) die Gleichungen: 
Yale ner Yalan) = ya. Sue DS 
sowie aus (25) die Beziehungen 
Y% (&n-ı) — Fa (In—ı) =, 21 \ 
Ya) ra) =ftat+ Atari) 
Durch Einsetzen von (30) und (31) in (29) erhält man schließlich 
(1) a Ya-1) — fr = (re a (32) 


1(@.) = f (ün, Ya) —fr-ı — ie "a fin “oe Men = En ne ). 
Während also die Funktion g(z) nach der eben gemachten Voraussetzung in den drei ersten 
Intervallen identisch verschwindet, besitzt sie nachher einen unstetigen, zackenförmigen Verlauf, 
indem sie im »-ten Intervall zwischen den Werten 


Biete! (ma) 0, la) der lee . (34) 
variiert. In Bild 1 ist die Funktion g (x) für ein Beispiel ausgerechnet. Die zrindickende (34) 
können sofort angegeben werden, da sich die vierten Differenzen unmittelbar aus dem Schema (23) 
ablesen lassen. 

Es lassen sich nun aber noch weitere Schlüsse ziehen über den Verlauf von g(«). Es läßt 


& 


al Zt ab m ZN ana 


. (81). 


und 


erstens monoton ansteigt und zweitens nach oben konkav verläuft, daß also die beiden ersten 
Ableitungen von |q(z)| stets positiv sind (vgl. Bild 1). Damit können wir aber eine sehr einfache 
Majorante g(x) für |q(«)| bilden. Diese ergibt sich dadurch, daß man die einzelnen Kurven- 
stücke, aus denen |g(z)| besteht, durch die Sehne zwischen ihren Endpunkten ersetzt. Die auf 
diese Weise entstehende, aus lauter geraden Strecken bestehende Funktion (x) ist, ebenso 
wie |g(@)| in den Intervallendpunkten unstetig. Sie ist in Bild 1 gestrichelt eingetragen. ‚Ge- 
mäß iR kann 9(x) äußerst einfach berechnet werden. Nach der obigen Bemerkung hat man 
dann also eine Fehlerabschätzung in der Form 


I —_ Fade Tpak 
| <iln)= ee,” ge wnde 


Als untere Grenze kann bei allen Integralen in (35) der Wert x, gewählt werden, da nach ı 
gemachlen Voraussetzung ö(x) in den drei ersten Intervallen ijenliseh verschwindet 
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Es handelt sich nun vor allem noch darum zu untersuchen, unter welchen Bedingungen |g (=) 

in jedem einzelnen Intervall monoton ansteigt und nach oben konkav verläuft. Dazu entwickeln 
wir g(&) im n-ten Intervall um den Anfangspunkt nach Potenzen der kleinen Größe h, was nach der 
auf S. 359 gemachten Voraussetzung sicher möglich ist bis zur 5. Potenz. Zur Aufstellung dieser 
Reihe nehmen wir zunächst an, daß y* (x) für &,_ı <x® <x„_ıgenau mit der exakten Lösung (x) 
übereinstimme. Ferner bezeichnen wir die Werte der Ableitungen dieser letzteren Funktion 


Sense a duich.y, 29a, 2: De ist 
Y (2) = 57 mn an 1! 1 Ym—ı Un—1 Aa 9 ER, mi Eli 3! are, wm 4. sur (36). 
Weiterhin ergibt sich die Entwicklung von Y% (2) = Fn-ı(«) nach (25). Unter der gemachten 
Voraussetzung bestehen die Gleichungen: 2 
SAm9R Rn. ,.; 2ER, 
JhAa= Fin Yn = Yn—ı— T! Yn—1 > ar YUd1 mi 1 Yyıt --: 
2h „ ZZ sh? un 
ing == Yu = Yn_ı — II Yn—1 fg 9 BE 1 Yn—1 +... CR Ih (37), 
% ‘ ?, h „ vr h3 777 
f me ln Un 1 71 Yn—1 ER 9] YA, = 1 Ya se 
) Jn-ı — Yn-1 
r aus denen sich sofort für die in (25) auftretenden Differenzen die Entwicklungen 
- Ye RR? >. h3 yar h* 
Afn-ı —— hin Gr 57 Yn—1 + 6 Ya 1 Dr3 Yn—ı + = A 
42 un ZZ 7h* h? 
Afti-— Bypa— hy 51 ae T ren 
A vn 3ht 553 
A — h? Uni FOR Yn-ı + — 4 Ynmı ee 


ergeben. Setzt man diese a in (25) ein, so erhält man durch leichte Rechnung 
Ir) HABE, nsE 11 Dr Er a" Fe uni Fa oa ae un 
z 


Du ; 5) 
TER) | Br 1 Ynı ni Un—1 oe mit 54 2 


22 Die Reihen (36) bis (39) sowie auch die im folgenden zu benutzenden Reihen setzen sich, wenn 
dr J(z, y) eine analytische Funktion ist, ins Unendliche fort. Andernfalls sind sie in bekannter Weise 
) durch ein Restglied abgeschlossen zudenken. Nach der Voraussetzung auf S. 359 kann angenommen 
ne werden, daß dies frühestens nach der 5. Potenz der Fallist. Der eckige Klammerausdruck in (29) 
BR ni: ergibt sich dann als Differenz von (36) und (39) ii in der Form 


Er -[Y% (@) — y'(#)] = ht yn_ı [% Uni 2 54 er mit > 54 wa) | 


. (40). 
+1 B ER Va Beat: ne... 

: Ya 10 y—1l — 9 N— 5) 24 N— 120 

; Dr R ee ergibt sich durch Integration von (40), wenn wie bisher die Glieder bis zur 5. Ordnung 


ke angeschrieben werden, 


1 11 
[LYR(a) — y(«)] = Sn — ya E Un + 72) Un + = Un Haar 120 NE DEN 41). 


andelt sich nun noch darum, die Differenz der beiden ersten Terme in der zweiten Zeile in (29) 
immen. Man erhält zunächst 


Sf Yia)) )—f(» ya) = en 5 ua) EL) ZIIE " 2. (A). 
A r ‚Abkürzung 


anne), 


Bi: 
E: 
Ba 
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und schreibt wiederum nur die Terme bis zur 5. Ordnung, so nimmt (42) gemäß (41) die Form an 


fl. Y3a)) —f(a ya) 
23 Te 110 1 1 
—=—h Yn—ı ,, = Un—1 + 72 Un—1 + 16 Un—1 + 120 Un—1 + 


Damit ergibt sich nun auf Grund von I; und (43) für Ri im n-ten Intervall die Entwicklung 


1 A 
q(&) = h*ı ls NE: 54 ER y Un— ir 5 Un— ) 


of (43). 


r 


3 5 1 
— h Ds ie Un—1 Rn 5 n-ı — 94 Un — 10 ui) (44). 


v [of 1 15% 1 1 
She en \ 2" tt ging “) es 
a Det 


Das ist die gesuchte Reihenentwicklung für q(x). Da h eine kleine Größe ist, wird der Verlauf 
von q(x) wesentlich bestimmt durch das erste, h* enthaltende Glied in (44). Wenn wir zunächst 
die höheren Glieder vernachlässigen (ihr Einfluß wird nächstens noch genauer untersucht werden), 
so gilt für q(x) unter Berücksichtigung von EN die sehr angenäherte Darstellung 


3 1 4 
a a = A Hal Un—1 +4 z Um 2 Una J4 u) | ER ER a0n (45). (> 
In Sn DEmas 1 | 


Diese Funktion ist aber ihrem Betrage nach stets monoton ansteigend und nach oben konkav. 
Für die Herleitung der Näherung (45) wurde über die Funktion y*(x) die Voraussetzung 
gemacht, daß sie für m 4 <® <x„_ı genau mit (x) übereinstimme. Dies ist zunächst nach 
Voraussetzung für n—4 der Fall. Infolgedessen gilt im 4. Intervall für g(z) die genaue Dar- 
stellung (44). Aber schon im nächsten Intervall ist die Voraussetzung nicht mehr erfüllt. Es 
soll nun untersucht werden, was diese Ungenauigkeit für einen Einfluß auf den Verlauf von g(&) 
im 5. Intervall hat. Wir verstehen dazu für diese Betrachtung (und nur für diese) im 5. Intervall 
unter y(x) nicht wie bisher die genaue Lösung des Problems, vielmehr diejenige Lösung von (1), 
die an Stelle von (2) durch die Anfangsbedingung %(2,) = yi bestimmt wird, was zulässig ist, 
wenn wir q(x) gemäß (29) bilden. Im Unterschied zum 4. Intervall stimmt dann aber y(x) in 
den drei vorangehenden Intervallendpunkten nicht mehr mit den Werten y} überein, unter- ’ 
scheidet sich vielmehr von diesen um Beträge 5. Ordnung, und das gleiche gilt auch für die 3 
Funktion y'(x) = f(x, y) hinsichtlich der Werte ff. Bei den Größen des Schemas (23) sowie bei 
den sämtlichen Formeln (37) bis (44) treten somit für n—5 Änderungen 5. Ordnung ein. Genau 
die gleichen Schlüsse lassen sich auch für höhere Werte von n ziehen, wobei wieder die Definition 
von %(x) entsprechend abzuändern ist. Insbesondere werden also in (44) für n>4 noch weitere 
Glieder 5. und höherer Ordnung auftreten. Ihre Gestalt wird stark mit den einzelnen 
Extrapolationsschritten wechseln, jedoch ergibt sich aus der Natur des Verfahrens, daß ihre 
Bedeutung gegenüber derjenigen der Glieder 4. Ordnung im allgemeinen stark zurücktreten wird. 


Die genannte Voraussetzung kommt also, ebenso wie die Ersetzung von (44) durch (45) 
auf die Vernachlässigung höherer Terme in den Entwicklungen für 1(«) heraus. Dies kann auf 
unsere Betrachtungen höchstens insofern einen Einfluß haben als evtl. in ganz extremen Fällen 


der für die Herleitung von ö(x) als wesentlich vorausgesetztemonoton ansteigende und nachoben 
konkave Verlauf von g(x) geändert werden kann. Die gleiche Bemerkung gilt auch hinsichtlich 
der auf 5.359 gemachten Annahme über die Differenzierbarkeit der Funktion f(z2,y). Auf der 
anderen Seite gibt es wieder zahlreiche Gründe dafür, daß g(x) auch dann eine positive erste und 
zweite Ableitung hat, wenn die Glieder 5. Ordnung stark vertreten sind. Z.B. sind diese letzteren 
rationale Funktionen von ähnlicher Art wie die Glieder 4. Ordnung. Weiterhin genügt q(&) 
mit großer Annäherung den in späteren Fällen exakt geltenden Gleichungen (66), in denen 
jedoch n durch n—1 und die Differenz durch A*f% zu ersetzen ist. Schließlich kann g(e), ie 
aus der Herleitung von (44) hervorgeht, auch bei Berücksichtigung der höheren Glieder mit 
äußerster Annäherung als ganze rationale Funktion 4. Grades betrachtet werden und hat 
also stets einen recht einfachen Verlauf. Zusammenfassend kann also gesagt werden, daß g(@) 
in jedem Intervall mit äußerster Wahrscheinlichkeit stets monoton ansteigend und nach oben 
konkav verläuft. AS 
Wenn auch eine Änderung dieses Sachverhaltes in praktischen Fällen kaum jemals zu w 
erwarten ist, soistes doch prinzipiell wertvoll, eine Kontrollmöglichkeit über den Verlauf von na) BT 


ash Bin u 
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zu haben, wodurch dann das ganze Verfahren auf eine vollständig strenge Grundlage gestellt 

wird. Eine solche Kontrolle kann aber in der einfachsten Weise jederzeit gemäß (29) erfolgen, 

indem man g(%) für einige geeignete Zwischenwerte berechnet. Dabei werden stets ganz wenige 

solche Werte in einigen wenigen geeigneten Intervallen genügen, auf Grund von welchen dann 

der Verlauf von g(&) sofort abgeschätzt werden kann. Z.B. wird es oft vorteilhat sein, den 

Wert ET ee 2) in der Mitte des n. Intervalles zu bestimmen, wobei zur Abkürzung 
» 2 . 2 


CT ı=%—— gesetzt ist. Dazu müssen gemäß (29) zuerst die Werte 


Re 2 
- EN ea, RE (A 
2) | es 
gebildet werden. Der erstere ergibt sich gemäß (28) aus dem Schema (23) in der Form 
se I 1 AARTEk- 
tn u are | rad 
der zweite gemäß (26) und (25) in der Gestalt 
1 - 3-2 er | 
vn le — ©; Afnı+t 54 fr in 16 Aafzsı EEE NE U Ber (48). 
) Damit kann man nun sofort gemäß (29) die Größe 
Sl * EAFL 
ER ER v4) y_ı a EEE A 


berechnet werden. Andere Werte von q(x) ergeben sich genau entsprechend, jedoch dürfte die 
Anwendung von (49) auch in den ungünstigsten Fällen genügen. 


4. Praktische Durchführung der Fehlerabschätzung in graphischer und 
| tabellarischer Form 


Die entwickelte Fehlerabschätzung wird im wesentlichen durch Formel (35) gegeben, in 
welcher alle Funktionen der rechten Seite leicht berechnet oder abgeschätzt werden können. 
Man kann diese Formel in verschiedener Weise verwenden. Der einfachste Weg ist der graphische, 
der zuerst kurz besprochen werden soll. Da die Funktion 9(x) ihrer zackenförmigen Gestalt 
wegen nicht besonders geeignet ist, ersetzen wir sie nochmals und zwar durch eine geeignete 
mittlere Majorante 7(x), d.h. durch eine möglichst regelmäßig verlaufende Kurve, für die die 
Beziehung 


oa lıaa........2...., (50) 


besteht. Eine solche Majorante 9(x) kann stets mit größter Leichtigkeit angegeben werden (vgl. 
ö) “ Bild 1). 5(x) kann nun wieder in (22) an Stelle von q(z) eingesetzt werden, wodurch man für, 


den Fehler eine’ Majorante ö(«) erhält. Es ergibt sich somit die Fehlerabschätzung 

= — /pds JS pds 

r lö(z)| <ö(x) = e * Is ED. ea 
R z : &s 

$ Die Schranke (2) kann in einfachster Weise nach bekannten graphischen Methoden berechnet 
RUE oder abgeschätzt werden. Man vergleiche dazu das nachfolgende Beispiel. Noch einfacher kann 
_ man zum Ziel kommen, wenn man g(x) an Stelle von g(x) in die Fehlergleichung (21) einsetzt, 

_ das Richtungsfeld dieser Gleichung durch Berechnung einiger Werte bestimmt und dann die 


Kurve ö(x) näherungsweise aus dem Richtungsfeld abliest. 


Zur tabellarischen Durchführung der Fehlerabschätzung greifen wir auf Formel (35) zurück. 
Schreibt man neben das eigentliche Rechenschema (23) ein zweites Schema der Gestalt 


er, 
W 


a: y Qz 


b 04 d, e4 ga 
[77 
& Y b, Gz d, e g BR u EI ae (52) 
a 
> b; [en ds & Is 
173 


ka gear 20 m al Eee BE Ba WI er a ne se 
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dessen Größen definiert sind durch 


1%, ,* ih. 
= 57 dy (%3, %5) » Se ur (8 ) n>3 | 
n—1 b 
RL WERE (53), 
i=3 
2 n l 
dn = c,| AR] ’ 0, = a di, I = 
i=4 2 Cn 
so erkennt man leicht, daß mit außerordentlicher Annäherung die Gleichungen bestehen 
2 © 
Der n ©, 
| Ste > / pda 
- %a Fr - 
He ; RT 1(8) e dx, öl), aloe 
%z 
Für den Fehler ö(@,) am Ende des »-ten Intervalles gilt also gemäß (35) die Abschätzung 
|) < Alan) = 94. 2-20 ee re 


Eine kontinuierliche Fehlerschranke erhält man, indem man die einzelnen Werte g, durch gerade 
Linien verbindet oder durch eine Kurve interpoliert. Die Größen (53) sind außerordentlich einfach 
zu berechnen, speziell auch für Rechenmaschinen. 

Damit ist nun unser Ziel erreicht. Zusammenfassend kann man sich an die folgende Regel 
halten: Die Differentialgleichung (1) sei unter der Anfangsbedingung (2) nach Adams integriert. 
Von der Funktion f(x,y) wird dabei einzig vorausgesetzt, daß sie in einem Bereiche, der nur die 
Näherung y*(x) sowie die exakte Lösung y(x) im Inneren enthalten soll, stetig und stetig nach y 
diffenrenzierbar sei, was durch Abschätzung in jedem konkreten Fall sofort festgestellt werden kann. 
Ferner ist abzuschätzen, ob für den Fehler die Näherung (22) zulässig ist. was ebenfalls in elementarer 
Weise geschehen kann. Dann erhält man eine Schranke für den Fehlerbetrag |ö(x)|, wenn man in (22) 
an Stelle von q(x) irgendeine Majorante für |g(&)| einsetzt. 

Praktisch wird man als Majorante immer die Funktion q(x) verwenden, was auf die außer- 
ordentlich leicht zu berechnenden Abschätzungen (51) und (55) führt. Die Verwendung von q(&) 
beruht auf der Vernachlässigung höherer Glieder in gewissen Reihenentwicklungen für q(x). Obwohl 
sich kaum ein Fall denken läßt, in dem diese Vernachlässigungen nicht zulässig nn besteht die 
Möglichkeit, ihren Einfluß in einfachster Weise zu kontrollieren und nötigenfalls an Stelle von 1(%) 
eine andere Majorante einzuführen®). Die entwickelte Methode stellt also ein ebenso einfaches wie 
strenges Verfahren der Fehlerabschätzung dar. 

Diese Betrachtungen lassen sich noch in interessanter Weise weiterführen. Wenn die 
Majorante y(x) so gewählt wird, daß in (50) annähernd das Gleichheitszeichen steht, was sich 
ebenfalls stets leicht erreichen läßt, so ergibt sich aus (45), daß, wenn q(x) bzw.die Differenzen A'fs 
im ganzen Intervall immer das gleiche Vorzeichen haben ‘(was fast immer der Fall sein wird), 
auch der Fehler stets dieses gleiche Vorzeichen hat. Ferner muß sich in diesem Falle die Fehler 
schranke (51) zum Betrag des een Saas stets angenähert verhalten wie 


0: De I a Ei 
2 tar * 360 (56). = 
55: 09 Re 


l Re +, u) 


Die Fehlerschranke ist N im genannten Falle gleichen Vorzeichens im allgemeinen (d.h., wenn 
wieder die in der obigen Regel ausgesprochenen Vernachlässigungen zulässig sind) stets etwa 
das Anderthalbfache des wirklichen Fehlers. Das gleiche Verhältnis ergibt sich stets automatisch Fe > 
für die Abschätzung (55). Voraussetzung bei den sämtlichen obigen Betrachtungen ist, daß das ri 
Verfahren mit so vielen Stellen durchgeführt wurde, daß die Abrund unge Rai gegen IR 
über den übrigen Fehlern nicht merklich in Betracht kommen. Br. we 


5. Berechnung eines Beispiels) Bu 


Die Gleichung y’—=2—y mit der Anfangsbedingung y(0)—0 soll im Intervall. 
nach Adams integriert werden und es soll der dabei begangene Fehler ab 


°) Praktisch kann man sich wohl auch in den ungünstigsten Fällen auf die Betrachtu 
gemäß (47) — (49) beschränken, von denen bei rechnerischer Durchführung einige als 


6) Es a absichtlich eine ganz einfache Differential ARHABE gewählt, damit 
der exakten Lösung sofort nachkontrolliert werden kann. je 


Z. angew. Math. Mech. 
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Die Lösung kann in der Umgebung des Anfangspunktes sehr rasch und genau durch 
Potenzreihen berechnet werden. Man errechnet auf diese Weise die Werte: 


y(— 0,1) = 0,0051709 
y(0,1) = 0,0048374 


Wir wählen dementsprechend 2, = — 0,1 und = 0,1. Um mit der bisherigen Bezeichnungs- 
weise im Einklang zu bleiben, muß man sich dann nur die ursprüngliche Anfangsbedingung durch 
die ihr völlig äquivalente erste 
Gleichung (57) ersetzt denken. Das 


y (0,0) = 0,0000000 | 


y (0,2) = 0,0187308 | °  ° ° © ° 57). 


lee 
Schema (23) nimmt dann die e- mt 0353| 0000 IL Be 
stalt (58) an, wodurch die gegebene N a) 


Gleichung integriert wird. Ferner 
wurde, obwohl dies zur Durchfüh- 
rung der Fehlerabschätzung nicht 
notwendig ist, der Verlauf der zu- 
gehörigen Funktion |g(z)| ermit- 
telt. Das Resultat ist in Bild 1 auf- 


0,570* 


X = 1,0 


Ko -01 X 0,2 


| gezeichnet. In den drei ersten EN 
Intervallen verschwindet diese : 
) Funktion identisch, in den folgenden erkennt man dagegen deutlich den monoton ansteigen- 


den Charakter von |g(x)| und seiner Ableitung. Der nicht ganz regelmäßige Verlauf der Zacken 


ist auf den Einfluß der Glieder höherer Ordnung zurückzuführen. 


%, yi fa Afk Arfa fr Acfn 
—0,1 0,0051709 —-0,1051709 
0,1051709 
0,0 0,0000000 0,0000000 —0,0100083 
0,0951626 0,0009523 
0,1 0,0048374 0,0951626 ; —0,0090560 —0,0000935 
0,0861066 0,0008588 
0,2 0,0187308 0,1812692 —0,0081972 —0,0000750 
0,0779094 0,0007838 
0,3 0,0408214 0,2591786 —0,0074134 —0,0000792 
0,0704960 0,0007046 
0,4 0,0703254 0,3296746 —0,0067088 —0,0000654 
0,0637872 0,0003692 (58) 
0,5 0,1065382 0,3934618 —0,0060696 —0,0000625 
. | 0,0577176 0,0005767 
0,6 0,1488206 0,4511794 —0,0054929 —0,0000529 
0,0522247 0,0005238 
0,7 0,1965959 0,5034041 —0,0049691 —0,0000516 
2 0,0472556 0,0004722 
) 0,8 0,2493403 0,5506597 —0,0044969 —0,0000441 
nu - b t 0,0427587 0,0004281 
. 0,9 0,3065816 0,5934184 —10,0040688 
0,0386899 
f 1,0 0.3678917 0,6321083 


ware 


4 


Do Genauer Wert: 4(1,0) = 0,36787944 


a "Anschließend an die Integration kann nun der dabei begangene Fehler abgeschätzt werden. 
Dies soll zunächst in der auf S. 363 erwähnten graphischen Form geschehen. Der Gang dieser 

einfachen Rechnung ist unmittelbar aus Bild 2 zu ersehen. Zunächst sind die Werte Id4f „| (in 
_  10:mal vergrößertem Maßstab) als Punkte eingetragen. Weiterhin isteine zugehörige F unktion 4(@) 
soeingetragen, daß sie gemäß der Bemerkung auf S. 363 möglichst etwa in halber Höhe dieser 


ınkte verläuft. Damit kann nun der Ausdruck (51) für (x) gebildet werden, wie das in Bild 2 
ittweise auf graphischem’ Wege durchgeführt ist, wobei der Übersicht wegen einige Zwischen- 
_ weggelassen sind, die man vorteilhafterweise auch graphisch aufgetragen wird. Für 
bt sich speziell ö(1)= 1,85 - 10°, während aus dem Rechenschema (58) für den 
en Fehler der Wert ö(1)—=— 1,23 - 10° hervorgeht. Der Fehler hat also, wie auf 
rkt, tatsächlich das gleiche Vorzeichen wie die Differenzen Are; und die ge- 


hranke ist etwa das Anderthalbfache des Fehlerbetrages. 
} 25 
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Zur Durchführung der Fehlerabschätzung in tabellarischer Form ist das Schema (52) 
gemäß (53) auszurechnen. Als Resultat ergibt sich: 


0,5 

1 0,05 1,05 0,0000982 0,0000982 0,0000468 
1,0 

} 0,15 1,16 0,0000870 0,0001852 0,0000798 
1, 

0,25 1,28 0,0001015 0,0002867 . 0,000112 
1, 


E: 0,35 1,41 0,0000916 0,0003783 0,000134 - 
» A ieh 
0,45 1,57 0,0000981 0,0004764 0,000152 Dr 


1,0 

m 0,55 1,73 0,0000915 0,0005679 0,000164 

Ar 0,65 1,92 0,0000990 0,0006669 0,000173 
N 0,75 2,12 0,0000933 0,0007602 0,000180 

1, 


Gemäß (55) stellt die letzte Kolonne eine Ma- 
jorante für den Fehler dar. Speziell erhält man 
für 2—=1 die Schranke g,1= 1,80 - 10”. Diese 
ist also wieder etwa das Anderthalbfache des Ü 
zugehörigen Fehlerbetrages |ö(1)| = 1,23 - 10%. 
Schließlich kann man sich noch, falls man dies 
für nötig hält, sowohl bei der graphischen wie 
bei der tabellarischen Form über die Zulässigkeit 
der Annahmen vergewissern. Z.B. können ge- 
mäß (49) einige der Größen qg ,„ berechnet 
ee 
werden. Die resultierenden Werte sind aus Bild 1 
zu ersehen. 


6. Fehlerabsehätzung bei mittelbarer 


“a . Extrapolation 


Die bisherigen Entwicklungen können fat 
unverändert auch auf diesen Fall übertragen 
werden. Das Verfahren von Adams werde 
x7 70” wieder durchgeführt vermittelst Polynomen vom 

Bild 2 4. Grade, wobei wir auch hier annehmen kön- 
nen, daß die angenäherten Werte in den drei 
ersten Intervallen nicht merklich von den genauen abweichen. Der Unterschied gegenüber der 
unmittelbaren Extrapolation besteht darin, dann daß für n >3 der Wert (27) nicht als endgültiger 
sondern nur als erster Näherungswert betrachtet wird, der noch iterativ verbessert wird?). 
(27) sei jetzt dementsprechend durch %7,o bezeichnet sowie der zugehörige Wert f(z,, Yn,0) 
durch fx,o und die Differenzen durch A af}, 0A ee 0s =: +. Alle diese Größen werden wieder 
in das Differenzenschema (23) eingetragen. Die iterative Verbesserung besteht darin, daß man 
vermittelst der obigen Größen den neuen, besseren Näherungswert 


X 01 %=0,2 


1 - 1 - 1 .- Pi 
Ynı —— Yn,0 En h If 0 D7 Afn,o — 'D AR — 94 nu BR ae. EAGBR (60) 


bildet, der in gleicher Weise weiter verbessert wird bis sich die neu entstehenden Werte nicht. = 
mehr ändern. Diese endgültige Näherung, die der weiteren Rechnung zugrunde gelegt wird, sei 
wie früher durch yı bezeichnet sowie die zugehörigen Größen im Schema (23) durch fx, A Ara FAR 5 
Auch die iterativen Verbesserungen bis zu den endgültigen Näherungen werden im Schema 3) 
eingetragen. in 
Zum Näherungswert 4% gehört wiederum eine Näherungsfunktion Yz(z) mit y: (= 
welche letztere ihrerseits gemäß (29) den Verlauf von g(x) im n. Intervall bestimmt. 
lassen sich zunächst wieder entsprechend (34) einige exakte Aussagen machen. N ähr a; 


?) Vgl. die unter !) genannten Werke sowie die unter ?) an dritter Stelle genannte 


AR 
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G1.(30) und die erste Gl.(31) ihre Gültigkeit unverändert beibehalten, geht die zweite G1.(31) 
der Natur des Iterationsverfahrens entsprechend®) über in 


| Yn (&) = Fr (&,) SE le ih). 
so daß an Stelle von (34) die Gleichungen treten 
EN 0 Ur Be 2 REED 


g(«) verläuft jetzt also stetig (vgl. Bild 3), während g’(&) in den Intervallendpunkten unstetig ist. 

Die Funktion Y7(z) ist eine ganze rationale Funktion 4. Grades, von der bis jetzt 
immer nur der Teil innerhalb des n.Intervalles in Betracht gezogen wurde. Im folgenden soll 
diese Einschränkung aufgehoben und insbesondere der Verlauf von Y}(x) im n+1. Intervall 
betrachtet werden. Wir behaupten, daß dann die Beziehung besteht 


le ER RR ER 


wobei Y7+1,0(«) entsprechend der obigen Bezeichnungsweise die zumersten Näherungswert X 1,0 

gehörende Näherungsfunktion für das n+1. Intervall bedeutet. In der Tat folgt (63) wieder 

unmittelbar aus der Natur des Iterationsverfahrens, indem die Ableitungen der beiden Polynome 

in (63), d.h. in unserer Bezeichnungs- 9 

weise die beiden Polynome 3. Grades 

Fn-ı(2) und F,„,o(x) durch dieselben 
’) vier Punkte bestimmt werden. 


n. Intervall 


Gleichzeitig mit Y}(z) kann natür- 
lich gemäß (29) auch der für das n. Inter- 
' vall definierte Teil von g(x) über dieses 
Intervall hinaus fortgesetzt werden. Wir 
bezeichnen die auf diese Weise ent- 
stehende Funktion durch Q,(x) (Bild 3). 
q(x) entsteht dann also aus den Funk- 
tionen Q,(x) dadurch, daß man von je- 
dem:Q,(z) nur den innerhalb des n. Intervalles verlaufenden Teil in Betracht zieht. Aus (63) 
folgt ferner sofort 


An (X) 


1x) 9X) 1X) 
Bild 3 


AS N N N: 


Nun ergibt sich, wie soeben bemerkt, die Funktion Y7;ı1,0(&) aus den früher berechneten 
Werten des Differenzenschemas in genau gleicher Weise wie dies im Falle der unmittelbaren 
Extrapolation für die dort betrachtete Funktion Y};ı(x) der Fall war. Folglich können wir 
auch die Eigenschaften dieser letzteren Funktion auf den jetzigen Fall übertragen. Insbesondere 
gilt gemäß (34): 


On+1,0(&n41) = A: fa0so TR LE EN EN SR ee (65). 
Demnach genügt Q,(x) den drei Gleichungen: 
) an), Mn), man) Afürno >: > >: (66). 


. Die ersten beiden dieser Gleichungen sind direkte Wiederholungen von (62), während sich die 
‚dritte gemäß (64) aus (65) ergibt. 2 

Neben den exakten Gl.(66) lassen sich wie früher noch weitere Schlüsse ziehen über den 
Verlauf von g(z) bzw. Q,(x) im n. Intervall. Gemäß (64) und (44) besitzt Q,(x) um den End- 
punkt des n. Intervalles die Entwicklung 


a Wär ee 19,4 
= hr (E are u)t- nr TEE SER (67). 


Berücksichtigen wir in dieser wie früher nur das angegebene erste Glied, so erhält man unter 


ER 

E 4, Beachtung von (66) für q(z) im n. Intervall die (45) entsprechende Näherung 

ar & 1 KR a 

® a) = Afkino (4 ur Eee! «i) Dee (68). 
. HAIE, 1<w<0 


Damit ist unser Ziel im wesentlichen erreicht. Nachdem der Verlauf von g(x) in den ein- 
zelnen Intervallen gemäß (68) wieder in einfacher Weise außerordentlich angenähert aus dem 
_  Differenzenschema bestimmt werden kann, kann man wieder eine Majorante für |q(x)| sowie 

eine zugehörige mittlere Majorante bestimmen, dergestalt daß die zugehörigen Ausdrücke (35) 
_ und (51) Schranken für den Fehlerbetrag darstellen. Noch einfacher ist es jedoch, den jetzigen 
I KR 


a EEE a as; 


Ya 
; 
k- 
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Fall auf den früheren zurückzuführen, indem die Integrale über die beiden Klammerausdrücke 
in (68) und (45) im n-ten Intervall sich im Mittel zueinander verhalten wie 


0 
ra 11 1 1 
2 2 3 Fri 4 
IG+- tet) 19 1 
2 ER LLUNN ee 


in Tr 
n 1 1 241 132105 


1 
er en N 
Ir ur 54% Het u) au 


Dazu definieren wir auch im jetzigen Fall eine Funktion 9(z). Sie soll, ähnlich wie früher, da- 
durch charakterisiert sein, daß sie für die drei ersten Intervalle identisch verschwindet, während 
sie für n>3 am Anfang des n-ten Intervalles verschwindet und von da an in gerader Linie zum 


’ II 
Wert >41 Men 
jetzt 9(x) keine absolute Majorante mehr für |g(z)|. Weiterhin kann dann die unstetige Funktion 


4(x) wieder durch eine stetige mittlere Majorante y(x) ersetzt werden. Zusammenfassend gelangt 
man dann zu folgender Aussage: 


am Ende des n-ten Intervalles ansteigt. Im Gegensatz zu früher ist aber 


Die auf 8.364 angegebene Regel kann in allen Stücken auch auf den jetzigen Fall der mittelbaren 
Extrapolation übertragen werden. Dabei sind einzig die Differenzen A'f% durch die Größen 


19 — - Ar BE Dr 
BT A’fkı,o (angenähert — S 2 fe) zu ersetzen, bzw. es ist die Definition von G(x) in der ü 
eben angegebenen Weise abzuändern. 


Auch die Möglichkeiten der strengen Kontrolle können ohne alle Schwierigkeiten über- 
nommen werden, wie das jedoch nicht im einzelnen ausgeführt werden soll. Ebenso kann die 
Rechnung natürlich wieder graphisch oder tabellarisch durchgeführt werden. Da gemäß der 
obigen Regel der frühere Index n durch n+1 zu ersetzen ist, ist im jetzigen Fall eine Differenz 
mehr zu berechnen, falls man sich nicht mit einer Extrapolation begnügen will, wie das praktisch 
immer ausreichen wird. Schließlich kann man auch jetzt schließen, daß, wenn die Differenzen 
Afeo immer das gleiche Vorzeichen haben, die Fehlerschranke stets etwa das anderthalb- 
fache des Betrages des wirklichen Fehlers aus- 
macht. Jedoch ist das Vorzeichen des letzteren 
jetzt demjenigen der Differenzen entgegengesetzt 
entsprechend dem negativen Vorzeichen in (69). 

Erwähnt sei ferner noch die folgende in- 
teressante Aussage, die sich aus den obigen Aus- 
führungen bei Annahme der zugrunde gelegte 
Näherungen unmittelbar ergibt: 

Das Verfahren von Adams ist bei mittel- 
barer Extrapolation stets.rd. I3mal genauer als bei 
unmittelbarer. 


7. Beispiel zur Fehlerabsehätzung bei 
mittelbarer Extrapolation 


Das gleiche Beispiel, das unter 5. durch 
unmittelbare Extrapolation berechnet wurde, soll 
jetzt unter Verwendung der mittelbaren berech- 
net werden. 

Wir verwenden die gleichen Ausgangs- 
REN RO 410 * werte. Dann nimmt das Di hen die 
Gestalt (70) an. Die endgültigen Näherungen 
ergeben sich dabei durch zweimalige Iteration®). 
Die letzte Kolonne stellt die Differenzen A* fx,. dar, die nun unmittelbar zur Fehlerabschätzung. Fa 
verwendet werden können. Diese kann z.B. wieder graphisch gemäß Bild 4 erfolgen, welch 


letzteres genau Bild 2 entspricht. Für @=1 erhält man speziell ö(1)=1,3-10, während dr 
wirkliche Fehler gemäß (70) den Wert ö( 1)=0,8 10” hat. Der Fehler hat also, wie das nach 


Bild 4 


®) Esist dabei nicht nötig, alle ARE SEBR Differenzen zu berechn ie d in (70 d Üb ae 
halber geschehenist. Der Rechnungssang läßt sich vi lmeh tli ü Ei N 4 u g 
Methoden der praktischen en 3. Aufl. San, 350352, ee) Br ne Val an, fa MAN 


i 
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ya fa: Afki Aa WER Af%,0 
—0,1 0,0051709 —0,1051709 
0,1051709 
0,0 0,0000000.  0,0000000 —0,0100083 i 
0,0951626 0,0009523 
0,1 0,0048374 0,0951626 —0,0090560 —0,0000935 
0,0861066 0,0008588 
0,2 0,0187308 0,1812692 —0,0081972 
0,0779094 0,0008623 
0,3 0,0408214 0,2591786 —0,0081937 
0,0779129 0,0008622 
0,0408179 0,2591821 —0,0081938 —0,0000860 
0,0779128 0,0007762 
0,0408180 0,2591820 —0,0074176 
0,0704952 0,0007795 
0,4 0,0703228 0,3296772 —0,0074143 
0,0704985 0,0007793 
0,0703195 0,3296805 —0,0074145 —0,0000763 
0,0704983 0,0007030 
0,0703197 0,3296803 —0,0067115 
0,0637868 0,0007059 
0,5 0,1065329 0,3934671 —0,0067086 
0,0637897 0,0007058 
: 0,1065300 0,3934700 —0,0067087 —0,0000701 
. ) 0,0637896 0,0006357 
0,1065301 0,3934699  : —0,0060730 
0,0577166 0,0006383 
0,6 0,1488135 0,4511865 —0,0060704 
0,0577192 0,0006382 (70) 
0,1488109 0,4511891 —0,0060705 —0,0000627 
0,0577191 0,0005755 
0,1488110 0,4511890 —0,0054950 
0,0522241 0,0005779 
0,7 0,1965869 0,5034131 —0,0054926 
0,0522265 0,0005778 
0,1965845 0,5034155 —0,0054927 —0,0000572 
0,0522264 0,0005206 
0,1965846 0,5034154 —0,0049721 
0,0472543 0,0005228 
0,8 0,2493303 0,5506697 —0,0049699 
0,0472565 0,0005227 
 0,2493281 0,5506719 —0,0049700 —0,0000516 
0,0472564 0,0004711 
0,2493282 0,5506718 —0,0044989 
0,0427575 0,0004730 
0,9 0,3065707 0,5934293 —0,0044970 
0,0427594 0,0004729 
0,3065688 0,5934312 i —0,0044971 —0,0000465 
1 0,0427595 0,0004264 
0,3065688 0,5934311 —0,0040707 
0,0386886 0,0004282 
) 1,0 0,3678803 0,6321197 —0,0040689. 
0,0386904 
0,3678785 0,6321215 
0,3678786 0,6321214 


0 a 2 


By 


Genauer Wert: y(1,0) = 0,36787944. 


8. Schlußbemerkungen und Zusammenfassung 


u En; der vorstehend entwickelten Methode der Fehlerabschätzung werden gegenüber den 
ren Arbeiten zwei neue Gedanken benutzt. Erstens wird zur Bestimmung des Fehlers die 
leichung (21) bzw. ihre Lösung (22) verwendet, wobei man praktisch stets mit Leichtigkeit 

zen Bo ob diese Danerung zulässig ist. Zweitens wird die wesentliche Funktion 22) 


unseren Überlegungen sein muß, tatsächlich das entgegengesetzte Vorzeichen wie die Differenzen 
A 7 TE o, und die Schranke ist wieder rd. das Anderthalbfache des Fehlerbetrages. 
, findet man auch die auf S.28 ausgesprochene Regel bestätigt, wonach der Fehler bei mittelbarer 
Extrapolation rd. 13mal kleiner ist als bei unmittelbarer. 


Schließlich 
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für den Fehler erhält. Bei unmittelbarer Extrapolation lassen sich solche Majoranten sofort 
angeben auf Grund der vier folgenden Eigenschaften von g(&): 

1. g(z) genügt den Gl.(34); 

2. |g(«)| ist in jedem Intervall in erster Näherung monoton ansteigend und nach oben 
konkav; 

3. g(x&) läßt sich gemäß (29) in einfacher Weise berechnen oder abschätzen; 

4. q(x) kann, auch bei Berücksichtigung von höheren Gliedern, mit äußerster Annäherung 
als ganze rationale Funktion vierten Grades betrachtet werden. 

Praktisch wird man als Majorante immer die Funktion 9(x) verwenden können, die sich 
in ganz elementarer Weise aus dem Differenzenschema ablesen läßt. Man gelangt so zu den 
außerordentlich leicht zu berechnenden Äbschätzungen (51) und (55). 

Bei mittelbarer Extrapolation kommt man durch ganz analoge Überlegungen zum Ziel, 
wobei sich noch als bemerkenswertes Resultat ergibt, daß sich dieser scheinbar viel kompliziertere 
Fall in einfacher Weise auf den ersten zurückführen läßt. 


Eingegangen am 22, Dezember 1950. 


Über die konforme Abbildung gewisser rechtwinkliger 
Achtecke 


Von N. Bundscherer in Landshut (Bayern) 


Durch Überlagerung zweier mit Einschnitten versehener Rechtecke, deren Abbildung auf einfache Funk- 
tionen zurückgeführt werden kann, gelangt man zu gewissen rechtwinkligen Achtecken. Die Parameterwerte 
des Abbildungsintegrals lassen sich aus den Achtecksseiten berechnen. Spiegelung dieser (speziell auch 
ausgearteter) Achtecke liefert eine Anzahl abbildbarer Bereiche. 


By swperposition of two notched rectangles that can be represented by relatively simple functions, the 
author gives Ihe conformal representation of cerlain rectangular octagons. The parameter values of Ihe 
integral that gives Ihe represenlalion are computed from the sides of the oclagon. By continued reflection 
of the oclagons (especially of degenerale ones), a great number of representable regions is generated. 


Par superposition de deux rectangles enlailles, dont la reprösentation conforme se fait au moyen de 
fonctions simples, on peul repr£senter cerlains octogones reciangulaires. Les valeurs parametriques de 
l’integrale de representation sont caleules & l’aide des cötes d’octogone. Par reflewion de ces oclogones, aussi 
d’octogones degeneres, on oblient quelgues domaines repr£ösentables. 

Hano:kenue AByX IIPAMOYTOABHMKOB C HanpesaMmH, KOHPOPMHOE OTOÖPakeHHe KOTOPEIX 
UPOHSBONHTCH IPH MOMOIM IIPOCTEX PYHEUHH, MAeT HEKOTOPBIe IPAMOYIOABHLIE BOCBMNUY- 
TOABHHKY. SHAYeCHNUA TAPAMETPOB OTOÖPpasKamıeroO HHTErPanma OIpenensiotca u8 CTOPOH 
BOCBMHYTOJIBHURA. AHalmTmyeckoe IPONOMKEHHe ITUX BOCBMHUYTOJIBHHKOB AaAeT TEeNIBIH PAA 
KOHPOPMHO OTOÖPpaskaeMmblIx OÖNAacTeh. 


Einleitung 


Ist & das Innere eines endlichen, geradlinigen, geschlossenen Polygons der 2-Ebene, das 
in den Ecken 4A, #—=1,2,..., n) die Innenwinkel a, besitzt, so bildet nach Schwarz- 
Christoffel ein bestimmter Zweig der Funktion: 


1 


n 
:(w)=0) Ho@—a)®7 do .12.07. 2 
v=1 
Wa 
wo Ü und (” geeignet gewählte Konstanten sind, die obere w-Halbebene $ konform auf & ab!). 
C und ©’ bewirken nur eine Ähnlichkeitstransformation, sind also unwesentliche Konstanten. 
Dagegen besteht in der Auffindung der den Punkten A, in der w-Ebene Ense „Ver- 
zweigungspunkte «a, die eigentliche Schwierigkeit. 
Da keine allgemeine Methode die Berechnung der Verzweigungspunkte aus der Gestalt 
eines vorgegebenen Polygons ermöglicht, ist die Lösung dieses ‚Parameterproblems‘‘ wenigstens 


für spezielle Klassen von Polygonen anzustreben. Setzen wir voraus, daß die Innenwinkel &, x nur 
4 r 


3 N 
die Werte 5 und Di annehmen, so erhalten wir ‚‚Treppenpolygone‘; das Abbildungsintegral 


der Gl.(1)ist in diesem Fall vom hyperelliptischen Typ. S. Bergmann konnte die Halbebene 9 
auf das Innere einer gewissen Klasse von Treppensechsecken abbilden .) [1]. In der vorliegenden 


Arbeit wird durch eine ähnliche uetboge die Abbildung von 9 auf eine Klasse von Treppen- er 


1) Siehe z. B. [6], S. 218. 


?) Unter „‚Abbildung‘‘ soll stets konforme Abbildung verstanden werden. ' 
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achtecken durchgeführt. Zwischen den Seitenlängen dieser Achtecke bestehen fünf Gleichungen, 

so daß noch drei Seiten (oder Gleichungen für die Seiten) beliebig vorgeschrieben werden können. 

_ Für diese Achtecke läßt sich das Parameterproblem lösen. Die Abbildung z(w) kann durch 

schrittweise Transformation mittels leicht übersehbarer Funktionen erreicht werden; dadurch 

wird die praktische Schwierigkeit der numerischen Auswertung des hyperelliptischen Abbildungs- 
integrals wesentlich vermindert. 


$1. Vorbemerkungen 
1. Die bekannte Transformation ?) 


y2 
=u+— 
TR 


bildet die obere w-Halbebene so auf die volle «-Ebene ab, daß die reelle Achse in den von 2r über 
oonach —2r verlaufenden reellen Verzweigungsschnitt der x-Ebene, der Halbkreis w=r ei? (mit 
0 <sp<n)in die Strecke —2r <x<2r übergeht. Zwei auf den beiden Ufern des Schnittes 
einander gegenüberliegende Punkte sind spiegelbildlich zu dieser Strecke gelegen; daher ent- 
sprechen ihnen nach dem Schwarzschen -Spiegelungsprinzip zwei (reelle) Punkte w, und ,, 
die der Bedingung, 
N ee 3 ee ee 17 (2) 
genügen. 
) 2. Durch einen bestimmten Zweig der Funktion 


(mit O<k<1) 


NE do 
et y eb A RE 
| Ya—r) 1 —Ra}) 


‚wird die obere w-Halbebene auf einen Rechtecksbereich der &-Ebene abgebildet*). Den Modul %, 
den mit ihm durch die Gleichung k®+%’?—=1 verbundenen konjugierten Modul % und die 


Periodizitätsmoduln 
1 1/k 
I j 2 und ıK'= RL 200 22 VATER 
JYa—?) (1—k?t2) J ya—t2) (1—K2t) 
9 i 


[G 


| :K 
können wir mit Hilfe der 9-Funktionen durch das (rein imaginäre) Periodenverhältnis 7 = - 


ausdrücken. Dabei ist 0 < G < ©?°). (Da hier ausschließlich die 9-Nullwerte verwendet wer- 


? i den, geben wir nur das Periodenverhältnis als Argument an.) Es ist ®) 
je | 9%(e) 3) 
er ked)=—5; ra IE TEUER EN DEN Be 3); 
Si. | re) ae) o 
) | Koa-iane; Kd=—iräne) 2...) 
2 wobei 
Nr Ben; Po F) j oo © 
BD ri dee. 
fr Nn= —ox Nn=—» Nn=—» 


und n- Br 0 ist. Ferner-gilt?): 
Se RR, (1) —iröse) TER WATTE (6a), 


A 3-1) ir) ER rin RS (CHE 
| »l-)=- 70:0) Be. NG; 


Sa) Vol. (31,8. 9088. . 
4) Vgl. 2. B. (#], S. 296ff. } 
4 5 & —0 bzw. — =, so artet das Rechteck in einen unendlich langen Parallelstreifen aus. 
B. [4], 8.202 und 14. 
4, S. 236, Gl. (5°) mit vo= 0. 


r 


y u 
” 


wählen. Es ist 
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also wegen (3) und (4): na | | Be B 
ter. when] 0 
1 1 EN AR al: ” { 4 
K Ki 2) Kr), K’ (& 2) —=Kß) 
, I 
Da = reell vorausgesetzt wird, geht für — i7 = nee X nach Null; wegen Ö)wrd E 


als, —0, u, =l, = k=0, 
Das Periodenverhältnis ist eine eindeutige Funktion des Moduls; mit der aus (6c) folgenden 


Beziehung: kK’ = - - (— ir) 9 (rd) = 5 0: e n) und 9,(0)—=0 bekommen wir daher 


für k=0: K=-5; Ko; FE— 5; KK=o; kK=0; kK'=0 (8a) 


und, indem wir r durch = ersetzen, gemäß (7): i 2 


fürh=1: K=-w; K=-5; KE=0; KR=0; kK-o; IE—5 (8b). 
Wie sich aus (5) leicht ableiten läßt, ist für jedes 7 
HS 1% d(T)21 ni 


Für die doppelten Periodenverhältnisse gelten die Formeln®): 
DE LHOEE TG 
282) = 0). 


Lösen wir diese Gleichungen nach 92(r) und 9%(7) auf und ersetzen 7 durch > so erhalten wir 


32) = 0) +90), 
(2) = 


$?. Die Durchführung der Transformation z(w) Fge> i 


Die obere Halbebene 9 der w-Ebene soll durch eine analytische Funktion 2(w) auf das Ru 
Innere & eines Treppenachtecks A, B, C, D, E,F,@, H der z-Ebene konform abgebildet werden. r% 
Die (reellen) Verzweigungspunkte der w-Ebene bezeichnen wir mit a...h. 


Ohne Beschränkung der Allgemeinheit können wir RE a ee 
In, a F;: BE a ENDE 
voraussetzen. Die anderen Parameterwerte c, d, e, f, h wollen wir weiter unten noch Zysckna j 


—0<e<d<e<f<O<h<il IRRE Ma RR 12). ; 
Durch die lineare Funktion | Br 
w—1 
w—d 
'  transformieren wir zunächst die w-Ebene auf eine Ehe (vgl. Bild 1). Die Verzweigungspunkte De 
in der v-Ebene seien a’... h’. TERN Be 
Dann ist \ 1 N; BR 
a —(, b’==1, ug d=.00, ... eh ar reg L % 2 


Sodann soll $ durch 


. RN I ie DE RAN 


+3 er x RK “ Te Rn > 


De N \ Ne, a? Men 4 5 . 
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so auf eine volle x,-Ebene abgebildet werden, daß die Punkte a und h, d und e, c und f bzw. 
 b’ und c’, e’ und Ah’, f’ und g’ paarweise zusammenfallen. Nach Gl.(2) ist dies nur möglich, wenn 
die Verzweigungspunkte die Gleichungen erfüllen: 


SSR Rau u ER N SE ER 2 (16a), 
N a NEE A (16b) 
Ri und (wegen ’=1):—=eh’=f’g, was wir mittels (13) und (11) umformen können in 

| £ c—1l e—1iıh-—1 


PER SERIE] . NA ee DIE FESTE AT Pa (16c), 
c—1 f—-1l1l 
| Or Er TR Ne TE SE (16d) 
Be. & Durch diese Gleichungen sollen die noch freien Parameterwerte c, d, e, f, h a 
werden. 
Für r, und r, gilt: 
Bd au c—1 
= 16% . n= nasd Be 1 gar a hte hen An Ah iu a Fan (17). 
ie Lösen wir (16a), (16b) und (16c) nach den Veränderlichen d und e auf, so erhalten wir die drei 
) Bedingungen: 
wi d(i+e) | _dre & 
5 Be; .- Ge — ey ’ = ———- 1m \ h= de a (18), 


die den folgenden Transformationen zugrundeliegen. Wie man durch Einsetzen bestätigt, 
' erfüllen sie auch die Gl. (16d). 
j Die Größen d und e sollen nur der Einschränkung 


T N RENTEN EUER DR A Tg (19) 
T unterworfen sein. Es ist leicht einzusehen, daß dann auch die Ungleichungen (12) gelten. 
n Setzen wir (17) in (14) bzw. (15) ein, so folgt unter Berücksichtigung der Gl.(18) und (13) 
Er < __ w®+de 
Be a 
Br: -. „. _ 4d—e) ®+1—de _ d(d—e) (w—1)®+ (1—de) (w—d)? 
ME d(d—e)v 2 d(d—e) (w—1) (w—d) 
Sg Aut die x,- bzw.. eye üben wir die Transformationen aus: 
ee 1 21(a) , a %—%,(b). 
a er mal) no) 
Wegen z(J)=1+d=1+h und z(d=d-+e wird 
Be: | w—(l-+h)w+h 
BR ' | & w—(d-+e) w w+de 
(w—1) (w—h) 
2 a er A A LTE aa) AT Te SEE ER Tre st Iran 
are (20a) 
; 1—2de+d? | 2d—e—d?e 


a = ng, Seke - ald-e)vs d—e) und 2%,(g) = %(0) ae 


{ „_ Ute 
Sr er 
U ? BE HS] 1-e uw +2] 
4 1-te 


BEUTE er wW—6 
} ze w(w—f) 


= /(e—4) (e—9) , 
sl NV e—1)c—h) 


N A 
Fr 


wi 


IE EFER EINEN: 


re 


Br Dr 
a u 


ED 


_ das wir durch w ausdrücken wollen. Unter Berücksichtigung von (18) erhalten wir nach einer 


LS A 
* r 


5 Ä ) b s ‘ u s Yr f 
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was mit (18): 


e—d | 
= ee en. ; (21a) | 
ergibt, und entsprechend Kine 
En U-ogele—f) , s 
3 Tim (1—f)(e—0) a - 
woraus 
RE Ei... (21b) 


TE 4 und dechab One 


folgt. Wegen (19) ist (1+d)(1—e)=1—de—(e—d)>0, also Fuer: 


offenbar ist auch O<k,<1. 
Die s-Ebene bilden wir durch 


| Yo?) 7 —kı 0?) 


auf eine 2,-Ebene, die t-Ebene durch“, 


t l nt Re 
m Te: 2) | | | © 


auf eine 2,-Ebene ab. Dabei sei 


1—d Bu : 

41 YA-+a)d-Ee).Q „u ee (23a), “ 

% = YLEB) (Fe) nn ee NE (23b) h 

und C eine willkürliche positive Konstante?). Er I 
Zunächst betrachten wir das Integral: F 


3 | ds g 
 NaZs)d—Re) . 


einfachen Rechnung: 


Nele dee) (1—e) w 
(w—d) (w—e) I; a 
si lreder (m Sn ee —h) . | Be 
BE, REEL I Be u 
e e+d?e-+.de? 
Sa Le nr 
TaRp 25 1 dena (vd) wzeor ne 


1a I -WU—E), (we) W—P) 


(1—de)? (w—d) (w—e) 


und 
we —d, N 
2sds= (1—d) (1—e)- a | | en 
Es ist also BR * , BE .. 
| 2sds IC a Pe 
Blauen ans) re GE % 
IE NENNE 
; 1—.de m ; 
G= nn 
ag +a(+e) 
Wir setzen j 


he nn (112.2) u 
Über den Zweig der Wurzel werden wir sogleich verfügen. | 
°) Die Zweckmäßigkeit dieses Ansatzes wird in der Schlußbemerkung dieses P a 


ge %s b e 3 j > 
BAR: , f 
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w- Ebene | v- Ebene 


s?-Ebene | t?-Ebene 


v2 Ei Er En HE Al 4 CG Die 6% 
a. u s-Ebene t-Ebene 


H,\| 


N? 


Az B2 C2 BD; 
2,-Ebene 
z-Ebene 
ua ı 
ne ee, N are En 
w(w— . w(w— 
De IR: W—2ontef,, 
| Ra 1=e ln sr 


Wir verwenden den Zweig der Wurzel, für welchen C, positiv wird: 


' Durch Überlagerung von 2,(w) und 2,(w) bekommen wir endgültig: 


.L Tr a a nur a 5 N 
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also 


1 2tdt (w—2cw+.cf) dw 
JS, a N => C, 5 5 
fi rya—P) RR) m. (w—.c) (w—d) (w—e) (w—f) (w—h) 


1771 
-. v8 oder, in d und e ausgedrückt, (,— lie 5 


1—c 

Es ist also: | 
ER j (w—de) dw ; 

| o(@=1) @—0) (@—4) (oe) (o—f)(o—1) 


“| (w—2cw+cf) dw 
J \o(@—1) (u—e) (u—4) (o—e) (of) (oh) 


ai 0, 0 = an C und damit wegen (16a) und (16b): 


6 01 (w®—de) + 0,0, —2 cw+ cf) = [o®—(d+de) w+d: de] CO = (w—d) (w—h) C. 


EL) Re a A N 


(u) = [ [a O,(@ —de) + 8,0, (0 —2 satef)ldw 
Yo(w—1)(w—e)(w—d) (w—e) (w—f) (o—h) 


(w--d) (w—h) dv 
—=( | nn m 5 
ee re ee 


\(w—d) (w—h) dw 
z(w) = 0 u Le 1 26). 
= ee a 


Damit haben wir — wie es sein muß — die Schwarz-Christoffelsche Formel für unsere _ 
Achtecksabbildung erhalten. Die einzelnen Transformationen sind in Bild 1 veranschaulicht). 
Sie ergeben sich leicht aus den Abbildungsgleichungen. Aus (26) erkennt man, daß D und H 
die einspringenden Ecken des Achtecksbereiches & sind. Die Schlitze der z,- und 2,-Ebene 
können infolge der Bedingungen (23a) und (23b) für a, und a, keine von den Ecken D und H a 
von & ausgehenden Einschnitte erzeugen. [Vgl. im Gegensatz hierzu G1.(42), $ 6,4 und Bild 14!] 4 


83. Die Längen der Achteckseiten 


K,,tK’ und K,,iK, seien die Periodizitätsmoduln der elliptischen Integrale J,(s) bzw. ul) 
der Gl.(24a) und (24b). Dann erhalten wir aus Gl. (22a) und (22b) die VerZweigUDEERUnE FE in: ee 
der z,;- und z,-Ebene und aus (25) die Polygonecken in der z-Ebene (s. Tabelle 1). So 


Tabellel 
Ne. 
A 0 0 
B ıKı & Ey 
c %(Kı +iK}) &%K, 
D &,3&K1 2 
E DICK, %,(2K,tiK;) 
F %(—K,ı +iK)) &(K,+ik,) 
G 0,8; &%(K,+iK,) 
H &,iK, 


m 2 ‚ Pr 5 r PR Kr 
Bus nr IF FR 2, 5 a au ak aa Di 


E A ; 
UN 


ee ; BT FL ef: & air , $ » f 
NY: der ! 
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3 , Aus der Tabelle gewinnen wir für die Längen der Achteckseiten die Gleichungen: 
2 N f 
R Malen &ı Kı, + %K, 

BC Le ’ 
ne . (27); Fi 
{ CED=uK,—%K,") N Br 
DE = %Ks 
j ; ferner wird 
| AB=EF 
‘ re er a ar (28) ; 

k COD=GH 
i DE = HA 
a dies besagt, daß die Achtecke in bezug auf einen Mittelpunkt symmetrisch sind. 
- Die G1.(27) können wir noch einfach umformen. 
; Nach (23) ist! — an oder, da aus (21a) und (21b) die Gleichungen _ 
2 } 2 es 
y 
kı+k 
4 e= —k,; FREE EB Be N FRE 29 
r® - BIER A) 
Ab folgen, 
Nr u WAL, 1+k, 30 
Fa VE a A a - (80). 
R: Setzen wir diesen Ausdruck in (27) ein, so erhalten wir mit einem positiven Proportionalitäts- 
R Y R 1 —k, 
ER; faktor 0 = — —: 
5 “x 
Er h S | edAB=(+k) Kı+(1—k,)K,, 
y% eBO=(1+k)K;, 
AUT ERE e0D=(1+k)K,—(1—k,)K,, 
Br: en gDE = (1— KR: 


| ni Nach Gl. a a (4) drücken wir die er in d- Penkiibucn aus. Mit oe*—= rt wird 


N AB [le + le) + [e) — 9a), 
TE RR TEL e* Bi= iR) + Ra]. 
‚0* CD=[%M)+ % MG) —[%) — v6)» 
o* DE = — 17) [9 (8) — 9 (t2)] - 


al). wel). me) mer) 


i—ı[®), world), mer), uroe(e) Be 
2 ar DR: 2 
| liese Gleichungen durch (10) vereinfachen: 
ABLE HR), 
RR 0.2. (32) 


CD HA) — Bl), 
w o*DE= — 2474 (72) 


F all &ı Br, ist, folgt sowohl aus ar Tatsache, daß D eine einspringende Ecke 
She), mit @ ii 


Pur 


_ und durch Division dieser beiden Gleichungen: 
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Dafür können wir auch schreiben: 
oAB—=KI+ KR}, 
oeBO=2 
eCD=Ki—W' 
oeDE—2% Kr 

Mit (6a) finden wir aus (32): 

e*[BC— (AB— CD] = 2% (& =) — 292(7), 

was wegen (9) >0 ist; d.h.: 

BC zAB—CD. 
Die vier Gleichungen (32) bzw. (33) enthalten drei unabhängige Veränderliche, nämlich zf, 73 


und o. Eliminieren wir diese, so bekommen wir zu den Gl.(28) noch eine weitere Bedingung 


für die Achteckseiten. Nach (33) ist 


* Tt 
M-5R)=-S(ABHOD ......l00., (34), 
F T 
i ME- SU) S(AB-OD .:. 22.000000 (35) 
und mit (31): 
a IB ORT (36), 
Rn ABECD 
K; 2DE 
—2ıd = I u In A art 
| K, AB-—-CD er) 
Setzen wir (36) in (34) und (37) in (35) ein, so erhalten wir 
na BO —)= 7 ( DE 
DR AB GD — 9% a ea I 
2 ii AB+ACOD 242 Dr ZN BESCD 


a EEE 


AB—CD 
Durch Fe: von (6) können wir diesen Ausdruck auf verschiedene Art umformen: 
ER br [i et: a Be 
| BORm) B 
r [ Sn En B20D 
B—( 


d2 [ AD 2 
E 

9 [ A Br 

Nr BO 

» [ 2 2 D 

‚DE | 

Sind drei Seitenlängen (oder Kombinationen wie AB+ CD) vorgehen, so läß 
abhängige Größe aus einer dieser vier Gleichungen ern Ez 2) 
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Eee 


Um die geometrische Bedeutung der Bedingung (38) zu erkennen, kann man aus den 
Gleichungen (33) eine Anzahl von Achtecken berechnen. Es ist am zweckmäßigsten, hierbei kf? 
und %3® als unabhängige Veränderliche zu verwenden, da die Größen K und K’ als Funktionen 
von %k® in Zahlentafeln angegeben sind’). Der Faktor o bewirkt nur eine Ähnlichkeitstrans- 
formation; er soll so angenommen werden, daß BÜ=1 wird. 


Aus (Sb) und (33) folgt 


für 1: eAB=ge0D=Kt und oDE=0 RI 
Tabelle 2 
BE ER EHEN TREE Ir E27 ARCHE 
Kr? A ERS 
RR ) 0,007 | 0,270 | 0,700 | 0,900 | 0,970 | 0,992 | 0,9993 1 
5 L 

0 0 0 0 0 0 0 0 0 CD 
0 1,000 0 ) ) ) ) N) ) 0 DE 
0,000 | 0,001 | 0,018 | 0,056 | 0,098 | 0,133 | 0,159 | o,1s6 | 0,2038 | ©D 
0,007 © | 0,997 | 0,484 | 0,242 | 0,132 | 0,071 | 0,036 | 0,011 ) DE 
0,024 | 0,024 | 0,056 | 0,123 | 0,196 | 0,257 | 0,304 | 0,352 | 0,382 | OD 
3 0,270 © | 1,759 | 0,854 | 0,428 | 0,232 | 0,125 | 0,064 | 0,019 ) DE 
) 0,147 | 0,148 | 0,188 | 0,274 | 0,368 | 0,446 | 0,506 | 0,567 | 0,606 | OD 
0,700 © | 2,251 | 1,084 | 0,548 | 0,298 | 0,160 | 0,083 | 0,024 ) DE 
0,312 | 0,313 | 0,356 | 0,447 |- 0,547 | 0,630 | 0,694 | 0,758 | 0,799 cD 
0,900 © 1.2393 | 1,163 | 0,582 | 0,316 | 0,170 | 0,087 | 0,026 0 DE 
0,501 | 0,502 | 0,545 | 0,638 | 0,739 | 0,824 | 0,889 | 0,955 | 0,997 | oD 
0,970 o 12438 | 1,184 | 0,593 | 0,322 | 0,173 | 0,089 | 0,026 ) DE 
0,710 | 0,711 | 0,754 | 0,848 | 0,950 | 1,035 | 1,101 | 1,167 1,209 cD 
0,992 © | 2452 | 1,191 | 0,596 | 0,324 | 0,174 | 0,089 | 0,026 N) DE 
0,901 | 0,902 | 0,946 | 1,040 |. 1,142 | 1,227 | 1,293 | 1,359 | 1,01 | 0D 
0,9976 © | 2,455 | 1,193 | 0,597 | 0,324 | 0,174 | 0,090 | 0,026 ) DE 
1,097 | 1,098. | 1,142 | 1,236 | 1,338 | 1,423 | 1,489 | 1,555 | 1,597 | 0D 
0,9993 oo | 2,456 | 1,193 | 0,597 | 0,324 | 0,174 | 0,090 | 0,026 N) DE 
1,297 | 1,298 | 1,341 | 1,435 | 1,537 | 1,623 | 1,688 | 1,755 | 1,797 | 0D 
0,9998 & 2,457 | 1,193 | 0,598 | 0,325 | 0,174 | 0,090 | 0,026 0 DE 
h 1,517 | 1,518 | 1,562 | 1,656 | 1,758 | 1,844 | 1,909 | 1,975 | 2,017 | oD 
0,99995 er 2,457 | 1,193 | 0,598 | 0,325 | 0,174 | 0,090 | 0,026 0 DE 
1,774 | 1,774 | 1,818 | 1,912 | 2,014 | 2,100 | 2,165 | 2,231 | 2274 | oD 
' 0,999990 © | 2,457 | 1,193 | 0,598 | 0,325 | 0,174 | 0,090 | 0,026 ) DE 
) 2,030 | 2,031 | 2,074 | 2,168 | 2,270 | 2,356 | 2.421 | 2,487 | 2,530 | OD 
0,999938 x 2,457 | 1,193 | 0,598 | 0,325 | 0,174 | 0,090 | 0,026 0 DE 
oo oo oo oo oo oo oo oo oo CD 
1 o 1! 2457 | 1,193 | 0,598 | 0,325 ! 0,174 | 0,090 | 0,026 ) DE 


In Tabelle2 sind die Seitenlängen CDund DE für verschiedene Werte von k%? und k%° angegeben. 
Die zugehörige Seite AB läßt sich aus der Gleichung: 


en) 


- * 
leicht bestimmen, da ja wegen (39) der zu kj” gehörige Wert n jeweils in der letzten Spalte - 


‚auftritt. Mit BO=1 und GI. (28) erhält man also für jedes angegebene Wertepaar k#’, k#* ein 
+  Treppenachteck. 


= Im &-y-Koordinatensystem des Bildes 2 ist = DE, y=(CD. Entsprechend der Tabelle 2 
liefert jedes Paar kf*, k3° einen Punkt. Die Punkte kf? = const wurden durch Kurven verbunden, 
a ; 2 


22) Siehe [2], 8. 6lfk. 2 
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’ 


die [für %$’—=1 gemäß (39)] die y-Achse in 


treffen und (da für ?—=0 DE= x wird) zur «-Achse parallele Asymptoten besitzen. Diese 
sind am rechten Rand der Abbildung an- 
gemerkt. 

Der Ecke B entspricht der Punkt (—1; 0), 
der Ecke (der Nullpunkt; D(0; y) kann auf 
der y-Achse beliebig angenommen werden, die 
Abszisse von E(#; y) ist ebenfalls willkürlich. 


Mit (28), EF= AB, erhalten wir aus (40): 
OCD+EF=2Y%. 
Dadurch ist der Punkt F(x; yr) bestimmt; 
* 


Cd-fy ] 


Ki aREE 
denn 9, =-— kann man für jeden einem 


Punkt E entsprechenden Wert von kr’ durch 
graphische Interpolation näherungsweise. er- 
halten. In Bild2 ist für ein spezielles Bei- 
spiel der Streckenzug ÜDEF eingezeichnet. 

Unterscheidet sich kf* nur wenig von 1, 
so ist mit großer Genauigkeit 2 Ki =, also 
wegen BC=1 und Gl.(33): 


Deu 


Halten wir nun %3* fest, während %f” 1 ver- 
änderlich ist, so bekommen wir aus (33) und (41): 


Benne ak 2ky Kr 


UaESzeN 
Bild 2 KR. lE 
n9=— 0Dr Ba 


Dies besagt: Ist % größer als eine gewisse Schranke und & = const, so ist die Differenz %— Y 
unabhängig von %f” d.h.: 

Für hinreichend große Werte von %f* gehen die Kurven k#’ — const durch Parallelver- 
schiebung auseinander hervor. Aus Tabelle 2 ersieht man, daß es bei einer Genauigkeit auf 
drei Dezimalstellen genügt, die Schranke 1?—0,9998 bzw. y=1,8 zu wählen. 


Bild 2 soll nur die Form der auf unsere Methode abbildbaren Achtecke erkennen are F 
Für die Auffindung der einem vorgegebenen Polygon entsprechenden Moduln k, und ky ist sie 
ohne Bedeutung (vgl. $5). 


e Es läßt sich zeigen, daß die durch die beiden orthogonalen Halbkreise | 
R=uw=Yde e” (<yp<a) EC =w=c+Y(e—1)(e—d): €? (0 <p <a) 


und die reelle Achse der w-Ebene gebildeten Kreisbogendreiecke in die vier durch die Strecken P, ®, 
und M,N, bzw. P,Q, und M,N, getrennten Teilrechtecke der z,- bzw. 2,-Ebene übergehen. Bei 
der praktischen Durchführung der Transformation genügt es daher, eines dieser Dreiecke durch 
2, (w) bzw. 2,(w) auf je ein Rechteck der 2,- und 2,-Ebene abzubilden. Die Abbildung der Enası 
Halbebene 9 erhält man durch Spiegelung an den entsprechenden Strecken. | (i 


Däs Prinzip unserer Achtecksabbildung bestand darin, daß wir die w-Halbebene 9 auf zwei 
mit je zwei geradlinigen Schlitzen versehene Rechtecke der z,- bzw. 2,-Ebene abbildeten und diese 
überlagerten. Die Seitenlängen des Achtecks mußten hierbei gewisse Gleichungen erfüllen. 
Jedem Verzweigungspunkt der w-Ebene entsprach je ein Eckpunkt der beiden geschlitz 
Rechtecke. Diese Zuordnung erscheint willkürlich. Zunächst ist anzunehmen, daß bei geeig 
Änderung der Rechteckseiten durch zyklische Vertauschung der Buchstaben in einer 
Ebenen — gegebenenfalls bei gleichzeitiger Drehung der Rechtecke um den Wink: 
Überlagerung wieder Treppenachtecke liefert, deren Seitenlängen jedoch anderen ] 
genügen. Bei ausführlicher Erörterung dieser Frage erkennt man aber, daß sich auf di 
außer einem Spezialfall (7,—=7,) der besprochenen Achtecksabbildung nur die z 
Achtecken BL BelbLldlichen Polygone ergeben. Wir erhalten also keine wesent; 


1) eh . 


ZZ. angew. Ma 
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$ 4. Ausartungen 
1. Nimmt einer der beiden Moduln A} und %$ der Gl.(31) den Wert Null oder 1 an, so artet 


das Achteck aus. Unter Berücksichtigung von (8a) und (8b) ergeben sich aus (33) und (28) 
‚die in Tabelle 3 zusammengestellten vier Möglichkeiten. 


N : Tabelle 3. 

a) it 0 bit-ı | our=o0 | s#=1 
a a T 9 x 
eAB=-eHEF | G+tar'Kt > Kt+7 Kr 
eBC=ofF@ 0 1 aRV 2.Kr 
wu wu T J 
e0D=-e0H | Z—-kYK © Kr K% 
oDE=oHA 2ER Kr 2kt Kr oo ) 


Der Fall d) ist unbedeutend, da er offenbar nur Rechtecke liefert. In den übrigen Fällen 
können wir jeweils eine Länge gleich 1 setzen. Von den anderen endlichen Längen ist nur eine 
unabhängig, da ja nur mehr ein Modul veränderlich ist. 


Y 

> Y 

A 
Mn 
re 
Baader 
4% A 
i Gun — £ 9 05 0 05 
r£ ® ber Bild 3a Bild 8b Bild 8e 


Tabelle4a:k*=0; GH+AB=1 


ze* 0 | 0007 | 0,270 | 0,200 | 0,900 | 0,970 | 0,002 | 0,0093 | 1 
Ei s=@H | 0 0,001 | 0,044 | 0,138 | 0,241 | 0,326 | 0,392 | 0,158 | 0,500 
ya HA | oo 2,457 | 1,193 | 0,598 | 0,325 | 0,174 | 0,090 | 0,026 0 


Tabelle4b:k*=1; BO=1 


| 0,007 | 0,270 | 0,200 | 0,900 | 0,070 | 0,992 | 0,0003 | 1 
0,998 | 0,911 | 0,724 | 0,519 | 0,348 | 0,217 | 0,084 0 
2,457 | 1,193 | 0,598 | 0,325 | 0,174 | 0,090 | 0,026 0 


'Tabelle4c:k*=0; AB—-OD=1 


0,007 | 0,270 | 0,500 | 0,700 | 0,000 | 0,970 | 0,992 | 0,0976 | ı 
| 0,001 | 0,033 | 0,090 | 0,161 | 0,321 | 0,505 | 0712 | 0,902 | ® 
1,465 | 0,397 | 0,150 | 0,001 | 0,026 | 0,008 | 0,002 | 0,001 | 0 


en 4a, 4b und 4c geben etliche Wertepaare an, durch die je ein ausgeartetes 
ımt ist. In den Bildern 3a, 3b und 3c sind mit Hilfe der den Tabellen zu entneh- 
die Kurven y(x) gezeichnet. Ihre geometrische Bedeutung für die Gestalt der 
Br | | | 26 


le ZU 
a = 


de 


wa 


“ 
2 
“ 
» 
vl 
. 
x 


a a 
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ausgearteten Achtecke ist aus den Bildern 4a, 4b und 4c zu erkennen: Jede dieser Kurven ist Bi 
der geometrische Ort für eine bestimmte Polygonecke, nämlich im Fall a) für 7, im Fall b) für E Se 
und im Fall ce) für (©. 


Bild 4a Bild 4b Bild 4c 


2. Bei der von Bergmann behandelten Abbildung von Treppensechsecken (A’, B’, 0", 
D’, E’, F’) ergeben sich die Seitenlängen ®): 


en oe DE=Kh—K,, 
PR | eEPF=—Kı+ KR 
ol’ D — DIR, BD 


wobei, falls D’E’ >0 ist, die einspringende Ecke in E’ liegt. Auch hier gibt es wieder grund- 
sätzlich vier Möglichkeiten von Ausartungen, die wir in Tabelle 5 zusammenstellen. 


Tabelle5 

)h=0 | )h=1 I )hR=0 | )h-l 
a = T+R © &tz 
oB’C’ Stk, oo &,\+> oo 
o0’D’ © 7 2& 2K, 
eDP | 5-, Ba ee 
oB’F ze ne Ri +5 
oF’A’ 2K, 2K, 7 0 


In den Fällen &) und ö) bekommen wir nur ausgeartete Vierecke, deren Abbildung nichts ; 


Neues darstellt. Ersetzen wir bei ß) Le durch Ki RL durch R,, so erhalten wir dieselben Seiten- 
längen wie im Fall y), nur in. der entgegengesetzten Reihenfolge. Dies bedeutet, daß die durch p) 
bzw. y) bestimmten Polygone durch Spiegelung auseinander hervorgehen. Mit B=4, (= 
D’—(, E’—D erweist sich y) als mit dem Fall c) (k$—=0) der oben besprochenen ee 
Achtecke identisch, wie man durch Vergleich mit Tabelle 3 erkennt. er, 


$ 5. Die praktische Durchführung der Abbildung 


Um die Abbildung der w-Halbebene 9 auf ein Treppenachteck, das den Be (28) 
und (38) genügt!*), praktisch durchzuführen, lösen wir zunächst das Parameterproblem: Aus den 
gegebenen Seitenlängen erhalten wir mit den Gl. (36) und (37) die Größen K,/Kı und K//RK,. Inder 
Zahlentafel [2], in der K’/K als Funktion von A? angegeben ist, finden wir k? und k2. (29) liefert ER 
die Verzweigungspunkte d und e, aus denen wir durch (18) c, fund h bestimmen können, Damit sız 
sind sämtliche Verzweigungspunkte bekannt. K/ ist als Funktion von %? ebenfalls der Tafel [2] 
zu entnehmen; aus (27) läßt sich somit «,, aus (30) % und aus (23b) die (für die Transtorzma ron 
nicht erforderliche) Konstante C der Abbildungsgleichung (26) berechnen. > A: 

“ Die Transformation eines Punktes w* von 9 in die z-Ebene führen wir schrittweise durch 
Mit Gl.(20a) und (20b) erhalten wir s(w*) und t(w*) und sodann aus (22a) bzw. (22 b) ‚die > 


Hr 


Punkte 2, (w*) und 2,(w*). Da keine Tafeln der Euren InlEsseie für Komp BLEND ente 


) Siehe [1], 8. 324, wo wir & durch 1/0, K, durch A und R, durch R, ersetzen. 

14) Ob ein vorgegebenes Treppenachteck außer (28) auch die Bedingung (38) erfüllt, k 
in $3 dargelegt wurde, an Bild 2 einfach feststellen, falls keine Be a an di 
gestellt werden. “, 


w es 
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existieren, wird man die Abbildung 2, (s) folgendermaßen durchführen: Man überträgt durch die 
Umkehrung s(2,) = sn (k,/z,) ein orthogonales Kurvennetz der z,- auf die s-Ebene und bestimmt 
durch graphische Interpolation das Bild von s(w*) in der z,-Ebene. Das Entsprechende gilt für 
die Transformation 2,(t). G1.(25) liefert schließlich den gesuchten Punkt z(w*) im Achtecks- 
bereich ©. 

Ist das vorgegebene Polygon ausgeartet, so fallen, wie mit Hilfe von (29) und (18) leicht 
nachzuweisen ist, gewisse Verzweigungspunkte zusammen. Die Zwischenabbildungen verein- 
fachen sich hierbei. Gehen wir beim Fall c) (s. $4) von der v-Ebene aus, so ergeben sich dieselben 
Transformationsgleichungen wie in dem damit identischen Fall y) von ausgearteten Sechsecken, 


$ 6. Praktische Anwendungen und Erweiterung der bisherigen Ergebnisse 

1. Wir können die w-Halbebene 9 durch ein elliptisches Integral auf ein Rechteck Rt einer 
£=&--in-Ebene mit zu den Koordinatenachsen parallelen Seiten so abbilden, daß vier vor- 
geschriebene Punkte der reellen Achse in die Ecken von Rübergehen. Da die Transformation z(%) 
bekannt ist, besteht somit die Möglichkeit, den Achtecksbereich ® durch eine Funktion 
&(2) = &(w(z)) so auf N abzubilden, daß die Eckpunkte des Rechtecks vier beliebig vorgegebenen 
Randpunkten S,, S,, 8;, 8, von ®entsprechen. Betrachten wir in der ö-Ebene eine Parallelströ- 
mungder Geschwindigkeit 1 in Richtung der negativen &-Achse, so ist &(z) das komplexe Strömungs- 
potential einer Strömung der z-Ebene 5). Den Potentiallinien &—=const entsprechen in & wieder 
Potentiallinien, den Stromlinien „—=const wieder Stromlinien. Insbesondere sind die einzelnen 
den Seiten des’ Rechtecks R zugeordneten Randstücke von & Potential- bzw. Stromlinien. Es 
dürfen daher je zwei Randstücke als Potential- bzw. Stromlinien vorgeschrieben werden. Fallen 
zwei oder zweimal je zwei der Randpunkte S,, S,, S,, 8, zusammen, so artet N in einen einseitig 
bzw. beiderseits unbegrenzten Parallelstreifen aus; sind drei dieser Punkte identisch, so ent- 
spricht dem Bereich & ein Quadrant der &-Ebene. 

Durch Spiegelungen erhalten wir das komplexe Strömungspotential in komplizierteren, 
z. T. mehrfach zusammenhängenden Bereichen. Die Abbildungen können auf die Strömungslehre, 
die Elektrotechnik, teilweise auch auf Probleme der Wärmeleitung und der Torsion angewandt 
werden. Freilich muß stets vorausgesetzt werden, daß das Problem als eben betrachtet werden 
kann. 

2. In den Bildern 5 bis $ ist eine Anzahl von Bereichen zusammengestellt, für die wir durch 
Spiegelung der Achtecke & die Strom- und Potentiallinien ermitteln können. Selbstverständlich 
ist es in jedem Fall möglich, Strom- und Potentiallinien zu vertauschen. Da den vier Ecken 
von X nicht unbedingt Eckpunkte von & entsprechen müssen, können im folgenden vielfach 
Teile von Symmetrielinien als Stromlinien vorgeschrieben werden, wie dies aus einigen Bildern 
(öd, 6e, 6f, 6g, 6i, 8c) zu ersehen ist. Die einzelnen Figuren sind leicht verständlich. Sie ergeben 
sich durch einfache, mehrmalige oder unendliche Spiegelung der jeweils ersten Figur (a), welche 
das besprochene Achteck bzw. eine seiner Ausartungen darstellt. Je nach der Zuordnung der 
Strom- und Potentiallinien lassen einige Bilder (öf, 5g, 6d, 6e, 6k, 61, 7b, 7c, Sc) verschiedene 
Deutungen zu. So kann Bild 5f als symmetrische Umströmung der rechteckigen Kerne oder 
auch als Zirkulationsströmung um jeden einzelnen Kern mit abwechselndem Richtungssinn auf- 
gefaßt werden. Eine Abbildung, die sich von de nur dadurch unterscheidet, daß die Einkerbungen 
einander gegenüberliegen, kann mit Hilfe von Treppensechsecken durchgeführt werden *). 

3. Wenden wir auf die z-Ebene die Exponentialfunktion an, so gelangen wir zu den in den 

_ Bildern 9 bis 12 angegebenen Bereichen. 

Legen wir in Bild 5a A in den Nullpunkt und setzen BC+-DE=n, so liefert die Trans- 
formation Z=e2einen Bereich von der Art des Bildes 9a und in Verbindung mit einer Spiegelung 
die Figur (b). Entsprechend bekommen wir durch Z= e!*und mehrmalige Spiegelung einen innen 
und außen mit je n Kerben versehenen Kreisringbereich. Dieser ist für n—2 in (c) angegeben. 

In Bild 6a sei GH + AB=n; Bliege im Nullpunkt, A auf der reellen Achse. Dann liefert 
Z =e-i2 die durch die Kurve des Bildes 10a begrenzte obere Halbebene. Durch Spiegelung 
ergibt sich die symmetrische Strömung bzw. Zirkulationsströmung um einen mit einer Aus- 

"kerbung versehenen Kreis (b) und durch Inversion an diesem das Innere der Kurve (c). Der 
Wirbelpunkt dieser Strömung kann auf der Symmetrielinie beliebig angenommen werden. Aus (a) 


B* bekommen wir durch Ziehen der n. Wurzel den Winkelraum (d), durch Spiegelung die Zirku- 
_  Jationsströmung um die Zahnradkurve (e) — im Fall n=2 (f) auch die Umströmung in einer 


der beiden Symmetrierichtungen — und durch Inversion von (e) am inneren Kreis den endlichen 
"Bereich (g) mit einem Wirbel im Mittelpunkt. Setzen wir (e) nach innen analytisch fort, so 


ja wer 
15) Vgl. [3], 8.63 ff. 
16) Vgl. [7], 8. 129ff. 
26* 
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erhalten wir, falls wir in (a) den Kreisbogen AB als Potentiallinie, die übrigen Randstücke als 
Stromlinien betrachten, die Quellströmung (h) um symmetrisch angeordnete Hindernisse, speziell, 
indem wir (e) durch (b) ersetzen, um ein Hindernis (i). Ist dagegen in (a) außer A B auch der 


c 
D 
Be 
ne | ©, ‚8 
b c 


a 
Bild 9 


Se 
PERERaGE 


L k 
. Bild 10 


N 


28 


Bild 11 


J 


ee, b 


yon @ nach Unendlich führende Teil der reellen Achse eine Potentiallinie, so stellt (h) die Zirku- 
lationsströmung mit abwechselndem Richtungssinn um jedes einzelne Hindernis dar (deren 
Anzahl in diesem Fall gerade sein muß). Entsprechend bekommen wir aus (d) die Umströmung 
eines Kreisringsektors in einem Winkelraum (k). 
Bild 11a zeigt den Bereich, den wir durch Anwendung der Exponentialfunktion auf Bild 7a 
BE Be. Run RT De 
Eu erhalten, wenn dort BC+ DE=nist. Ist speziell BÜ=DE= 5 ‚so ergibt sich durch Spiege- 


ve N STR i 
Jung die Figur (b). Durch Radizieren und Spiegeln gelangen wir von (a) zu (c), (d) und ähn- 


7 ER 


u 3 
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Ist in Bild 8a A B=rund liegt Aim Nullpunkt, so liefert Z= e* den Bereich des Bildes 12a. 
Durch Radizieren und Spiegeln bekommen wir z.B. (b) und (c), aus (b) durch Inversion (d) 


und aus (c) durch analytische Fortsetzung am inneren Kreis (e). Ist speziell in Bild 8a C cD= 5 ; 


so erhalten wir aus Bild 12a durch Spiegelung (f) und hieraus durch Inversion (g). Aus Bild 8f 
folgt mit Z= ei die Strömung (h). 


u x 
= e © Gi 


Bild 12 


h 


Auch auf die von Bergmann behandelten Treppensechsecke (Bild 13a) können wir die 
e-Funktion ausüben. Wir bekommen so die Figur (b), durch Spiegelung (c) und hieraus durch 
Inversion bzw. analytische Fortsetzung am äußeren Kreis bzw. am inneren Kreisbogen (d), (e) 
und (f). Ziehen wir aus (b) z.B. die 4. Wurzel, so erhalten wir nach mehrmaliger Spiegelung 
statt (c) das Bild (g). Entsprechend lassen sich die Figuren (d), (e) und (f) verallgemeinern. 


4. Alle gezeigten Abbildungen sind (außer speziellen Symmetriebeziehungen) noch einer “ 
gewissen Einschränkung unterworfen; diese entspricht der a (38) für die Seitenlängen 
der Treppenachtecke. 


Bild 13 


Pr i 
Unsere Zusammenstellung der konformen Abbildungen (Bild 5 bis 13), die 
ausführlich besprochenen Methode der Treppenachtecke (bzw. Treppensechsecke) 
werden können, läßt sich durch Anwendung anderer Funktionen, z.B. der trig 
der hyperbolischen, der linearen Funktionen oder der Potenzen, Haliebiu erweil 
meinen ergeben sich hierbei ziemlich komplizierte Figuren. 2 


te 


'Z. angew. Math. Mech. 
Bd. 31 Nr. 11/12 Nov./Dez. 1951 


Kleine Mitteilungen 387 


Die Achtecksabbildung können wir noch etwas abändern: Lassen wir die G1.(23 a) und (23b) 
für die (reellen) Größen a, und a, fallen, so geht (26) über in einen Ausdruck von der Form 


(w—Y) (w—) dw 


F E 


Al 


0 —— 
| \o(@—1) (—0) (0-4) (@—e) (@—f) (@—h) — 


(42), 


wo yund ö von a, und a, abhängige reelle Konstanten sind. 
Aus der Schwarz-Christoffelschen Formel folgt, daß 4 H 


das durch (42) gegebene Polygon der z-Ebene in den 


A B 


Ecken A, B, 0, D, E, F,@,H jeweils den Innenwinkel x/2, Bild 14 

in zwei weiteren Ecken den Winkel 2x besitzt. Wegen der 

Symmetrie der Rechtecke der z,- und der z2,-Ebene liegen diese beiden Eckpunkte symmetrisch 
zum Mittelpunkt des Polygons. Dieses ist also von der in Bild 14 angegebenen Gestalt. 
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KLEINE MITTEILUNGEN 


Bemerkungen zur Konvergenz des Iterations- 
verfahrens bei linearen Gleichungssystemen!). 


I. Einleitungund Programm 


Das lineare Gleichungssystem, dessen Lösung ge- 
sucht ist, sei gegeben in der Form: 


Be re a Re Sıhh 


T 
Hierbei bedeutet r = ) den gesuchten Lösungs- 
In 


& ; 
vektor, e = (Ü) einen vorgegebenen Vektor und A 


: 27 
eine vorgegebene quadratische Matrix mit n Zeilen 
f und Spalten. Sämtliche Elemente von c und X seien 
) reelle konstante Zahlen. Um die Lösbarkeit von (1) 
für beliebiges c zu sichern, wollen wir voraussetzen, 
daß der Rang von U — & gleich n sei. 
Die Iteration des Gleichungssystems werde in der 


‚folgenden Weise: 
Im limtt: 2... (2) 


durchgeführt. Bezeichnet man die gesuchte Lösung 
von (1) mit r, und führt die Abweichungen der Ite- 
 rationen von dieser Lösung: dtm = Im — to in (2) ein, 
eo wird &us (2): 
# ms FU dem 2...) 


Für die Konvergenz des Iterationsverfahrens ist 
dann bekanntlich notwendig und hinreichend, daß der 
 Absolutbetrag sämtlicher Eigenwerte von W kleiner 
als 1 ist. 

Die Anwendung dieses Konvergenzkriteriums setzt 
also die Bestimmung des absolut größten Eigenwertes 

von X, den wir A, nennen wollen, voraus. Diese erfolgt 
bekanntlich aus der Berechnung aufeinanderfolgender 
 — Potenzen Av, Av+1, Yv+2, ... der Matrix X. Hier- 
 beimuß man, wenn außer dem Eigenwert mit größtem 


jiese Arbeit ist entstanden im Anschluß an gemeinsame 

or E. Krahn und dem Verfasser aus dem 
‚1949. 
67 ar 


\ 


AND 


Absolutbetrag noch Eigenwerte mit fast ebenso großen 
Absolutbeträgen auftreten, » ziemlich groß wählen, 
um A, bestimmen und damit die Frage nach der Kon- 
vergenz des Iterationsverfahrens entscheiden zu kön- 
nen. Außerdem kommt erschwerend hinzu, daß sich 
die Berechnung von A, verschieden gestaltet, je nach- 
dem ob X einen reellen Eigenwert oder ein konjugiert 
komplexes Eigenwertpaar mit größtem Absolutbetrag 
besitzt (ganz abgesehen von den Fällen, in denen das 
Auftreten mehrerer Eigenwerte bzw. Eigenwertpaare 
von gleichem größten Absolutbetrag die Berechnung 
von A, noch weiter kompliziert). 

Aus diesen Gründen kann die Bestimmung von }, 
zur Untersuchung der Konvergenz des Iterations- 
verfahrens zuweilen ziemlich umständlich werden. 

Außer der oben genannten notwendigen und hin- 
reichenden Bedingung für die Konvergenz des Ite- 
rationsverfahrens ?) sind in der Literatur zahlreiche 
einfache, aber nur hinreichende Konvergenzkriterien 
bekannt. 

Im folgenden — in III — sollen einfache notwendige 
und hinreichende Bedingungen für spezielle, praktisch 
wichtige Arten der Konvergenz des Iterationsverfah- 
rens aufgestellt werden. Durch diese speziellen Arten 
der Konvergenz werden, wie sich zeigen wird, alle 
überhaupt möglichen Fälle von Konvergenz erfaßt. 
Diese im folgenden aufzustellenden Konvergenz- 
bedingungen werden die Kenntnis der Eigenwerte von 
U nicht voraussetzen. 

Ehe wir zu diesen Bedingungen kommen, wollen wir 
in II einige anschaulich-geometrische Betrachtungen 
über den Konvergenzprozeß anstellen. 


II. Anschauliche Darstellung des 
Konvergenzprozesses 
Der Grad und die Art der Annäherung der Ite- 
rationen fm für großes m an die Lösung r, hängt 
2) C. E. Berry: Ann, Math. Statistics 16 (1945), p. 398/400 


stellt die mit der obigen Bedingung äquivalente notwendige und 
hinreichende Konvergenzbedingung: lim A? = 0 (0=Nullmatrix) 


v0 
auf, die indessen praktisch zum Nachweis der Konvergenz bzw, 
Divergenz kaum verwendbar sein dürfte, 
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ab ehe 


wesentlich von A, ab. Wir setzen in II. stets jA,|<]1, 
d.h. Konvergenz des Iterationsverfahrens voraus. Die 
absolut genommen nächst größeren Eigenwerte von 
A bezeichnen wir mit A,, A,, usw. Mehrfache Eigen- 
werte werden ihrer Vielfachheit gemäß gezählt; bei 
mehreren verschiedenen Eigenwerten mit gleichem 
Absolutbetrag erfolgt die Numerierung im Sinne 
wachsender Winkelargumente. 

Im folgenden beschränken wir uns auf einige wenige 
häufig vorkommende Möglichkeiten für A,, da diese 
bereits zur Erläuterung unserer Betrachtungen in III 
genügen werden. 


1. Fall:A,sei ein einfacher reeller Eigenwert von 
X, und es sei |A,| > |A,|. 

Dann gilt bekanntlich für hinreichend großes m 
(in asymptotischem Sinne): 


Irm+ı = hı' im 
oder wegen: (3): 
(X — 4, &)dım = 0 v% . (5). 


Aus 4 ist ersichtlich, daß für hinreichend großes m 
die Abweichung dr, von der Lösung r, sich mit jedem 
Schritt auf das |A,|-fache verkleinert (|A,|<1). 
Aus (4) folgt weiter, daß für hinreichend großes m die 
Annäherung an die Lösung mehr und mehr auf einer 
Geraden erfolgt, und zwar für 1, >0 von einer Seite 
her, für A,<0 oszillierend. Aus (5) folgt, daß diese 
Gerade die Richtung des zu A, gehörigen Eigenvektors 
y, von W hat, also von dem Ausgangsvektor des Ite- 
rationsverfahrens völlig unabhängig ist. 

Uns interessiert für das folgende lediglich, daß von 
einem hinreichend großen m an in diesem Falle die 
Konvergenz stets monoton erfolgt, d.h. daß für hin- 
reichend großes m stets |dr,,.,|< |dtm| ist. 


...(4) 


2.Fall:A, und), = —A, seien einfache reelle 
Eigenwerte von X, und es sei |A,| = |A,!>|A;|. 
Dann gilt entsprechend (in asymptotischem Sinne): 


Imre =}: dim . + (6) 
oder wegen (3): 


(A? — 2: E)dim = 0 TI: 


Aus (6) und (7) folgt: Da der Rang von A? — AE€ 
jetzt n — 2ist, hat WA? zweizu den Eigenwerten/? =} 
gehörige linear unabhängige Eigenvektoren y, und 4,. 
Die Annäherungen an die Lösung t, erfolgen für hin- 
reichend großes m mehr und mehr auf zwei Geraden, 
diein der durch i), und 4), aufgespannteri Ebene liegen. 
Und zwar wechseln die Näherungswerte mit jedem 
Schritt die Gerade, so daß dym, drnız» Hmıas »*- 
auf einer und dr„ı1 Ulnız mis, --- auf der 
anderen Geraden liegen. Die Lage dieser Geraden in 
der durch die beiden Eigenvektoren bestimmten Ebene 
hängt von dem Ausgangsvektor des Iterationsver- 
fahrens ab. Beim Übergang von einer Näherung dry 
zur übernächsten dy,,, ı „, die auf der gleichen Geraden 
wie dim liegt, bzw. von dy,,.; zu dr,,„ı.; nimmt die 
Abweichung von der Lösung auf den }3 = }}-fachen 
Teil ab. 

Dagegen gestatten die Gl. (6) und (7) keine Aussage 
über die Änderung der Abweichung beim Übergang 
von dm zu dy,,., bzw. von dr. zu dio, d.h. 
beim Überwechseln von einer Geraden auf die andere. 
Hierbei kann es nun durchaus passieren, daß 
trotz konvergentem Iterationsverfahren stets z.B. 
drmzıl> Idim|, Idrmzsl> Idimzo| usw. ist). 

LEN 

3) Beispiel: Y= ie ü mit den Eigenwerten Aa = + 1. 

5 


Geht man z.B. von dt, = (1) aus, so liefert das Iterationsver- 
fahren: E ; 


drı | dt, | dia | dis | 
1 | 0,3125 0,2500 0,078125 \ 
n. 2,7500 0,2500 0,687500 | .... 


Offensichtlich ist |dr,| >|drzıl, »[dx,]| > |dr;| usw. 
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Natürlich muß dann |dr„.s|< drysıl dimrs 
< Idrnzsl usw. sein, 


Für uns ist lediglich wesentlich, daß in diesem Falle 
die Konvergenz bei noch «so großem m nicht not- 
wendig monoton zu sein braucht. 


3. Fall:}, =a+bi undd,=a— bi sei ein ein- 
faches komplexes Eigenwertpaar von W, und es sei 
4,| = |},| >]A,|. Dann gilt für drei aufeinander- 
folgende Abweichungen dm, dr,ı1> 4%,.,9 für hin- 
reichend großes m bekanntlich (in asymptotischem 
Sinne) die Beziehung: 


Uimig = — (a + B)dim + 20 - dry‘) (8): 


Aus ihr folgt, daß von hinreichend großem m an alle 
Näherungen sich mehr und mehr einer festen Ebene 
durch den Lösungspunkt x, nähern. Diese Ebene er- 
gibt sich aus den beiden linear unabhängigen Lösungs- 
vektoren ), und }), von : 3 


(A-1,H9A-,99=0....0), 


wie man durch Umformung von (8) unter Benutzung 
von (3) erkennt. 


Verbindet man die aufeinanderfolgenden Nähe- 
rungen Im, Eıı> Imra --- geradlinig miteinander, 
so erhält man ein Polygon, das sich spiralartig der 
Lösung r, nähert, wobei sich die einzelnen Polygon- } 
strecken immer mehr der durch 4, und }), ausgespann- 
ten Ebene nähern. 


Um zu untersuchen, ob die Konvergenz in diesem 
Falle monoton sein muß, bilden wir nach (8): 


Ins = — (+ B)dinıı + 20 Armin» 
In+z = — 2a(a? + d?)drm + (3a? — bB)dy„ıı- 
Setzt man abkürzend: 4 
— (a? + b’)dym = gm, 20 dr 11 = dmti> 
so wird: | 
Ilmya — am + äm+ı> | 
a? — b2- . (10). 


3 
Um+3 = 20 mt, " dmti J 
Aus (10) ist ersichtlich, was uns hier allein inter- 
essiert, daß auch in diesem Falle die Konvergenz des 
Iterationsverfahrens bei beliebig großem m durchaus 
nicht moncton zu sein braucht’). 


III. Spezielle Arten von Konvergenz 


Die Ausführungen in II haben ergeben, daß nicht 
jede Konvergenz monoton zu sein braucht. Von prak- 
tischem Interesse ist nun offenbar die Frage: Wie kann KAN 
man aus der Matrix U schließen, ob für sie das Ite- , 
rationsverfahren monoton konvergiert oder nicht ? 


Wir wollen uns hier indessen nicht auf die monotone 
Konvergenz beschränken, sondern unseren Unter- 
suchungen eine etwas allgemeinere Art der Konver- N 
genz zugrunde legen. Und zwar wollen wir frrdern, 
daß wir bei der Iteration (2), von einem beliebigen 
Vektor 1, = to + dr, ausgehend, erst nach 9 Ite- 
rationsschritten immer notwendig zu einem Vektor 
%p+1 = Lo + dr,;ı kommen, der von der Lösung x, 
weniger abweicht als ,, d.h. für den |dr,, 
ist. Ist diese Forderung erfüllt, so wollen 
das Iterationsverfahren sei für die Matr 
gent nach p-ter Art. Für p=1 (Konv 
RE HS a 


% Vgl.z.B. L. A. Pipes: Journal of the Fr 
225 (1938), p. 437/454. ABER AN 
°) Eine grobe Abschätzung liefert aus ( 
sicherlich nicht monoton ist, wenn — a? + 
la] >% und ®< 30° —2|a] ist. Diese 
füllt z.B, für 4 = 3: b=1l. 
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Art) erhält man hieraus den Fall der durchweg mono- 
tonen Konvergenz. 


Für die Konvergenz p-ter Art des Iterationsver- 
fahrens läßt sich nun leicht ein notwendiges und hin- 
reichendes Kriterium angeben. Aus |dr, ‚| <|dr;| 
oder dro+1 -dr„41<drı dr, folgt nämlich unter 
Verwendung von (3): drj - (AP) - AP - dr, <dr; dr, 
oder: 


dz(E — (AP YAr)dı, >0°%).... (1) 


Die linke Seite von (11) stellt eine quadratische 
Form in den Komponenten von dr, dar. Diese soll 
. nach (11) positiv definit sein. Dazu ist bekanntlich 
notwendig und hinreichend, daß die Abschnittsdeter- 
minanten von & — (X?)’- X? sämtlich positiv sind, 
Notwendig und hinreichend für Konver- 
genz. p-ter Art ist also, daß sämtliche Ab- 
schnittsdeterminanten von € — (AP)’ - XP posi- 
tiv sind. Speziell liegt durchweg monotone 
Konvergenz des Iterationsverfahrens dann 
und nur dann’ vor, wenn alle Abschnitts- 
determinanten von E—W-XW positiv sind. 


Daß Konvergenz p-ter Art für jede natürliche Zahl 
p wirklich auftreten kann, sieht man an dem Beispiel: 


) 5 | 
el (& 0) (k ganze Zahl > — 1). 
220 
Hier wird 
E— (UP — ( ET ee cn: ‘) 
0 1 


Die für Konvergenz p-ter Art notwendige und hin- 
reichende Bedingung ergibt hier: 


(Der + ermeern<ı, 
Sie ist offenbar nur erfüllt für 2P —2(k+1)>1 
ı 
oder p>k+3,.d.h. für die Matrix YW= % 2) 


ist das Iterationsverfahren nach p=k-+ 2.ter Art 
konvergent. Läßt man k die Zahlenreihe —1,0, +1, 
+2, +3, +4, ... durchlaufen, so erhält man in 
unserem Beispiel also Konvergenz lter, 2ter, 3ter, 
4ter, 5ter, 6ter, ... Art. 


Es soll jetzt noch gezeigt werden, daß die Konver- 


= genz p-ter Art (p=1, 2,3, ...) sämtliche überhaupt 
möglichen Fälle der Konvergenz des Iterationsver- 

= fahrens umfaßt, d.h. daß jede Konvergenz schlecht- 
) - hin des Iterationsverfahrens eine Konvergenz p-ter 


Art für ein ganz bestimmtes p in dem von uns defi- 
x  nierten Sinne ist ?). 
Aus der notwendigen und hinreichenden Konver- 
_ genzbedingung lim A’ — 0 des Iterationsverfahrens 


A x vo | 
folgt nämlich: Bei Konvergenz des Iterationsverfah- 
 rens kann man eine hinreichend große ganze positive 
erg ‚stets so angeben, daß jedes Element von AM 
dem Betrage nach < — wird, wobei & eine belie- 
ne positive Zahl und n die Zahl der Zeilen oder 
al von X bedeutet. Dann aber ist jedes Element 
er I)". AM dem Betrage nach kleiner als e. Es 
iegen dann in der Matrix € — (XM)’AM die 
‚der Hauptdiagonale von €, so daß sich die Ab- 
ittsdeterminanten von E — (AM) - AM um so 
ron 1 unterscheiden werden, je kleiner & ist. 
ügend großes M sind also, Konvergenz des 
ationsverfahrens vorausgesetzt, sicher alle Ab- 


R - 
ee 
ee 


rich soll die transponierte Matrix bzw. den trans- 
ch so 
tor bedeuten. Bike 
t mich E. Krahn aufmerksam gemacht. Von ihm 
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schnittsdeterminanten von & — (UMy’'AM positiv. 
Ist die kleinste Zahl M, für die das zutrifft, M = M, 
so liegt offenbar Konvergenz 9 = M,-ter Art vor. 

Unsere Erkenntnis, daß jede Konvergenz des Ite- 
rationsverfahrens eine Konvergenz p-ter Art für ein 
ganz bestimmtes p darstellt, ermöglicht es uns nun 
sofort, jedes konvergente Iterationsverfahren auf ein 
solches, das nach erster Art konvergiert, zurückzu- 
führen. 

Hat man nämlich erkannt, daß das durch (2) er- 
klärte konvergente Iterationsverfahren nach p-ter Art 
konvergiert ®), so konvergiert offenbar das auf 


mt BEP Pre LACHe ange- 


wandte Iterationsverfahren nach erster Art. 

Wir haben in unseren obigen Ausführungen die 
Forderung der monotonen Konvergenz (jeder Ite- 
rationsschritt führt näber an die Lösung heran als der 
vorhergehende) durch die Einführung der Konvergenz 
p-ter Art abgeschwächt. Man kann nun natürlich 
auch umgekehrt die Forderung der monotonen Kon- 
vergenz verschärfen, indem man verlangt, daß die Ab- 
weichung von der Lösung r, nach jedem Schritt einen 


I 


kleinen Bruchteil NV der Abweichung vor diesem 


Schritt nicht übersteigen darf. Die Bedingung dafür 
1 


lautet offenbar: dr, ,1' Imıı< ws dr, dtm oder 
bei Benutzung von (3): 
1 
dr), 7 E— Wa)-din SOSE F 2010: 


Soll also nach jedem Iterationsschritt die Abweichung 


1 
von der Lösung r, stets kleiner als x tel der Ab- 


weichung vor diesem Schritt sein, so ist dafür not- 
wendig und hinreichend, daß sämtliche Abschnitts- 


1 
determinanten der Matrix; E — WX positiv sind. 


Dies kann man — Konvergenz des Iterationsver- 
fahrens vorausgesetzt — stets erreichen, wenn man 
die Matrix X durch die Matrix AM für genügend großes 
(ganzes positives) M ersetzt, d.h. an Stelle des Glei- 
chungssystems = Xr + c das System = AMr-+ d 
mit d = VWM—1c + YM—2c+ ... +NXc + c iteriert. 

Abschließend sei noch bemerkt, daß die in III an- 
gestellten Betrachtungen über die Konvergenz p-ter 
Art des Iterationsverfahrens und das dafür aufge- 
stellte Konvergenzkriterium nicht auf lineare Glei- 
chungssysteme beschränkt sind, sondern sich auf be- 
liebige Gleichungssysteme: x» = @ (&,, %g, - ++, Zn) 
w=1,2,...n) übertragen lassen, bei denen die 
Funktionen a, (x,, &3, -.-, 2n) in der Umgebung der 
gesuchten Lösung mindestens einmal stetig diffe- 


renzierbare Funktionen von &,, £3, -. -, 2n sind. Die 


Matrix in (3) hat man lediglich durch die Funktional- 


dd, oa \ 
= O2, Oxn 
matrix: %=| : : | im Lösungspunkt tr, 
dan da, 
on 


0% 
1 = 
zu ersetzen. Liegt der Näherungswert r, genügend 


"nahe bei r,, so wird man X, durch W,, d.h. die im 


Punkte r, gebildete Funktionalmatrix ersetzen kön- 
nen. Man hat dann das obige Abschnittsdetermi- 

. . 2 a 
nanten-Kriterium auf die Matrizen E — (UP) A an- 
zuwenden, um die spezielle Art der Konvergenz zu 
erkennen. 


Frankfurt/Main. E. Fehlberg. 


s) Die Frage, wie man Konvergenz des Iterationsverfahrens 
überhaupt erreichen kann, lassen wir hier beiseite, da sie in der 
Literatur zur Genüge behandelt ist (vgl. z.B. die Beschreibung 
des Rungeschen Verfahrensin Runge-König: Numerisches Rech- 
nen. Berlin 1924. p.1837). 
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Zur genäherten Auflösung algebraischer 
Gleichungen'!). 


1. Eine Näherungsmethode zur analytischen 
Fortsetzung 


Wendet man das Newtonsche Verfahren oder die 
regula falsi zur Bestimmung komplexer Wurzeln 
algebraischer Gleichungen an, bereitet das Auffinden 
einer ersten Näherung Schwierigkeiten. Im folgenden 
wird eine Methode beschrieben, die diese graphisch 
mittels genäherter analytischer Fortsetzung findet. 


Durch w=f(z) werde eine konforme Abbildung der 
z-Ebene auf die w-Ebene vermittelt. Bekannt seien 
die Bildpunkte der Punkte 9m =r: em + mp) 
(m=0, +1, +#23,..,0<p<n). Gesucht sind 
die Bildpunkte der Punkte „m = r: !M+mArnp 
(= +1, +2, ...). Ohne Einschränkung der Allge- 
meinheit kann man r=1, 9,=0 und f(1)=1 voraus- 
setzen. Es handelt sich also darum, das Bild eines 
Quadratnetzes der Z-Ebene (& =1og 2) zu finden. ‚Es 
werden die Bezeichnungen a = f(e'?), b=f(eT"?), 
c=f(e?) und d=f(e”?) eingeführt. Die erste Nähe- 
rung c, bzw. d, für c bzw. d erhält man durch lineare 
Fortsetzung mit Benutzung von « und b. Aus Bild 1 
ergibt sich 


a—b 
2i 


ee 
ee N -- 


: -tga& 
sin 


(1). 


cosp—e ? 


2 
a a 
sin p 


ern udn > di 


Die Konstruktion von c, und d, kann man sich er- 
leichtern, wenn man ein aus durchsichtigem Koor- 
dinatenpapier geschnittenes Viereck benutzt, das aus 
zwei gleichschenkligen Dreiecken mit den Basis- 
winkeln x und ß zusammengesetzt ist. 


z-Ebene 


Bild 1 


Die zweiten Näherungen c, und d, erhält man durch 
Fortsetzung eines Polynoms w=&,+ &% + Ag2?, 


1) Diese Arbeit stellt einen Teilder von der Technischen Hoch- 
schule Dresden 1951 genehmigten Diplomarbeit des Verfassers dar, 
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das mit f(z)in den Punkten z—= e'?,z=1lunde=e"'P ; 
übereinstimmt, Führt man: A = 5 5 BER 1 und > 
F 
Er a: ß ein, erhält man 3 
A-cosp + B-sino ä 

&ı — ar 
2 sin?5 ; 

(2). 

A cos 5+ B sin | 
=— — 
B F 


Bild 2 . ER, = 


Wegen fll)=1istx, +0, +0,=1,also 
a-1l= a, (1) + 5, (et). 


„e? —cosp 


. ’ 
sin gp 


4—l=4A+B 
erhält man : ze 
g— = elO AND 
29 _ 2. BR 
wobei g(p) FE 289° und 
4 sin®5- 005 | os 
ö(p)=A cos 5 + Bsin gist, Die Konst 
in Bild2 durchgeführt. Man wird dab 
Fällen A cos? mit genügender C 


ersetzen können. Entsprechend ist 


Da 
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unabhängig von f(z) ist, kann man sich die Kon- 
struktion von ö ersparen, wenn d bekannt ist. Ersetzt 
man nämlich in (4) d, durch d, ergibt sich 


Ne Pla a). ide) 


als Näherung für c. Ihr Fehler ist ebenfalls klein von 
dritter Ordnung. Entsprechend ist 


‚ Meute PioZe)...rlab). 


Die Fortsetzung wird demnach in der Weise durch- 
geführt, daß man ausgehend von den Bildpunkten der 
Punkte 4m =r’e ya} nach (3) die Bildpunkte 
der Punkte m =r IM +mm+p und daraus nach 
(5a) die Bildpunkte der Punkte z,,—= r: ed M +mp) +29 
konstruiert. So fährt man fort. Entsprechend geht 
man bei der Konstruktion der Bildpunkte der Punkte 
zn, mr MtmM—RP yor, In der folgenden 
Tabelle 'sind a, ß, g(9), g(—P), e®? und €?P für 
einige Werte von p zusammengestellt. Bei Berech- 
nung der Winkel x und $ kann man die aus (1) sich 
ergebende Beziehung « —ß = benutzen. 


) 
5° | 47°28’ | 42°28° | 2,19 | 1,84 | 1,19 |. 0,84 
10° | 49°517 | 39°51° | 2,40 | 1,69 | 1,42°| 0,71 
15° | 52°10’ | 37°10’ | 2,64 | 1,56 } 1,69 | 0,59 
20° | 54°25’ | 34°25’ | 2,91 | 1,45 | 2,01 | 0,50 
2571 56°357 | 31°35° | 3,22/| 1,35 | 2,39 | 0,42 
30° | 58°42’ | 28°42’ | 3,58 | 1,26 | 2,85 | 0,35 
221%° | 55°31’ | 33° 17 | 3,06 | 1,39 | 2,19 | 0,46 
4,81 | 0,21 


45° | 64°33’ | 19°33’ | 5,01 | 1,04 
B 


2. Anwendung zur Lösung algebraischer 
Gleichungen 


Will man jetzt die beschriebene Methode zur Lösung 
einer algebraischen Gleichung mit reellen Koeffi- 
zienten verwenden, wird man 2,=r:e'”, also den 
Ausgangspunkt für die Fortsetzung, nicht beliebig 
wählen. Es soll vorausgesetzt werden, daß die Glei- 
chung keine reellen Wurzeln hat. Schreibt man sie in 
der Form ; 

n—1 


I): + 2 AQmz”=0 (n gerade) . (6), 
_ 20 = 


va so hat f (2) mindestens ein reelles positives Minimum. 
) mit der Abseisse &. Sei also f(&$)=n und f(&) = 0 
’ (& reell, n positiv). Nähert man jetzt f(z) durch 


Ka—ntlf’Oe@-9°....M 


an, so sind die Nullstellen von h(z) grobe Näherungen 
für die Nullstellen von f(z2). Man wählt also 

20 = & +iV2n/f’’ (€) und bestimmt f(z,)entwederrech- 
nerisch mit Hilfe des erweiterten Hornerschen Sche- 
mas?) oder graphisch durch Vektoraddition. Die Kon- 
| struktion durch Vektoraddition läßt sich für die Bild- 

punkte eines Kreises |z| = const leicht mit Hilfe fest 
 eingestellter Zirkel durchführen. Wenn |f(z,)| schon 
genügend klein ist, kann man sofort mitdem Newton- 
‚schen Verfahren oder der regula falsi weiterrechnen. 
Andernfalls setzt man f(z) unter Verwendung von 
ildpunkten des Kreises |z|= |z,| wie oben be- 
eben in Richtung auf den Nullpunkt fort. In 
meisten Fällen werden ein oder zwei Schritte 

en, um eine brauchbare Ausgangsnäherung für 
Rechnung zu finden, wenn man die Schritt- 
t zu klein wählt. Manchmal kann man 
ortsetzung ganz ersparen, wenn man be- 


z, Z.angew. Math. Mech. 20 (1940), $. 235—236. 
d 


achtet, daß man sich in der Umgebung eines 
stationären Punktes von f(z) befindet. Das Bild eines 
2=£ einschließenden Kreises K der z-Ebene hat dann 
eine w=n einschließende Schlinge. In Bild 3 ist der 
Fall gezeichnet, daß f’(z) in z=£ eine einfache Null- 


stelle hat. Die Bildkurven der Kreise K, und K, sind 


mit X] und K, bezeichnet. Aus der Lage von f(z,) 
zum Punkt w=0 kann man mit Hilfe der in Bild 3 
dargestellten typischen Gestalt der Schlinge leicht eine 
bessere Näherung z, schätzen. 


Y z-Fbene 
a K2 
w-Ebene 
K; 
ni 
: An) 
Bild 3 


Hat die vorgelegte Gleichung komplexe Koeffi- 
zienten, kann man nicht in der beschriebenen Weise 
verfahren. Es lassen sich aber auch hier leicht obere 
und untere Schranken für die Absolutbeträge der 
Wurzeln ergeben. Ist 


n 


iz) = am m = 0 


m=0 


die Gleichung, so kann z=t keine Wurzel sein, wenn 


n—1 
| 77 ie | > | >: dm jun 


m=0 


gilt. Wegen 


n—1 n—1 Fr n—1 
2 mir =, [amt |sıM P> | am | 
m=0 m=0 De) 
mit 
allerlei 
| 
ist das sicher erfüllt, wenn 
Ss | | 
Am 
een 3 
Pr M | an | 


ist. © bzw. ”/C ist also die gesuchte obere Schranke. 
Eine untere Schranke ist das Reziproke einer oberen 
Schranke für die Wurzeln der Gleichung 


n m 

=) ad 

> dn-mr —=0, 
m=0 


an 


BR e En 
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deren Wurzeln ja die Reziproken der Wurzeln der 
ursprünglichen Gleichung sind. Mit 


[ao | 
3 | am | 
m—1 


ist also entweder D(für D<1) oder */D (für D>1) 
die gesuchte untere Schranke. Andere Schranken sind 
bei L. Bieberbach, Vorlesungen über Algebra 
(Leipzig und Berlin or S. Gr) angegeben: 


Bee Sys 
fz)=0, wobei N > lam| ( ma en)srAls 
Radius des a wird man etwa das 
geometrische Mittel aus oberer und unterer Schranke 
nehmen, da man dann zur Fortsetzung von f(z) über 
das Kreisringgebiet nach beiden Seiten etwa gleichviel 
Schritte braucht. Man kann die Ausgangspunkte 
wieder durch Vektoraddition gewinnen. Zu Beginn 
wird man die Maschenweite » nicht zu klein wählen. 

Verwendet man zur weiteren Verbessernng der 
Wurzeln die regula falsi, kann man ebenfalls graphisch 
vorgehen): Sind z, und z, zwei Näherungen für eine 
Wurzel, bestimmt man zunächst rechnerisch w, =f(z}) 
und w;=f(2,) und überträgt dann das Dreieck mit 
den Eckpunkten 0, w, und w, ähnlich in die z-Ebene. 
Die w=0 entsprechende Ecke ist dann die neue 
Näherung z;. 


D= 


3. Beispiel 
Als Beispiel soll 


+ (14 4)22 42 (— 2,54 081)24+3 + 2i=0 


behandelt werden. Aus |a,| = 1, |a,| = 1,41, |a,| = 2,62 
und |a,| = 3,61 findet man 7,6 "und 0,7 als Schranken 
für die Absolutbeträge der Wurzeln nach der ersten 
- Methode; 4,6 und 0,5 nach, der zweiten Methode. Als 


V 


Radius des Ausgangskreises ist r—=2 gewählt. Es ist 
dann asr2—8, la,r?| — 5,66, arg (4,7?) — — 161°25° 
laır| = 5,24 und arg (a,r) = 45°. Man erkennt aus 
einer überschlägigen Zeichnung, daß das Bild des 
Kreises „1? | =2 sich für Argumente in der Nähe 
von 180° dem Nullpunkt nähert. In Bild 4 sind des- 
halb die durch Vektoraddition gewonnenen Bildpunkte 
der Punkte ,—2e*/, „ —2e'” und 2, = 2e"i”I® 
eingezeichnet. 

Zwischen w, und w, hat die Bildkurve eine auf den 
Nullpunktzielende Ausbuchtung. Es ist deswegen noch 
der Punkt w,,, als Bild des Punktes 2,,, = Ze?! 12 
eingeschaltet. Mit Hilfe der Punkte w,, w;,, und w, 
ist dann das Bild des Punktes ?,,, = Dar ANNE 


®) Fr, A. Willers, Methoden der praktischen Analysis. Berlin 
1950, 8.256. 
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— 2,6. 0? i®12 näherungsweise wie oben beschrieben 
konstruiert. Aus seiner Lage zum Nullpunkt kann 
man ohne weiteres entnehmen, daß 2, = rye'?° mit 
r,—=2,5 und 9, = 190° eine brauchbare Ausgangs- 


näherung ist (2, = — 2,46 — 0,43i). Durch Ver- 
besserung dieser Näherung findet man den genaueren 
Wert 2= — 2,4592 — 0,5542 i. 


Die angegebene Näherungsmethodezur analytischen 
Fortsetzung läßt sich natürlich auch bei der Bestim- 
mung von Nullstellen anderer analytischer Funktionen 
verwenden. 


Dresden. K.-H. Bachmann, 
Zur Stabilität des Differenzenverfahrens bei 
der Stabschwingungsgleichung. 

Bei der Gleichung für freie Querschwingungen ge- 
rader Stäbe konstanten Querschnitts im einfachsten 
Fall!) 

ou 02 u 


Ort Oy? 

(Auslenkung u aus der Ruhelage an der Stelle x zur 

Zeit y; K gegebene Konstante) werde ein rechteckiges 

Gitter mit den 'Maschenweiten A und Z und den 
y=yr—kl 


Gitterpunkten 
£ } (2) 
Gk=0, +1, +2...) 


v==jh; 

zugrundegelegt, und U;; sei ein Naherungeyert für 
u(xj, yr). Beim gewöhnlichen Differenzenverfahren 
ersetzt man (1) durch die Differenzengleichung 
Uj+2,3 — 4Uj41,% + 605,5 — 4U;-1,6 + Uj2,% 
h*ı 
Uj,r+41 — 295,6 + Oj,r-1 
+K a 0 


Diese Formel wird verwendet, um den Wert U5,2+1 
auf der Zeile y=(k-+-1)! auszurechnen, wenn die 
Werte auf den beiden vorangehenden Zeilen yi=kl 
und y = (k — 1)l bereits bestimmt sind. Mit der Ab- 


a 


(3). 


kürzung 
== RE (4) A 
Eu ar 
lautet (3) 
U;,8+1 = — U;,x-1 — 2U;42,% + 42041, 


"0. 
+(2— 62)U;,% + 42U;-ı, k—2Uj-, k. 
Wählt man 2 =1, so entsteht 
D;,8+1 = — Uj,8-1— Uj42,6+ 4 0541, 
— AU; RK +4U;—1,R —-U;-,r. 

Diese Formel ist für die praktische Rechnung un- 
brauchbar, weil kleine Fehler in den U-Werten, z. B. 
schon die "Abrundungstehler, sich so ungünstig fort- 
pflanzen, daß die dadurch hervorgerufenen Schwan- 
kungen die Rechnung ganz überdecken; das Verfahren 
ist „instabil“. Ein Zahlenbeispiel in dieser Zeitschrift 
bestätigt diese Erscheinung?). Herr K. Oswa- 
titsch?°) hat nun dasselbe Beispiel mit kleinerer 
Zeitmaschenweite durchgerechnet (und zwar mit dem 
Wert z=1/6) und dabei stabiles Verhalten beim 
Differenzenverfahren festgestellt. Es soll daher herr 
untersucht werden, bis zu welchem (günstigsten) 2- 
‚Wert das Differenzenverfahren sich stabil verhält 
Dies läßt sich nach der Theorie von J. v. Neu- 
mann) durchführen, hier aber auch nach 
elementareren Methode 5). Dabei legt man au 


2) Vgl. etwa Handbuch der nun Bd. VI, 
elastischen Körper. Berlin 1928 8.3 EM 

®2) L. Collatz: Über das Ditterenzenverfahre n b 
wertproblemen partieller Ditferentlalgleichung en 
Mech. 16 (1936), S. 239— 247, J 

») Die Rechnung von Herrn K. Os watits 
zu der ne Betrachtung. 

*) R. P, Eddy: Naval Ordnance Tabo 
Ba. x (1949), S. 232. Rn 

O0. Brien, M. Hyman, 8. ] 

N NEO Solutions of Partial Differe 
of Mathematics and Physics 29 (1951) 


Z. angew. Math. Mech. 
Bd. 31 Nr. 11/12 Nov./Dez. 1951 


Kleine Mitteilungen 


393 


eine spezielle Art von Anfangsrandwertaufgaben zu- 
grunde. So möge es sich hier etwa um einen beid- 
seitig gelenkig gelagerten Träger endlicher Länge mit 
den Anfangs- und Randbedingungen handeln: 


2,0), | 
Ge fr 0<e<a, 
92,0) n er 


wobei f,(z) und f,(x) gegebene Funktionen sind. Das 
Intervall< 0, a > werde in N gleiche Teile der Länge 


(6), 


h= 7 geteilt. Auf zwei benachbarten Zeilen, etwa auf 


den zuk=p,k=p-+.1 gehörigen Zeilen, seien die 
U, 2-Werte mit gewissen Fehlern &;,» bzw. &j,p-+1 be- 
haftet, und man fragt, ob die dadurch beim Rechnen 
nach der Gl. (5) bewirkte Fehlerverteilung ej, x für 
k — oo beschränkt bleibt oder nicht. 

Nun hat Gl.(5) die Partikularlösung 


wenn (wie man beim Einsetzen sofort nachrechnet) 
zwischen & und 7 der Zusammenhang besteht 
22(cs&— 1? =1— Con 418) 


Wählt man für & die Werte E= &,= T und sind 


+n, nach (8) zugehörige Werte von n, so ist auch 
NA iv = = 

2,2 Zain 7 (A, 42,6). 
y= 


UNSER 


eine Lösung von (5), welche den (6) entsprechenden 
Differenzen-Randbedingungen 
Vene 9,n  2Vo,n tt 91,E = On4l,% 
20 gr 01,0 (k=0,1,2,...) 


genügt und für k=p undk = p-+ 1 die Werte &,,p 


bzw. &,p+1 annimmt, sofern man die noch freien. 


. Konstanten A» und B, aus dieser Forderung be- 
stimmt. Es ist stets möglich, die A, und B, auf diese 
Weise zu ermitteln; setzt man nämlich 

N—l 


. vn Ni... je 
# we N’ 4,9417 2 Pyein Na 
: y= R v- 


wobei sich die a,, f, nach den bekannten Formeln der 
harmonischen Analyse ergeben, so erhält man für die 
4A, und B, die Gleichungen 
A,+ B,=0,; A," L Betr —=Bß,, 
, ) aus denen sich A, und B, eindeutig berechnen, da 
7,0 ist (mit cos &=#1, 20 ist nach (8) auch 
Cfin=+1 also Sinn--0 und damit die Deter- 
AN minante des Gleichungssystems + 0). Somit ist in 
F (9) 9,2 = &;,, und e, „bleibtfür k — oo bei beliebiger 
Fehlerverteilung, d.h. bei beliebigen A,, B,, genau 


E% dann beschränkt, wenn 
, Bi e. für 


a je” | = ı velsaNnmı 
ist. Die Gleichung s=1 — Coj n hat bei gegebenem 
reellen sim Streifen |Jmnj <r Lösungen n, die für 
© von O bis 2 anwachsendes s auf der imaginären Achse 
liegen und von O0 bis -+in wandern, und die für 
8>2 auf den Geraden Jmn= +” liegen. Für 
0<s<s2 ist dann |e?|=1 und für s>2 wird 
je”) >1. Die Stabilitätsforderung verlangt daher 
ee... 232(c08&,— 1? <S2 


W=1,.,N-1 (10). 


rhält eine einfache, für alle N gültige 
enn man den Nenner durch seinen Größt- 
ap | 4 


wert 4 ersetzt; die Stabilität ist daher gesichert für 


za | 
oder I<YRZ RER 


AN 
A| 


Die Differenzengleichung (5) nimmt dann für = = 
die Form an 


1 
Gm = ut Zn Yen 


(12). 
+49,,1, 520,424 40, 1,2 a 


Man darf also auch hier die Zeitmaschenweite nicht 
zu groß wählen. 
Hannover. L. Collatz. 


Numerische Werte des Fehlerintegrals für 
Ynz. 

‚ Zu den zahlreichen Anwendungen des Ga ußschen 
Fehlerintegrais \ 


En 
F(&) = En () 
Er Yn CR ve un r ı 

2 o 

gehört seit einigen Jahrzehnten auch die numerische 
Berechnung gewisser Werte von Zetafunktionen, dieins- 
besondere in der Kristallphysik (Berechnung der Gitter- 
energie periodischer Ionenanordnungen) gebrausht 
werden. Hierbei ist man zwar nicht an eindeutig be- 
stimmte Werte des Arguments gebunden, jedoch bietet 


die Verwendung der Werte r—= nz, wo n eine natür- 
liche Zahl bedeutet, mehrfache Vorteile. In neuerer 
Zeit wird das Fehlerintegral hierbei mit größerer Ge- 
nauigkeit benötigt, so daß die allgemeiner zugänglichen 
Tafeln dieser Funktion nicht ausreichen. 

Tabelle 1 


n F (Yan) 


0,98781 11178 15197 10 
0,99960 72494 11718 
0,99998 58555 97311 
0,99999 94648 35338 
0,99999 99791 73448 
0,99999 99991 74578 
0,99999 99999 66864 
0,99999 99999 98657 
0,99999 99999 99945 
0,99999 99999 99998 

zı 1,00000 00000 00000 
Für 2s<r<s7 beträgt der Fehler möglicherweise bis zu + 1 Ein- 
heit der letzten Stelle. 

Daher sollen die Werte im Folgenden mit etwa 
14 Dezimalstellen mitgeteilt werden: siehe Tabelle 1, 
deren Entstehung kurz erläutert werde: 

1. Aus (1) folgt durch Integration der Potenzreihe 
für die Exponentialfunktion die bekannte in der gan- 
zen Ebene konvergente Potenzreihe des Fehlerintegrals 

2 x? a5 Fo 
Pa= le ta - 3 en (2), 
die sich zur genaueren Berechnung im allgemeinen nur 
bei kleineren |x| eignet. Für «= /z erhält man: 


1 PR: E17 2 3 
zla=1-3t35- 307 + 0 9. 


Da die Werte von nk/k! für positive ganzzahlige k<< 16 
mit mindestens 18 Dezimalstellen tabuliert sind!) und 
sie für einige weitere k mit normalen Hilfsmitteln be- 
rechnet werden können, ergab sich diese Potenzreihe 
relativ einfach zu 


F (Yr) = 0,98781 11178 15197 1028. 
+4 


nt 
SO ID RUND - 


2. Die übrigen Werte von F(Ynx) für ganzzahlige 
n wurden durch Interpolation aus vorhandenen Tafeln 


1) K, Hayashi, Tafeln der Besselschen, Theta-, Kugel- 
und anderer Funktionen, Berlin 1930, Tafel IX, 8.3 und 5, 
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berechnet?). Für n = 2,3 und 4 wurde kubisch, für 
n—=5 und 6 quadratisch interpoliert. 

3. Für n= 17,8 und 9 konnte linear interpoliert 
werden. Dabei wurden die für n—= 8 und 9 vorhande- 
nen Schwankungen der ersten Differenz im Betrag von 
je einer Einheit der letzten Stelle (die Differenzen be- 
trugen 14,13,14 für n= 8 und 1,1,0 für n = 9), die 
offenbar eine Folge von Auf- oder Abrundungen 
darstellen, dazu verwendet, um diese Unsicher- 
heit der letzten Stelle nicht nur auszugleichen, 
sondern die Fehlergrenze unter den üblichen 
Wert von +0,5 der letzten Stelle herabzu- 
drücken. Voraussetzung hierbei ist, daß bei 
Erstellung der Tafeln die Auf- und Abrundun- 
gen exakt erfolgt sind. 


0) 


Tabelle 2 
© F(«) 
5,316 0,99999 99999 99944 
DT 0,99999 99999 99945 
5,318 0,99999 99999 99946 
5,319 0,99999 99999 99946 


Wegen 2, = Y9r= 5,317 362 geht man aus von 
dem in Tabelle2 wiedergegebenen Auszug aus der 
Tafel?). In Bild 1 sind die Werte der letzten Stellen 
eingetragen, wobei die Intervalle der wahren Werte 


946 


...945 


- _} 0.000362 


Ih 


7% 8 9 


5,315 6 


= i Bild 1 


von der Höhe je einer Einheit der 15. Stelle durch einen 


stärkeren Strich angedeutet sind. (Strecken 17 P; s 
j= 67,89). Da die Funktion in diesem Intervall 
praktisch geradlinig verläuft, muß sie durch eine alle 
4 Strecken P, 2 schneidende Gerade darstellbar sein. 


Aus der speziellen Lage dieser 4 Strecken zueinander 
ersieht man, daß der gesuchte Wert innerhalb des 


Fünfecks P5 P,P,Ps,P? liegen muß (P” sei die Mitte 


zwischen P, und P,), für 5,317 < x, < 5,318 also in 
dem Dreieck ER P,, so daß F (x,) in einem Intervall 
von der Länge Q-Q+= = (2 — 0,362) = 0,546 liegen 
muß und dem unteren Intervallende Q- mit 

Fix.) >(... 944,5 + 0,362) 10-15, 
"®) Tables of probability functions, I u. II, Condusted under 


the sponsorship of the National Bureau of Standards, New York 
1941, (15 Dezimalstellen mit Intervallen von 0,001), 


5,320 


also F(y9r) = 0,90?) 4513. 
+28 


Dabei bedeutet die Ziffer in Klammern hinter der auf _ 
das Komma folgenden 9 die Anzahl der dort vorhan- 
denen Neunen. 


Für n=8 vgl. das — schematisierte — Bild 2. 
Es ergibt sich F (8 r) = 0,901) 865 715. 
+19 


5,011 2 3 % 4 5 


Bild 2 


4. Anwendungen. Die angegebenen Werte ge- 
statten nunmehr eine genauere Berechnung gewisser 
Werte Epsteinscher Zetafunktionen, insbesondere 


z|o) (a=—4+ an Br Me] (4), 


wo r,(n) die Anzahl der Zerlegungen der natürlichen 
Zahl n in die Summe von 2 Quadraten bedeutet, also 
Fs(n)— Olten —19,6,7,1112; 

= 4 für n= 1,2,4, 8, 9; 

-8fürn=510. . 
In Verschärfung einer früheren Berechnung?) ergibt 
sich 


Z 0 | (0= —3,000%6 49200 020 . . (8) 


und 


00 
BI 
22 


ei 00 “. 
! x .—= — 1,61554 26267 128 
Br 
Letzteres hat die Bedeutung der alter- 
nierenden Summe der reziproken Ab- 
stände eines Gitterpunktes in einem quadratischen 
Gitter von allen übrigen nach wachsender + 
Entfernung geordnet. Dieser Wert wird als x 
Madelung- Konstante der alternierenden Netz- 
ebene mit großer Genauigkeit benötigt, um einen PR 
asymptotischen Ausdruck für die Gitterenergie nd 
licher, quadratischer Stücke einer solehen Netzebene EN 
zu erhalten ®), j sn 
Herr cand. math. H. Eitel, Darmstadt, hat sich 
freundlicherweiseder Mühe unterzogen, mir mit liebens- 
würdiger Zustimmung von Herrn A. Waltheraus 
einemim dortigen Institut vorhandenen Tabellenwerk?), 
die von mir benötigten Zahlenwerte mitzuteilen. Herr 
cand. math.H. Thiele, Berlin, hat mir bei Durch- 


Z.angew. Math. Mech, 


Buchbesprechungen 


395 


Bd. 31 Nr. 11/12 Nov./Dez. 1951 
7 


BUCHBESPRECHUNGEN 


L.Collatz (Prof. a. d. Techn. Hochschule Han- 
nover,, Numerische Behandlung von 
Differentialgleichungen. (Die Grund- 
lehren der Mathematischen Wissenschaften in Einzel- 
darstellungen. Bd. LX.) XIII + 458 S. m. 110 Abb. 
Berlin-Göttingen-Heidelberg 1951. Springer-Verlag. 
Preis geh. 45,— DM, geb. 48,— DM. 

Verf. gibt in dem vorliegenden Buch eine um- 
fassende Übersicht über die Methoden der näherungs- 
weisen Integration von Differentialgleichungen, wobei 
er insbesondere die neuerdings auf den verschiedensten 
Gebieten entwickelten Verfahren einheitlich dar- 
stellt. Er berücksichtigt in weitem Umfang insbe- 
sondere in den angelsächsischen Ländern in den letz- 
ten Jahren entwickelte Methoden, kann aber auch 
vielfach auf eigene Arbeiten zurückgreifen. Im ersten 
Kapitel behandelt er Anfangswertaufgaben bei ge- 
wöhnlichen Differentialgleichungen. Er gibt zunächst 
einfachere Summationsverfahren für Differential- 
Gleichungen erster Ordnung,dann auch fürGleichungen 
höherer Ordnung das. Verfahren von Runge- 
Kutta und die Verfahren der tabellarischen Inte- 
gration. Kapitel II bringt Randwertaufgaben bei 
gewöhnlichen Differentialgleichungen, und zwar wer- 
den zunächst die zum großen Teil vom Verfasser ent- 
wickelten verschiedenen Differenzenverfahren, dann 
die auf den direkten Methoden der Variationsrechnung 
beruhenden Verfahren und einige andere behandelt. 
Kapitel III beschäftigt sich mit Anfangs- und An- 
fangsrandwertaufgaben bei partiellen Differential- 
gleichungen, für die wieder verschiedene Differenzen- 
verfahren in Frage kommen. Hier wird auch auf die 
Charakteristikentheorie und ihre Verwendung zur 
numerischen Integration eingegangen. Kapitel IV 
bringt die Differenzenverfahren, das Ritzsche, das 
Trefftzsche Verfahren usw. zur Lösung von Randwert- 
problemen bei partiellen Differentialgleichungen und 
Kapitel V endlich Methoden zur genäherten Lösung 
von Integral- und Integro-Differentialgleichungen. 
Hier geht Verfasser auch kurz auf Funktional- 
gleichungen ein. . - 


Das Buch zeichnet sich durch umfangreiche Unter- 
suchungen über die Konvergenz der behandelten Ver- 
fahren aus. An vielen Stellen findet man sowohl re- 
kursive wie independente Fehlerabschätzungen. Die 
reiche Erfahrung des Verfassers in der Anwendung 
der in dem Buch dargebotenen Methoden findet ihren 
Niederschlag in zahlreichen Hinweisen und Winken 
für die praktische Durchführung, in Angaben über den 
Arbeitsaufwand, den die verschiedenen Verfahren 
zur Erreichung der gleichen Genauigkeit erfordern 
und vor allem in den zahlreichen instruktiven Bei- 
spielen, die zu jeder Methode durchgerechnet werden 
und die sich nochmals, den in dem betreffenden Ab- 
schnitt behandelten Stoff zusammenfassend, am 
Schluß eines jeden Kapitels nebst Lösungen finden. 

Das Buch schmückt ein Bild von Carl Runge. Mit 
Recht! Sein Inhalt baut vielfach auf Arbeiten Run- 
ges auf. Er hätte sicher Freude an dem ausgezeichne- 
ten Buch gehabt. 


Dresden. Willers. 


Dr.-Ing. Karl Klotter (o. Prof. a. d. Techn. Hoch- 
schule Karlsruhe), TechnischeSchwingungs- 
lehre. 1. Band: Einfache Schwinger 
und Schwingungsmeßgeräte. Zweite um- 
_ gearbeitete und ergänzte Auflage. 399 8. mit 360 Abb. 


Mm Berlin-Göttingen-Heidelberg 1951. Springer-Verlag. 


ste Auflage des 1. Bandes des Klotterschen 
rschien bereits im Jahre 1938 (s. Besprechung 
Math. Mech. Bd.18, 8.255). Verf. plant 

gabe eines dreibändigen Werkes, sah sich 
aus verlagstechnischen Gründen veranlaßt, 


noch vor Herausgabe des 2. Bandes eine Neuauflage 
des 1. Bandes zu bringen. Diese ist wesentlich 
umfangreicher als die erste Auflage ausgefallen. Ins- 
besondere wurde das Büchlein ‚‚Messung mechani- 
scher Schwingungen, Dynamik der Schwingungs- 
meßgeräte“, Berlin 1943, Springer-Verlag, das Verf. in- 
zwischen als eine Art Ergänzung des 1. Bandes ver- 
öffentlicht hatte, in ergänzter Form aufgenommen. 
Weiterhin ist die Aufnahme eines Kapitels über 
rheolineare Schwinger bemerkenswert, das sich mit 
solchen Schwingern befaßt, deren Differential- 
gleichungen veränderliche Koeffizienten aufweisen 
(Beispiele: Pendel im Schwerefeld, Pendel und 
elastischer Schwinger mit erschütterter Drehachse, 
Druckstab mit pulsierenden Längskräften, Torsions- 
schwingungen von Kurbelwellen u.a.); die fort- 
schreitende Technik hat insbesondere auf diesem 
Gebiete eine Fülle von Anregungen gebracht, die zu 
interessanten Stabilitätsbetrachtungen Veranlassung 
gegeben haben. Die Fachwelt wird die Aufnahme 
dieses Kapitels besonders begrüßen. Im übrigen war 
der Verf. bemüht, hinsichtlich Vollständigkeit und 
ae verschiedene Verbesserungen anzubringen, 
so daß nunmehr eine recht geschlossene Darstellung 
des Schwingers mit einem Freiheitsgrad vorliegt und 
die Streichung des Wortes ‚Einführung‘ im ur- 
sprünglichen Titel des Buches gerechtfertigt er- 
scheint. 


Dresden, H.Neuber. 


Edward R. van Driest, Die linearisierte Theorie 
der dreidimensionalen kompressiblen Unter- 
schallströmungunddieexperimentelle Unter- 
suchung von Rotationskörpern in einem ge- 
schlossenen Windkanal. (Mitt. Institut f. Aero- 
dynamik an der E. T. H. Zürich, herausgeg. von 
Prof. Dr. J. Ackeret, Nr. 16) 30 S. Zürich 1949. 
Verlag Leemann. Preis: 7,— Fr. 


Die Prandtl-Glauertsche Regel gestattet recht 
bequem und mit guter Annäherungdie ebene Strömung 
um schlanke Körper im Unterschallbereich durch 
zugeordnete inkompressible Strömungen wiederzu- 
geben. In der vorliegenden Abhandlung wird diese 
Prandtl-Glauertsche Regel auf rotationssymme- 
trische Strömungen verallgemeinert, wobei in die 
Umrechnungsvorschrift, wie schon früher durch 
B. Göthert festgestellt wurde, neben der Machschen 
Zahl aber auch der Schlankheitsgrad des Körpers 
eingeht. Die theoretischen Ergebnisse werden für 
Rotationshalbkörper experimentell geprüft, unter 
Berücksichtigung der Kanalkorrekturen, die wesent- 
lich sind. Die gute Bestätigung der Theorie ist be- 
merkenswert. 


Göttingen. W. Tollmien. 


Dr. sc. techn. Ambros P. Speise, Entwurf 
eines elektronischen Rechengerätes 
unter besonderer Berücksichtigung 
der Erfordernis eines minimalen Ma- 
terialaufwandes bei gegebener ma- 
thematischer Leistungsfähigkeit (Mit- 
teilungen aus dem Institut für angewandte Mathe- 
matik an der eidgenössischen Technischen Hoch- 
schule in Zürich, herausgegeben von Prof. Dr. 
E. Stiefel), IV+65S. in 4Abb., Zürich 1950. 


Es wird ein elektronisch arbeitender, programm- 
gesteuerter digitaler Rechenautomat beschrieben, 
der etwa 1000 Elektronenröhren und 300 elektro- 
magnetische Relais erfordert. Der Entwurf basiert auf 
den Erfahrungen, die Verfasser an amerikanischen 
Rechenautomaten gesammelt hat und berücksichtigt 
diesen gegenüber die wichtige Forderung nach mini- 
malem Materialaufwand. 


396 Eingegangene Bücher — Nachrichten 


: g' 
Bd.31 Nr.11/12 Nov./Dez. 1951 


Als Speicher für Zahlen und Befehle dient die 
bewährte Magnetophontrommel, deren Einzelspeicher 
mit einem Relais-System aufgerufen werden, Das 
elektronische Rechenwerk für die vier Grundrech- 
nungsarten ist durch die Verwendung des Dezimal- 
systems charakterisiert, die Ziffernstellen sind in 
Dualform (Tetraden) nach dem Aiken-System ver- 
schlüsselt und werden nacheinander verarbeitet. 
Die Zahlen- und Befehlsein- und -ausgabe erfolgt 
unmittelbar mit handelsüblichen Fernschreibappa- 


raten. Das Gerät soll zur Addition 0,8, zur Multi- 
plikation 30, zur Division 40 ms benötigen und gehört 
damit noch zur Klasse der schnellen Rechenauto- 
maten. Das Leistungsvolumen genügt bei 1200 
Zahlenspeichern mit einer Suchzeit von 16 ms für 
jedes mittelgroße mathematische Institut. 

Der Entwurf beschränkt sich auf die Diskussion 
des prinzipiellen Aufbaues der Maschine und bringt 
keine elektrischen Schaltbilder. 


f' 
Dresden. N.Joachim Lehmann. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten), 


Arthur Linder (Prof. a. d. Univ. Genf und der ETH. 
Zürich), Statistische Methoden für Na- 
turwissenschaftler, Mediziner und 
Ingenieure (Lehrbücher und Monographien aus 
dem Gebiete der exakten Wissenschaften. Mathema- 
tische Reihe Bd. 3). 2. erweiterte Auflage. 238 S. mit 
42 Abb. Basel 1951. Verlag Birkhäuser. Preis brosch. 
26 Fr., geb. 30 Fr. 


P. Finsler (o. Prof. a.d. Univ. Zürich), ÜberKur- 
venundFlächeninallgemeinen Räu- 
men. (Lehrbücher und Monographien aus dem Ge- 
biete der exakten Wissenschaften. Mathematische 
Reihe Bd.11.) XT + 160 S. Basel 1951. Verlag Birk- 
häuser. Preis brosch 12 Fr., geb. 14,80 Fr. 


R.von Mises (Prof. a. d. Harvard-Univ.), Wahr- 
scheinlichkeit, Statistik und Wahr- 
heit. Einführung in die neue Wahrscheinlichkeits- 
lehre und ihre Anwendung. 3. neubearbeitete Auflage. 
IX +278S. Wien 1951. Springer-Verlag. Preis 
brosch. 18,— DM. \ 


H.Beckert, Bemerkungen über dieVer- 
biegung hyperbolisch gekrümmter 
Flächenstücke. H.Sali,ÜberAbelsVer- 
allgemeinerung der binomischenFor- 
m eln. Berichte über die Verhandlungen der Sächsi- 
schen Akademie der Wissenschaften zu Leipzig. Mathe- 
matisch-naturwissenschaftliche Klasse. Band 98 Heft4. 
228. Berlin 1951. Akademie-Verlag. 


Dr. Kurt Reidemeister (o. Prof. a.d. Univ. Marburg), 
Einführung in die kombinatorische Topo- 
logie (Die Wissenschaft, herausgeg. von Prof. Dr. 
Wilhelm Westphal, Band 86). Unveränderter Nach- 
druck. XIT-+ 2098. Braunschweig 1951. Friedrich 
Vieweg & Sohn. Preis geb. 12,— DM. 


Dr. A. Ostrowski (o. Prof. a.d. Univ. Basel), Vor- 
lesungen über Differential- und Integral- 
rechnung. Band 2: Differentialrechnung auf dem 
Gebiete mehrerer Variablen. (Lehrbücher und Mono- 
graphien aus dem Gebiet der exakten Wissenschaften 
mathem. Reihe Bd. 5.) 482 S. mit 55 Abb. Basel 1951. 
Verlag Birkhäuser. Preis brosch. 63,— Schw. Fr., 
geb. 67,— Schw. Fr. 


_ R. Sauer (Prof. a. d. Techn. Hochschule München). 
Ecoulements des Fluides Compres- 
sibles. XVI-+ 307 S. mit 150 Abb. Paris et Liege 
1951. Librairie Polytechnique Ch. Beranger. Preis 
geb. 3.900 fr. 


Dr.L. Vietoris (o. Prof. a. d. Univ. Innsbruck) und 
Dr. Gustav Lochs, Vorlesungen über Dif- 
ferential-undIntegralrechnung.4l6S8. 
mit 159 Abb. Innsbruck 1951. Universitätsverlag 
Wagner. Preis geb. 30,— DM. | 


Colloque de Topologie (Espaces fibres) Tenu & 


Liege du 5 an 8 juin 1950, 1378. Liege, Paris 1951, 
Georges Thone, Masson & Cie. Preis 175 francs 
belges, 1.225 francs francais. BR 


NACHRICHTEN 


Am 2. Oktober 1951 beging Prof. Dr. Ludwig 


Prandtl sein 50jähriges Professoren-Jubiläum. 


Z. angew. Math. Mech. . 


Schriftleitung und Verlag sprechen dem hervor- 


ragenden Wissenschaftler, dem erfolgreichen For- _ 7 de 
scher zu diesem seltenen Jubiläum die herzlichsten BEN 


Glückwünsche aus und danken ihm, der diese Zeit- 
schrift mitbegründete, für die Förderung, die er 
jederzeit der ZaMM zuteil werden ließ und für 
seine Mitarbeit. 


Die nächste Tagung der Gesellschaft für angew. München: Prof. Dr. Auma nn 


Mathematik und Mechanik wird in der Pfingstwoche Universität München berufen, 
1952 in Braunschweig stattfinden. h WE 
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